Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,135,181 C→T L172L (TTG→TTA)  unknown ← Type IV secretion system protein VirD4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,135,1810CT100.0% 211.7 / NA 78L172L (TTG→TTA) unknownType IV secretion system protein VirD4
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (46/32);  total (46/32)

GTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGTGCTCGGTCGCGCATATCAATACAATCCAATGGATTAAATCGGTGCGTGTTACGGTCAAACGGGGCAAAGCGATAGACTTTTTGCCCCAATATTTTTTGACGGTATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACGAATTGCGCACCGATAAAGCGTAAAAGTTTTTTGCCGAGCTTGCCGACGATAATGCCGTTTTCTTCTTCCATCAAGCCCATTTCCCGCGCTTCCGAG  >  1/1134937‑1135421
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
gTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAatat                                                                                                                                                                                                                                             >  1:411677/1‑250 (MQ=255)
gTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAata                                                                                                                                                                                                                                              >  1:772374/1‑249 (MQ=255)
gTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATAtt                                                                                                                                                                                                                                            >  1:897884/1‑251 (MQ=255)
gTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATAtt                                                                                                                                                                                                                                            >  1:300752/1‑251 (MQ=255)
 tATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAatat                                                                                                                                                                                                                                             <  1:858618/249‑1 (MQ=255)
 tATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATAttt                                                                                                                                                                                                                                           <  1:286276/251‑1 (MQ=255)
         cGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTgtcga                                                                                                                                                                                                                                   >  1:84567/1‑247 (MQ=255)
                gCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATAttt                                                                                                                                                                                                                              <  1:356409/249‑1 (MQ=255)
                 cTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATTTAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTAc                                                                                                                                                                                                                            <  1:119774/250‑1 (MQ=255)
                         gcgcATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATa                                                                                                                                                                                                                    >  1:796497/1‑250 (MQ=255)
                          cgcATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACg                                                                                                                                                                                                                  <  1:170613/251‑1 (MQ=255)
                              tatCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAg                                                                                                                                                                                                              <  1:876069/251‑1 (MQ=255)
                                       gTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACg                                                                                                                                                                                                                  >  1:864890/1‑238 (MQ=255)
                                       gTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGttt                                                                                                                                                                                                     >  1:806529/1‑251 (MQ=255)
                                       gTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGttt                                                                                                                                                                                                     >  1:241414/1‑251 (MQ=255)
                                       gTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGttt                                                                                                                                                                                                     >  1:779473/1‑251 (MQ=255)
                                       gTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGAGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATAttt                                                                                                                                                                                                                                           >  1:157528/1‑213 (MQ=255)
                                             aaGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATg                                                                                                                                                                                               <  1:795126/251‑1 (MQ=255)
                                                      cTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTcaca                                                                                                                                                                                      >  1:636583/1‑251 (MQ=255)
                                                        aCAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAAc                                                                                                                                                                                    <  1:463565/251‑1 (MQ=255)
                                                           aaaaGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACcac                                                                                                                                                                                 <  1:222724/251‑1 (MQ=255)
                                                                aTTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGAt                                                                                                                                                                             <  1:407013/250‑1 (MQ=255)
                                                                     gtgCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTc                                                                                                                                                                        >  1:601535/1‑250 (MQ=255)
                                                                                        gTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGaa                                                                                                                                                    >  1:658709/1‑251 (MQ=255)
                                                                                              aTCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAAc                                                                                                                                                                                    >  1:489360/1‑213 (MQ=255)
                                                                                              aTCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACaa                                                                                                                                              <  1:559981/251‑1 (MQ=255)
                                                                                                     tCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTa                                                                                                                                        <  1:246059/250‑1 (MQ=255)
                                                                                                              aGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTAc                                                                                                                               <  1:527708/250‑1 (MQ=255)
                                                                                                               gCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAACATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTAc                                                                                                                               >  1:454256/1‑249 (MQ=255)
                                                                                                                cGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCt                                                                                                                             <  1:664566/250‑1 (MQ=255)
                                                                                                                  gATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATg                                                                                                                                                                                               <  1:353291/182‑1 (MQ=255)
                                                                                                                  gATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTcac                                                                                                                                                                                       <  1:177614/190‑1 (MQ=255)
                                                                                                                  gATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATg                                                                                                                                                                   <  1:761916/210‑1 (MQ=255)
                                                                                                                  gATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTcc                                                                                                                           >  1:769445/1‑250 (MQ=255)
                                                                                                                  gATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTcc                                                                                                                           >  1:646119/1‑250 (MQ=255)
                                                                                                                  gATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTcc                                                                                                                           >  1:429828/1‑250 (MQ=255)
                                                                                                                   atatAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACg                                                                                                                                                                                                                  >  1:40005/1‑162 (MQ=255)
                                                                                                                   atatAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCAt                                                                                                                                                                                                           >  1:715876/1‑169 (MQ=255)
                                                                                                                   atatAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGAt                                                                                                                                                                             <  1:709998/199‑1 (MQ=255)
                                                                                                                    tataAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCg                                                                                                                        >  1:146287/1‑251 (MQ=255)
                                                                                                                        aGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATAtt                                                                                                                                                                                                                               >  1:63099/1‑144 (MQ=255)
                                                                                                                         gCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATg                                                                                                                                                                                               <  1:630866/175‑1 (MQ=255)
                                                                                                                         gCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCg                                                                                                                   >  1:16287/1‑251 (MQ=255)
                                                                                                                                                             cgcccGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTAcc                                                                                                                              >  1:688824/5‑204 (MQ=255)
                                                                                                                                                                       cgcATATCAATACAATCCAGCGGACTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagttt                                                                     <  1:392692/251‑13 (MQ=255)
                                                                                                                                                                         cATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATg                                                                                                                                                                                               >  1:680844/1‑127 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAAc                                                                                                                                                                                    <  1:55126/136‑1 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATGACAATCCCTACGc                                                                                                                                     >  1:452677/1‑183 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATGACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttatta                                                                 >  1:426273/1‑235 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATcc                                                                                                                                           >  1:868365/1‑177 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTAcc                                                                                                                              >  1:886792/1‑190 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTAcc                                                                                                                              >  1:888076/1‑190 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTAc                                                                                                                               <  1:534007/189‑1 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttatta                                                                 >  1:827872/1‑235 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttatta                                                                 >  1:341298/1‑235 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttatta                                                                 >  1:558210/1‑235 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttatta                                                                 >  1:345639/1‑235 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaac                                                                               >  1:590676/1‑235 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaac                                                                               <  1:791628/237‑3 (MQ=255)
                                                                                                                                                                           tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaac                                                                               <  1:231387/237‑3 (MQ=255)
                                                                                                                                                                              cAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttattaccg                                                              >  1:674293/1‑232 (MQ=255)
                                                                                                                                                                              cAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGCTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttattaccg                                                              <  1:673921/251‑20 (MQ=255)
                                                                                                                                                                                aTACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTAcc                                                                                                                              <  1:663433/185‑1 (MQ=255)
                                                                                                                                                                                aTACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttattaccgag                                                            >  1:225216/1‑230 (MQ=255)
                                                                                                                                                                                            ggATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGAt                                                                                                                                                                             <  1:892073/126‑1 (MQ=255)
                                                                                                                                                                                             gATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCt                                                                                                                                    >  1:96235/1‑166 (MQ=255)
                                                                                                                                                                                             gATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATa                                               >  1:281345/1‑251 (MQ=255)
                                                                                                                                                                                                         gCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAAt                                                                                                                                                   >  1:109102/1‑139 (MQ=255)
                                                                                                                                                                                                         gCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttattac                                                                <  1:90083/222‑18 (MQ=255)
                                                                                                                                                                                                         gCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaac                                                                               >  1:495215/1‑205 (MQ=255)
                                                                                                                                                                                                           gtattACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttattaccg                                                              >  1:551262/4‑203 (MQ=255)
                                                                                                                                                                                                            tattACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGAATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGAAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGttttttt                                  >  1:611429/3‑246 (MQ=255)
                                                                                                                                                                                                                                  gCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaac                                                                               <  1:123490/182‑3 (MQ=255)
                                                                                                                                                                                                                                   cGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCAttttttgtac         <  1:410405/251‑8 (MQ=255)
                                                                                                                                                                                                                                   cGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCAttttttgta          >  1:71111/1‑244 (MQ=255)
                                                                                                                                                                                                                                     atataccttttGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCAttttttgtactt       >  1:877799/8‑242 (MQ=255)
                                                                                                                                                                                                                                          ccttttGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCAttttttgtacttctgat  <  1:163853/249‑15 (MQ=255)
                                                                                                                                                                                                                                           cttttGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttattac                                                                <  1:469414/187‑18 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
GTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGTGCTCGGTCGCGCATATCAATACAATCCAATGGATTAAATCGGTGCGTGTTACGGTCAAACGGGGCAAAGCGATAGACTTTTTGCCCCAATATTTTTTGACGGTATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACGAATTGCGCACCGATAAAGCGTAAAAGTTTTTTGCCGAGCTTGCCGACGATAATGCCGTTTTCTTCTTCCATCAAGCCCATTTCCCGCGCTTCCGAG  >  1/1134937‑1135421

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑