Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,213 A→G 10.3% S746G (AGT→GGT)  mngB → alpha‑mannosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,2130AG10.3% 190.0 / 2.1 75S746G (AGT→GGT) mngBalpha‑mannosidase
Reads supporting (aligned to +/- strand):  ref base A (45/21);  new base G (2/6);  total (47/28)
Fisher's exact test for biased strand distribution p-value = 4.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGATTTTACGGCTGTTTAATCCGGCT  >  NC_000913/769976‑770446
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
cTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAg                                                                                                                                                                                                                                          >  1:450529/1‑239 (MQ=255)
 tCAACTATGTTGCCTTACAGGAAGGGCGTAACTGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACGCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAg                                                                                                                                                                                                                                          <  1:1174244/238‑1 (MQ=255)
   gctgaTGGTGCGTTAGATGAATGGCTTGACGGCATTGCTGTGGTGAGCGAAGGGTTTCGTTATTTTGAACGCGTCGGTGAAGAGAAGATAACCTTTGCCGTTACCTTGCTGCGTGGCGTGGGCTTTCTGGTCAAAGAAGATCTGCTTTTAAGGCATGGGCGGCCTTCCGGAATTAAAATGCCAGTTCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCATAAGTTTATTGAGTta                                                                                                                                                                                                                                       <  2:425800/235‑1 (MQ=255)
   aaCTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTAAGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCCCAACTACGTGGTCTGCTTTCTTGTCGCCTAAAGTTATTGa                                                                                                                                                                                                                                           >  1:170389/1‑235 (MQ=255)
    aCTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTtat                                                                                                                                                                                                                                      >  1:131123/1‑239 (MQ=255)
    aCTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGAATCACAAATACGGGGCCTGCTTTCTTGGTGCCTAAGGTTATTGAGTtat                                                                                                                                                                                                                                      >  2:60302/1‑239 (MQ=255)
     cTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTtata                                                                                                                                                                                                                                     <  2:207449/239‑1 (MQ=255)
        tGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGTTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTATTGTCGCCTACGTTTACTGAGTTATAccc                                                                                                                                                                                                                                  >  1:300257/1‑238 (MQ=255)
          ttGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGt                                                                                                                                                                                                                                >  1:1170830/1‑239 (MQ=255)
             ccTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTAc                                                                                                                                                                                                                              >  1:563740/1‑238 (MQ=255)
             ccTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGGGAAGAGAAGAAAACCTTTGCCCTTACGTTTCTGCGTGGCGGGGGCTTACTGGGCAAAGAAGATTTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAACATGCCAGGCCCGTACTCACAACTACGGGGTCTGCTTTCTTTACGCCTCAGTTTATTGAGTTATACCGGGAc                                                                                                                                                                                                                              >  2:34668/1‑238 (MQ=255)
                       aGGGCGTAACGGCATGGCTGTCTTTAGCGACGGGCTACGGGAATTTTGAGTCATCGGTGAAGAGAAGTAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACGACTACGTGGTCTGCTTTCTTGTCGCCTTAGTTTATTGAGTTATACCCGTACGCCAAccgcc                                                                                                                                                                                                                    >  2:217280/1‑238 (MQ=255)
                       aGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAAccgcc                                                                                                                                                                                                                    >  1:46859/1‑238 (MQ=255)
                       aGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTAGCCTAAGTTTATTGAGTTATACCCGTACGCCAAccccc                                                                                                                                                                                                                    >  2:724663/1‑238 (MQ=255)
                                   cATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGGCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTc                                                                                                                                                                                                       >  1:164001/1‑239 (MQ=255)
                                     tGGCTGGGTTTAGCGAAGGGGTACGTGAATTTGAAGTCATCGGGGAAAAGAAGAAAAACTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTTCTTTTAAGGCATGGGCGGCCTTTGGGAATTAAAATGCCAGGCCCGGAATCGCAACTAGGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGGGTTAAACCGGTACGCCAACCGCCGGTGGTGTAGCTca                                                                                                                                                                                                      >  2:549496/1‑238 (MQ=255)
                                            cTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGgcgc                                                                                                                                                                                              >  1:713992/1‑239 (MQ=255)
                                                   gaggggTTAGGTGAGTTGGAGGTCAGGGGTGAAAAGAAGAAACCGTTTGCCATTACGTGGCGGGGGGGGGGGGGTTGAGGGGGAAAAGAAGACCGGCTTTAAAGGCCGGGGCGCCTTGCGGGAATTAAAATGCAGGCCCGGGACTCACAACTCGGTGGCCTGCTTTCTGGTCGCCTAAGTTTATGGAGTTATACCGGTCCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCa                                                                                                                                                                                         <  2:122645/234‑1 (MQ=255)
                                                       gtttagggtattttGAGGTCATCGGTGAAAAGAAAAAAACCTTTTCCATTTCTTGGCTGCGTGGCGTGGGTTTACTGGCAAAAGAAGATCTGCTTTTAAGGCCGGGGCGGCTTTCGGGAATTAAAATCCCAGTCCCGGATCCACACCTACGTGGTCTGCTTTCTTGTCGCATAAGTTTATTGAGTTATCCCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                                                                                                                                    <  2:76195/227‑1 (MQ=255)
                                                       gtttACGTGAATTTGAAGTCATCGGTGAAGAAAAAAAACCCTTTGCCATTACGTTGCTGGGTGGCGGGGGTTTACGGGGCAAAGAAGATCTGTTTTAAAGGCCTGGGCGGCCTTCGGAAATTAAAATGCCAGCCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                                                                                                                                    <  1:82639/236‑1 (MQ=255)
                                                       ggttgtgTGATTTTGAAGGCATGGGGAAAGAGAAGAAAACTTTGGCCTTTCTGTGGCGGCGGGCGGTGGTTTTAGGGGAAAAAAAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGATTTAAAATGCCAGTCCCGAACTAACAACTGCGTGGTCTGCTTTCTGGGGGCTTAGGTTTATTGGGTCATCCCGGTCGCCCACCCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                                                                                                                                    <  2:563740/232‑1 (MQ=255)
                                                       ggTTACGTGAATTTAAGGTCATGGGGGAAGATAAAAAAACCTTTGCCATTACGTTGCGGCGGGGCGTGGGTTTACTGGAAAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATAAAAACGCCAGTCCCGGACTCACACCTGCGTGGTCTGCTTTCTTGTCGCCTAAGTTTATGGAGTTATCCCGGTCCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                                                                                                                                    <  2:1172698/238‑1 (MQ=255)
                                                        gTTACGTGAATTTGAAGTCATCGGGGAAGAGAAGAAACCCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGTTTTTAAGGCCTGGGCGGCCTTCGGAAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATCCCGTTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                                                                                                                                    <  1:836977/237‑1 (MQ=255)
                                                            cGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTcc                                                                                                                                                                              >  1:1034838/1‑239 (MQ=255)
                                                               aattttGAAGTCATCGGTAAAAAGAAAAAACCTTTTGCTTTCTCGTGGCCGGGGGGCGGGGGTTTAAGGGGAAAAAAAGTCGTGCTTTAAAGCCCGGGGGGGCCTTCGGGATTTAAAATCCAAGTCCGGGACTAAAAACTGCGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATCCGGGCCCGCCATCCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAg                                                                                                                                                                            <  2:788029/235‑1 (MQ=255)
                                                                 tttttAAGTCAGGGGTAAAGAGAAGAATACTTGTCCCATTTAGTGGGGGCTTGAGGGGGTCTTGATGGGAAAAGAAGATCTTCGTTTCAGACCTGGGCTGCTGACGTCACTTAAAATGCCAGTCCCGTTCTCACAACTTGTTGGTCTGCTCTTTTTTCTCTTCATTTAATTGAGGTTTACCGCTACGCCAACCGCCGCTGGTGTAGCTCAGCAGTCGCGAGCCTTGCTGACTCCAGTAc                                                                                                                                                                         <  2:831299/238‑1 (MQ=255)
                                                                 aTTTGAAGTCATGGGGAAAGAAAAAAAAACTTTTGCCTTTACGTGGCTGCGTGGCGTGGGCTTACTGGAAAAAAAAGTTCTGCTTTAAAGTCCGGGGCGGCCTTCGGGATTTAAAAGGCCAGCCCGGGACTCACAACTCGGTGCTCTGCTTTTTGGTCGCCTAAGTTTAGTGAGTTATCCCGGTCCGCCACCCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTAc                                                                                                                                                                         <  2:154584/239‑1 (MQ=255)
                                                                 aTTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTa                                                                                                                                                                          >  2:1023220/1‑238 (MQ=255)
                                                                 aTTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTAc                                                                                                                                                                         >  1:283621/1‑239 (MQ=255)
                                                                           aTCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTaca                                                                                                                                                                >  1:298138/1‑238 (MQ=255)
                                                                           aTCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTaca                                                                                                                                                                >  1:154584/1‑238 (MQ=255)
                                                                                      aaagaaACCTGTTGCCATTCCGTGGCGGGGGGGGGTGGGCTTACGGGAAAAAAAAGATCTGCTTTTAAGGCCGGGGCGGCCTTGGGGAATAAAAATGCCAGCCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATCCCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATccc                                                                                                                                                      <  1:279614/236‑1 (MQ=255)
                                                                                               cTTTGCCATTACGCTGCTGCGTGGCGTGGGCTTACTGGGAAAAGAAGATCTGCTTTTAAAGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTACCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGATATACCGGTACGCCAACCGCCGCTGGTGTACCTTGGCAGGCACGAGCATGGCCGTCTCCAGCAACGTGTTAAACGAAAATACCAtgggatggga                                                                                                                                           >  2:622005/1‑239 (MQ=255)
                                                                                                        tACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGCATATCACCTTTTAAGTCCTGTGCGGCCTATGGGAATTAAAATGCCAGTCCCGTACTCTCAACTACGGGGTCTGCTTACTTGTCTCCTAAGTTGATTGAGTCATAGCGGTAGGCTAACCGCCGCTGGCATAGCTCAGCAGCTGCGAGCATGACCGCATACACTACAGCGTTACACCTAAATACCATGGACTGTGCTGACGcccc                                                                                                                                  >  2:51432/1‑236 (MQ=255)
                                                                                                        tACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCAAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCCca                                                                                                                                  >  2:171633/1‑239 (MQ=255)
                                                                                                                    tggcgtggGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTa                                                                                                                                                                          >  1:34154/1‑187 (MQ=255)
                                                                                                                    tggcgtggGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACACAATCCCCTGGGATGTGATGAAGCCCAACAAAGCCGGAt                                                                                                                      >  2:421399/1‑239 (MQ=255)
                                                                                                                    tggcgtggGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGAt                                                                                                                      >  1:674317/1‑239 (MQ=255)
                                                                                                                    tggcgtggGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGAt                                                                                                                      >  1:258006/1‑239 (MQ=255)
                                                                                                                       cgtggGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCa                                                                                                                   >  1:628693/1‑239 (MQ=255)
                                                                                                                          ggtcTTATTGGGCAAAGAAGATCGGCTTTTAAGGCCTGGGCGGCCTTGGGGAATAAAAATCCCAGTCCCGGACCCACAACTGCGTGGTCTGCTTTCTTGTCGCTTAATTTTATTGAGTTATCCCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAAc                                                                                                                 <  1:696166/235‑1 (MQ=255)
                                                                                                                               tcgtGGGAAAAGAGAATCCGCTTTTAAGGCCGGGGGGGCCTTCGGAAATAAAAATGCCAGTCCGGGACTCACACCTACGTGGTCTGCTTTCTGGTCGCTTAATTTTATGGAGTTATACCGGTCCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCg                                                                                                           <  2:132559/236‑1 (MQ=255)
                                                                                                                                tatgggacaaaaagaCCTGCGTTTAAGGCCCGGGGCGGCCTTGGGAATTTAAAGCCCCAGTCCCTGACTCACACCTAGGTGGCTCGTTTCTTGGTCGCCAAATTTAATGGATTTATCCCGGAGCGCAAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCg                                                                                                           <  2:1218981/228‑1 (MQ=255)
                                                                                                                                 cggggAAAAGAAAATCTGCTTTTAAGGGCTGGGCGGCCTTGGGGAATAAAAATGCCAGTCCCGGACCAACAACGCGGTGGTCTTCTTTCTGTTCGCCTAAGTTTATTGGGTTAACCCGGGCCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGa                                                                                                         <  2:298138/237‑1 (MQ=255)
                                                                                                                                 cTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGa                                                                                                         >  2:30493/1‑239 (MQ=255)
                                                                                                                                 cTGGGAAAAGAAGATCTGTTTTTAAGGCCGGGGCGCCTTTCGGAAATTAAAATGCCAGTCCCGGACTCACAACTCCGTGGTCTTCTTTCTTGCCGCCTAAGTTTATTGATTTATCCCGGTACCCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGa                                                                                                         <  1:512387/239‑1 (MQ=255)
                                                                                                                                      cAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGtt                                                                                                    >  1:931344/1‑239 (MQ=255)
                                                                                                                                           aacatttcttttaAAGGCCTGGCCGGCCTTCGGTAATTAAAATCCCAGTCCCGGACTCACAACTACGGGTTCTTTTTTTTTGTCGCTTAATTTTATGGAGTTATCCCGGTCCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGt                                                                                               <  2:71610/230‑1 (MQ=255)
                                                                                                                                              acttgtttttAATGCCTGGGCGTCCTTCGGTAATAAAAATGCCATTCCCGCACCCACATCTCCGTGGTCTCCTTTCTTGTCGCCTAAGTTTATGGGGTTATGCCGGTCCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGa                                                                                                         <  2:308020/220‑1 (MQ=255)
                                                                                                                                              aTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTttg                                                                                            >  2:205259/1‑239 (MQ=255)
                                                                                                                                                     tttAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGaaaa                                                                                     >  1:848883/1‑239 (MQ=255)
                                                                                                                                                     tttAAGGCCGGGGCGGCCTGGGGGAATAAAAATGACATGCCCGGACTCACAACTCGGGGGTCTGTTTTCTTGTCGCTTAATTTTATTGTGTTAAGCGGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGaaaa                                                                                     <  2:46859/239‑1 (MQ=255)
                                                                                                                                                     ttgaaGGCCGGGGCGGCTTTGGGAAATTAAAATGCAAGCCCCGCACTCACAACTACGTGGTCTTCTTTCTTGTCCCTGAAGTTTATTGGTTTATCCCGGTGCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGaaaa                                                                                     <  1:34668/236‑1 (MQ=255)
                                                                                                                                                        aaGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGt                                                                                               >  1:308020/1‑226 (MQ=255)
                                                                                                                                                              cTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGCCTCCAGCACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCa                                                                             >  2:998508/1‑239 (MQ=255)
                                                                                                                                                                 ggcggcCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTg                                                                          >  1:98333/1‑239 (MQ=255)
                                                                                                                                                                     gttgtggggAATTAAAATGCCACTCCCGGCCTCACAACTACGTGGTCTGTTTTCTTGTCTCCGAATTTTATTGAGTTATGCCGTTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGa                                                                       <  1:135594/232‑1 (MQ=255)
                                                                                                                                                                     ggcgtCGGCAATTAAAATGACAGTTCCCGGTCGCCCAGCTCCGTTGTTTTCTTGCTGTTCGTCTATGTTAATGAGGTTAACCGGGTCCGACAACGGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGa                                                                       <  2:258006/235‑1 (MQ=255)
                                                                                                                                                                       gttcgagtgatAAAACGCCTCTCCCGTACTCCACACTGGGTGGTCTTCTTTTTTTTCGCGTAATTATAGTGATTTATACCGGACCCCAAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGc                                                                    <  2:450529/229‑1 (MQ=255)
                                                                                                                                                                                ttAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAgc                                                           >  1:547998/1‑239 (MQ=255)
                                                                                                                                                                                       gcaaGTCCGGGACTCACAACTAGGTGTTCTTTTTTCTTGCCGCCTAAGTTTATTGAATTATCCCGGTACGCAAGCCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACt                                                      <  1:221643/234‑1 (MQ=255)
                                                                                                                                                                                                   aCTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGaa                                         >  1:46015/1‑239 (MQ=255)
                                                                                                                                                                                                       acaACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGAcc                                     >  2:1116143/1‑239 (MQ=255)
                                                                                                                                                                                                                   tCTGCTTTCTTGTCGCCTAAGTTTATTGGGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAGGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGa                         >  1:1109010/1‑239 (MQ=255)
                                                                                                                                                                                                                   tCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGa                         >  1:367192/1‑239 (MQ=255)
                                                                                                                                                                                                                   tCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGa                         >  1:27467/1‑239 (MQ=255)
                                                                                                                                                                                                                   tCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGa                         >  1:1185163/1‑239 (MQ=255)
                                                                                                                                                                                                                   tCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGGGa                         >  2:566491/1‑239 (MQ=255)
                                                                                                                                                                                                                     tGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAAAGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGAACTTAAGAAAGCGGAAGACCGACAAGAAGGGAtt                       >  2:1036405/1‑239 (MQ=255)
                                                                                                                                                                                                                      gCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTACGAAAGCTGAAGACCGACAAGAAGTGCttt                      >  2:865929/1‑239 (MQ=255)
                                                                                                                                                                                                                      gCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGAttt                      >  1:317764/1‑239 (MQ=255)
                                                                                                                                                                                                                       cTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGAtttt                     >  1:754666/1‑239 (MQ=255)
                                                                                                                                                                                                                            ttGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGCAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGATTTTACgg                 >  2:474000/1‑238 (MQ=255)
                                                                                                                                                                                                                             tGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGATTTTACg                  >  1:1180271/1‑236 (MQ=255)
                                                                                                                                                                                                                                     tAAGTTTTTGGGATTATCCCGGTACCCCACCCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGATTTTACGGCTGTTTAAt        <  2:931344/239‑1 (MQ=255)
                                                                                                                                                                                                                                           tATTGAGTTATGCGGGTCCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGATTTTACGGCTGTTTAATCCGGCt  <  1:421399/239‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
CTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACGTGCCGGAAAGTTATAGTTTGTTGAAAATGCCCCCAGTGGGATGCCTGATAAGCGCACTTAAGAAAGCTGAAGACCGACAAGAAGTGATTTTACGGCTGTTTAATCCGGCT  >  NC_000913/769976‑770446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: