Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,135,172 T→C K175K (AAA→AAG)  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,1720TC100.0% 221.9 / NA 77K175K (AAA→AAG) unknownType IV secretion system protein VirD4
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (47/30);  total (47/30)

GGCATAGGTCAGTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGTGCTCGGTCGCGCATATCAATACAATCCAATGGATTAAATCGGTGCGTGTTACGGTCAAACGGGGCAAAGCGATAGACTTTTTGCCCCAATATTTTTTGACGGTATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACGAATTGCGCACCGATAAAGCGTAAAAGTTTTTTGCCGAGCTTGCCGACGATAATGCCGTTTTCTTCTTCCATCAAGCCCATTTCCCGCGC  >  1/1134926‑1135414
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
ggCATAGGTCAGTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACctttt                                                                                                                                                                                                                                                >  1:162708/1‑251 (MQ=255)
ggCATAGGTCAGTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGCGCGTATTACGGTCAAACGGGGCAAAGCGATATACctttt                                                                                                                                                                                                                                                >  1:733331/1‑251 (MQ=255)
 gCATAGGTCAGTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTg                                                                                                                                                                                                                                               >  1:158817/1‑251 (MQ=255)
      ggTCAGTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGcccta                                                                                                                                                                                                                                          >  1:368887/1‑249 (MQ=255)
           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:779473/1‑251 (MQ=255)
                                                  gTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGttt                                                                                                                                                                                              >  1:241414/1‑251 (MQ=255)
                                                  gTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGttt                                                                                                                                                                                              >  1:806529/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:429828/1‑250 (MQ=255)
                                                                                                                             gATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTcc                                                                                                                    >  1:769445/1‑250 (MQ=255)
                                                                                                                             gATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGCGCCCGGTCGCGCATATCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTcc                                                                                                                    >  1:646119/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:888076/1‑190 (MQ=255)
                                                                                                                                                                                      tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTAcc                                                                                                                       >  1:886792/1‑190 (MQ=255)
                                                                                                                                                                                      tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTAc                                                                                                                        <  1:534007/189‑1 (MQ=255)
                                                                                                                                                                                      tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttatta                                                          >  1:827872/1‑235 (MQ=255)
                                                                                                                                                                                      tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttatta                                                          >  1:558210/1‑235 (MQ=255)
                                                                                                                                                                                      tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttatta                                                          >  1:341298/1‑235 (MQ=255)
                                                                                                                                                                                      tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaacggagcagtttatta                                                          >  1:345639/1‑235 (MQ=255)
                                                                                                                                                                                      tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaac                                                                        <  1:791628/237‑3 (MQ=255)
                                                                                                                                                                                      tatCAATACAATCCAGCGGATTAAACCGGTGCGTATTACGGTCAAACGGGGCAAAGCGATATACCTTTTGCCCTAATATTTTTTGTCGATATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATaaaac                                                                        >  1:590676/1‑235 (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)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
GGCATAGGTCAGTATCATATCGGTAAGGCTTTTTTGGCGCATATCCCACAGTAAAAAAGTGGTTGCTACAAAAAGATTCTGTGCCTGGTGCTGCCAAAAGTCTTTATCCGGGTCGCCGGTGAGCGGATATAAGCCCATAGCCTGATTTTTTAAATCGAGTTCGGCGCGTGCTCGGTCGCGCATATCAATACAATCCAATGGATTAAATCGGTGCGTGTTACGGTCAAACGGGGCAAAGCGATAGACTTTTTGCCCCAATATTTTTTGACGGTATTTACTGGTGATACGATAGCATTCGTTTTTAATGTCCTTCACAACCACAGATTGTTCCCATGCCAGCAAATTTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACGAATTGCGCACCGATAAAGCGTAAAAGTTTTTTGCCGAGCTTGCCGACGATAATGCCGTTTTCTTCTTCCATCAAGCCCATTTCCCGCGC  >  1/1134926‑1135414

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: ‑