Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,529,180 T→C 8.8% pseudogene (1279/2037 nt) rhsE → pseudogene, Rhs family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,529,1800TC 8.8% 323.4 / 16.2 122pseudogene (1279/2037 nt)rhsEpseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base T (45/66);  new base C (4/7);  total (49/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.76e-01

Only 100 of 123 total aligned reads displayed.

ACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAA  >  NC_000913/1528907‑1529457
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                      
aCGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTAcc                                                                                                                                                                                                                                                                                  <  1:1046481/280‑1 (MQ=32)
              ggTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAAATGGCCAGACAGGTGGAGCCGGAATACAAACCGGCGCGAAAAGTTCATTTCTACCACTTCGACCACCgg                                                                                                                                                                                                                                                                    >  1:528405/1‑280 (MQ=18)
                tCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGCCAGGCTGGCGGTAGAAATCCGGGCAGACCGCGTGCGCCGTGAAAGCCGGGCGTGGCTTGCGCAGGGCGGGCTGACGGTGGAGCAACGGGCCAGCACGGTGGAGCCGGAATACAAACCGGCGAGCAAAGTTCATTTTTACCACTCCGACCACGGGGGc                                                                                                                                                                                                                                                                 >  2:80964/1‑281 (MQ=11)
                   agaCAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTg                                                                                                                                                                                                                                                              >  1:477712/1‑281 (MQ=34)
                   agaCAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGGTGAGCAACTGGCCAGACAGGTGGAGCCGGACTACACACCGGCGCGAAAAGTGCATTTTTACCACTGCCACCACCGGGGcca                                                                                                                                                                                                                                                               >  2:879440/1‑279 (MQ=19)
                                 gAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGgcg                                                                                                                                                                                                                                                     <  1:730895/276‑1 (MQ=34)
                                   gCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGTAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCtcatc                                                                                                                                                                                                                                               <  1:678788/280‑1 (MQ=25)
                                   gCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCtcatca                                                                                                                                                                                                                                              >  1:1039333/1‑281 (MQ=35)
                                   gCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGGGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGGGGAGCCGGAATCCACACACGCGCGGAAACGTCATTTTTACGACTGCGACCACCGGGGCCCGCCGCTCGCGCtcaca                                                                                                                                                                                                                                               >  2:550093/1‑278 (MQ=11)
                                         aaaaGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGaa                                                                                                                                                                                                                                         <  2:506194/280‑1 (MQ=35)
                                           aaGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATCTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGa                                                                                                                                                                                                                                       <  1:580086/280‑1 (MQ=18)
                                           aaGCGCAGCGGCCACGCTTGGCAGCGACGCGCCGGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGa                                                                                                                                                                                                                                       <  1:622951/280‑1 (MQ=255)
                                                 aGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTTAGCAGTGAAAGCCGGGAGTGGCTTGCGCAGTGAGGGCTGCCGGTGGAGCGACTGGCCAGACAGGGGGAGTCGGAATACACACCGGCGCGAAAAGTTCACTTGTACAACTTCGACGACCGGATCCTGCGTCCGGCGCTCAGCAGCGAAGACGGccag                                                                                                                                                                                                                                >  2:717232/1‑278 (MQ=11)
                                                       gcAGCCTGGCAGCGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGc                                                                                                                                                                                                                           <  1:150941/280‑1 (MQ=25)
                                                             tGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACAACGGCGCGAAAAGTTCATTTTTACCACCGCGACCCCCGGGGCCTGCCGCTGGCGCTCCTCCGCGAAGACGGCAATACGGCGTGgccg                                                                                                                                                                                                                    >  2:882617/1‑279 (MQ=255)
                                                               gcaaggaCGCTCCAGCAGGAAAGAGGTGGGAACGCCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGcgcg                                                                                                                                                                                                                   <  1:171825/275‑1 (MQ=255)
                                                               gCAGAGCCGCTCCAGCAGGAAGGGAGGGGGAACGGCCACGGCGTGGTGTTCCCTGCTGAACTGGGGCGGCTGCTGGACGGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCAGGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGTAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGcgcg                                                                                                                                                                                                                   <  2:465386/280‑1 (MQ=255)
                                                               gCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGcgcg                                                                                                                                                                                                                   <  1:882617/280‑1 (MQ=37)
                                                                cAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCggg                                                                                                                                                                                                                 <  1:307566/281‑1 (MQ=38)
                                                                 aggaccacaccAGCAGGAAGGGAGGGGGAAGGGCCAAGGCGTGGTGTGCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGGGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCggg                                                                                                                                                                                                                 <  2:368438/271‑1 (MQ=255)
                                                                     aCGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAACGTTCATTTTTACCACTGCGACCACCGGGGCCCGCCGCTGGCCCTCATCAGCGAAAACGGCAATACGGCGt                                                                                                                                                                                                                         >  1:764020/1‑268 (MQ=255)
                                                                     aCGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGGGGAGCAACTGGCCAGACAGGTGGCGCCGGAATACACACCGGCGCGAAAAGATCATTTTTACCACTGCGACCACCGGGGACTGCGCCGGGCGCCCAGCAGCGAAGAGGGCAATCGGACGTGGGGCGGGGGGTa                                                                                                                                                                                                            >  2:771504/1‑281 (MQ=11)
                                                                         tCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                                                                                         >  1:368438/1‑280 (MQ=39)
                                                                         tCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCCCATCAGCGAAGGCGGCAATACGGCGGGGCGCGGGGAGTagtgt                                                                                                                                                                                                        >  2:857283/1‑277 (MQ=255)
                                                                                       aGTGAGAACGGCCACGGGGTGGTGGTCCCGGCTGTACTGGTGCGGCTGCTCGACAGGCTGGGGGCAGAAATCCGGGCAGACCGCGTGAGGAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGATGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCaa                                                                                                                                                                                             <  2:774393/278‑1 (MQ=255)
                                                                                       aGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCaa                                                                                                                                                                                             <  1:914106/278‑1 (MQ=39)
                                                                                         tGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGCAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGATCAACTGGCCCGATAGGTGGAGCCTGAATACACACCGGCGCGAAAAGTTCATTTCTGCCACTGCGACGACCGGGGACTGCCGCGGGCGCTCCTTATCGAAGACGCCAATAAGGGGTGGCGAGGGGGGTATGCTGAGTGGGGAAACGAga                                                                                                                                                                                        >  2:21416/1‑280 (MQ=11)
                                                                                         tGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGTAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTTCCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCCTCACCGAAGACGGTAATACGGCGTGGCGCGGGGCGGGTGCAGAATGGGGCAACCAGc                                                                                                                                                                                        >  1:1004168/1‑281 (MQ=14)
                                                                                                ggCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGTGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGACGGAATACAAACCGGCGCGAAAAGCTCATTTTTACCACTGCGACTACCGGGGCCTGCAGCTGGCGCGCATCAGCGCAGACGGCAATACGGCGCGGCGCGGGGAGGATGATGAAGGGGGCAAACAGATTaagg                                                                                                                                                                                  >  2:258645/1‑278 (MQ=11)
                                                                                                 gCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGACAACGGCAATAAGGCGGGGGGCGGGGGGTAGGAAGAATGGGCGAACACGCTTaaggg                                                                                                                                                                                 >  2:341128/1‑277 (MQ=11)
                                                                                                  ccACGGCGTTGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCGCTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                                                                                <  1:415207/280‑1 (MQ=255)
                                                                                                     cGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGCTGAATGGGGCAACCAGCTTAATgaggag                                                                                                                                                                             >  1:186411/1‑280 (MQ=39)
                                                                                                      ggggtggggTTCCCGGCGGAACTGGGGCGGCTGCTGGACAGGCTGGAGGAGGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGGGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                                                                                           <  2:477712/278‑1 (MQ=255)
                                                                                                      ggCTTGGTGTTCCCGGCTGAACTGGTGCGGCTCCTGAACCAGCTCAAGGAATCAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                                                                                           <  1:717232/281‑1 (MQ=255)
                                                                                                           ggtgTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGCAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑c                                                                                                                                                                     <  1:258645/281‑1 (MQ=39)
                                                                                                           ggtgTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGAAGGAAGAAATCCGGGCAGACCGCGTGGGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑c                                                                                                                                                                     <  1:195372/281‑1 (MQ=255)
                                                                                                            gggtTCCCGGATGGACTGGGGCGGCTGCGGGACGAGCGGGGTGAATAAGTACGGGGAGAACGCGTGAGCAGTGGAAGCCGTGCGTGTCTTGCGCGGGGCGGGACGACGGGGGAGGAACGTGCCAGGCAGTTGGAACCTGAATAGACACCGGCGCGAAAAGTTCATTGTTACCACTGCGACCCCCGGGGCCTGCCTCTGGAGCTCATCCGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑Ca                                                                                                                                                                    <  2:570420/279‑1 (MQ=14)
                                                                                                                    ggCTGAACTGGTGCGGCTGCTGGACAGGCTGGCGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAACGCCGGGCGTGGCTTGCGCAGTGCGGGCTGCCGGTGGAGAAACCGGCCAGAAGGGTGGAGCCGGAATGCACACCGGCGTGAACAGTTCATTTTTACCAGTGAGACCACCGGGGCAGGCCGCTGGCGCTCATAAGCGGAGAAGGCAATACGTCGTGGCGCGGGGGGTATGGTGGATGGGACAAGCAGAGTAAAGAGGCAACCTCG‑CAGTACCTGc                                                                                                                                                            >  2:446613/1‑281 (MQ=14)
                                                                                                                     gCTGATCGGGTGCTGCGGCTGGGCGTGCTGGAGGAAAGAATCGGGGCTGACCGGCTGAGCAGGCAAGGCCGGGCGGGGCTTGCGCAGTGCGGGCTGGCGGTGCAGACACTGGCCAGACAGGTGGCGCCGGACTGCAAACCGGCGCGATAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGc                                                                                                                                                            <  2:233077/280‑1 (MQ=14)
                                                                                                                     gCTGAAGTGGTGCGGCTTCTGGACAGGCTGGAGGAGGAAATCCGGGCAGACCGCGTGGGCAGTGAAAGCCGGGCGTGGCTTACGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGGGAACCCG‑CATCACCTGc                                                                                                                                                            <  1:742581/280‑1 (MQ=255)
                                                                                                                      cTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCAc                                                                                                                                                          <  2:745280/281‑1 (MQ=255)
                                                                                                                        gAACTGGGTCGGCTGCTGGACAGGCTGGAGGAAGAAAGCCGGGCAGACCGCGTGGGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAGGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCa                                                                                                                                                        <  1:953705/281‑1 (MQ=255)
                                                                                                                                      gctgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCCGGCGCTCATCAGCGAAGACGGCAATCCGGCGTGGCGTGGGGAGTCTGATGCATGGGGAAACCAGCTAAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGAACCGTCTGc                                                                                                                                          >  1:693583/1‑281 (MQ=255)
                                                                                                                                              ggCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCCGGGCCTGCCGCTGGCGCTCATCAGCTAAGACGGCAATACGCCGTGGCGCGGGGCGGATGATGAATGGGGCAACCAGCTTAAAGAGGGGAACCCG‑CATCCCCCGCACCAGCCGGACCGTCTGCCAGGgca                                                                                                                                   >  2:834292/1‑280 (MQ=14)
                                                                                                                                                 tGGAGGAAGAAATCCGGGCAGACCGCGTGAGCCGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTGCATTTTTACCACTGCGACCACCGGGCCCTGCCGCTGGCGCCCATCAGCGAATACGGCAAGACGGCGTGGCGCGGGGCGTACGCTGAACGTGGCAACCAGCTTAATGAGGAGAACCGG‑CAACACCTGCCCCAGCCGCAACGTCGGCCAGGgctggag                                                                                                                               >  2:110801/1‑276 (MQ=14)
                                                                                                                                                 ggggggaggacagccGGGGAGACCGCGTGAGCAGTGAAGGCGGGGGGTGGCTTGCGCAGTGCGGGCTGCGGGTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAACGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGAAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagca                                                                                                                                <  2:802807/267‑1 (MQ=11)
                                                                                                                                                    aggaAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCat                                                                                                                            >  1:68524/1‑281 (MQ=255)
                                                                                                                                                        agaaTTCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCGTGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagcagca                                                                                                                             <  2:68524/276‑1 (MQ=255)
                                                                                                                                                         gaaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCAGGCGCTCATCAGCGAAAACGGCAATACGGCGTGGCGCGGGGAGTATGCTGAATGCGGCAACCACCTTAATGAGGAGAACCCG‑CACCACCTGCACCCGCCGTACCGTCTGCCAGGGCAGCAGCAcgatga                                                                                                                       >  2:832839/1‑281 (MQ=255)
                                                                                                                                                          aaaTCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgag                                                                                                                      >  1:777906/1‑281 (MQ=255)
                                                                                                                                                               cGGGCAGACCGCGGGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGt                                                                                                                  <  1:879440/280‑1 (MQ=255)
                                                                                                                                                                   cAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTCCCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                             <  1:680589/281‑1 (MQ=255)
                                                                                                                                                                   cAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                             <  1:158344/281‑1 (MQ=255)
                                                                                                                                                                    aGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                            <  1:611343/281‑1 (MQ=255)
                                                                                                                                                                    aGACCGCGTGAGCAGTGAAAGCCGGGCGGGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                            <  1:21416/281‑1 (MQ=255)
                                                                                                                                                                      aCCGCCTGAGCCGGGCAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTg                                                                                                          <  1:446613/281‑1 (MQ=255)
                                                                                                                                                                        cgcgTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCCGGCGCTCATCAGCGAAGACGGCAATACGGCGGGGTGCGGGGAGTATGAAGAATGGGGCAACCAGACTAATGAGGGGAACCCG‑GAACACCTGCACCAGCCGGACCCGCTCCCCGGGCCGCAGCAGGATGGGGAGTCGGGGCGGt                                                                                                         >  2:287611/1‑280 (MQ=14)
                                                                                                                                                                                cAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCCGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGGCACAGCCGTACCGTCTCCCAGGGCAGCAGCATGGTGAGGAGTCGGGGGTGTACGATAAcc                                                                                                >  2:877858/1‑281 (MQ=255)
                                                                                                                                                                                  gTGAAAGCCGGGCGTGGCTTGCGCCGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGACGGAATACACACCGGCGCGAAAAGTTCCTTTTTACCACTGCGACCACCGGGGTCTGCCGCTGGCGCTCATCAGCTAAGGCGACAACACGCCGTGGCGCGGGGAGTCTGCTGACTGGGGCAACCAGCGTAATGAGGAGAAACAGACGCCACAT‑CACCAACGGTACCGTCAGCCAGGGCAGCAGCAATAGGAGGGATGGGGTCTGTACGCTAACCg                                                                                               >  2:778108/1‑280 (MQ=14)
                                                                                                                                                                                   tGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGt                                                                                              <  1:211810/280‑1 (MQ=255)
                                                                                                                                                                                            ggCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACt                                                                                    >  1:457476/1‑281 (MQ=255)
                                                                                                                                                                                                          gcaaggggCTGACGGTTGAGCAACTGGCCTGGCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCa                                                                       <  1:205874/275‑1 (MQ=255)
                                                                                                                                                                                                             gCGGGGTGACGGGGGAGCAACTGGCCGGAAAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                   <  1:80964/281‑1 (MQ=255)
                                                                                                                                                                                                             gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCCCTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAAAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                   <  1:287611/281‑1 (MQ=255)
                                                                                                                                                                                                             gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                   <  1:978656/281‑1 (MQ=255)
                                                                                                                                                                                                             gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAggg                                                                    >  2:359958/1‑280 (MQ=255)
                                                                                                                                                                                                              ggggCTGACGGGGGGGGAACGTGCCAGACAGGTGGGGCCGGAATGCACACCGGGGCGAGAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCGGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAGGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCGGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                   <  2:817558/279‑1 (MQ=17)
                                                                                                                                                                                                                    gACGGTGGAGCGACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagca                                                                                                                                <  1:24110/213‑1 (MQ=255)
                                                                                                                                                                                                                      cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGACCCGTTGCAGGGGCGGTatat                                                           >  1:772596/1‑280 (MQ=255)
                                                                                                                                                                                                                      cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGACGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCCGGTCTATAACCGGCACCGGCACTCCGATCCGTTGCAGGGGCGGTagat                                                           >  2:691169/1‑280 (MQ=255)
                                                                                                                                                                                                                      cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGgcgc                                                                                                                                                                                                                    <  2:502409/128‑1 (MQ=255)
                                                                                                                                                                                                                      cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGgcgc                                                                                                                                                                                                                    >  1:502409/1‑128 (MQ=255)
                                                                                                                                                                                                                      cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTatat                                                           >  1:974803/1‑280 (MQ=255)
                                                                                                                                                                                                                      cGGTGCAGCAACTGGCCAGACGGGGGGGGCCGGAATACACCCCGGCGCGACAAGTTCATTTTTACCACTGCGACCACCGGGGCCGGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTatat                                                           <  2:277783/280‑1 (MQ=255)
                                                                                                                                                                                                                          gggcgaaCGGGCCAGACAGGGGGAGCCGGAATACACACCGGGGCGAAAAGCTCAGCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCAcc                                                       <  2:1047476/275‑1 (MQ=11)
                                                                                                                                                                                                                           gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCAcccc                                                     >  1:50623/1‑281 (MQ=255)
                                                                                                                                                                                                                           aaGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCGGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCAccc                                                      <  1:634481/279‑1 (MQ=255)
                                                                                                                                                                                                                               aacgtgCCAGGAAGGTGGGGCCAGAATACAACCGGGGGCGAAAAGCTCTACTTTATCACGGCGACCACCGGCGGTGGCCGCTGGCGTTTAGCGGCGGAGACGGCAATACGGCGGGGCGCGGGGAGTAGGATGAATGGGGCAACCAGCTTAATGAGGAGAACACG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTAAGATCCGTTGCAGGGGCGGTATATCACCCCGGa                                                  <  2:748525/275‑1 (MQ=14)
                                                                                                                                                                                                                                aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CCTCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGACACCGGCACCACGATCCGTTGCAGGGGCGGCATATCACCGCGGAcc                                                >  1:395734/1‑281 (MQ=255)
                                                                                                                                                                                                                                    gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCTCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGa                                             <  1:391886/280‑1 (MQ=255)
                                                                                                                                                                                                                                    gCCAAACAGGTGGGGCCGGAATACACCCCGGGGCGAAAAGTTCATTTTTACCCCTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGa                                             <  2:457476/280‑1 (MQ=255)
                                                                                                                                                                                                                                     cagggagagTGGAGCCTGGATGCACCCCGGCGCGAAAGGCTCATCTTGATCACTGGTGCCCGCGGGGCCTGCCCCTGGCGCTTATCAGCGAAGACGGCGATACGGCGGGGCGCGGGGAGGATGATGAATGGGGCAACCAACTTAATGGGGATAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGAtt                                           <  2:1033778/273‑1 (MQ=14)
                                                                                                                                                                                                                                      cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTg                                          >  1:884699/1‑281 (MQ=255)
                                                                                                                                                                                                                                      cAGACAGGGGGAGCCGGAACACACCCCGGCGCGAAAATTTCATTTTTACCACGGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGAtt                                           <  1:913703/280‑1 (MQ=255)
                                                                                                                                                                                                                                       aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGAtt                                           <  1:607508/279‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCCGAAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGagag                                 <  1:550093/278‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCAGGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGAACCGTGGCAGGGGCGGTAGATAACCCCGGACCCGACTgg                                         >  1:1033778/1‑269 (MQ=14)
                                                                                                                                                                                                                                                     ggAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtgga                            <  2:884699/280‑1 (MQ=255)
                                                                                                                                                                                                                                                     aatacacacacCGGCGCGAAAAGTTGTTTTTACCCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGGCGGCAATACGGGGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCGGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtgga                            <  2:7143/275‑1 (MQ=17)
                                                                                                                                                                                                                                                         gacaaacCGAAGCGAAAAGTTCATTTTTACCACGGCGACCACCGGGGGCTGCCGCTGGCGCTCATCAGCGAAGGCGGCAATACGGGGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCATAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGAGGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtggatgga                        <  1:158094/279‑1 (MQ=25)
                                                                                                                                                                                                                                                          acgcacCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGaa                       <  2:275935/277‑1 (MQ=255)
                                                                                                                                                                                                                                                                cGGCGCGAAAAGTTCATTTTTACCACTGCGACCCCCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTa                 <  1:834292/280‑1 (MQ=255)
                                                                                                                                                                                                                                                                cGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCAGATTGGGTTGAGAGGGGGGTGGAATATGTAt                >  2:658935/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                 ggCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATCAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGCTGCGAGGTGGATGGAATATGTATc               >  1:576619/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                   cgcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGGATATCACACCGGAACCGACTGGGTTGAGAGGTGGATGGAATATGTATCAg             >  1:1034367/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                   cgcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCCTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAACAGCTTAATGCGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGACAGGGCAGCAGCACGATGAGGGGTCGGGGCTGCACTATAACCGTCACCGGCACTACGACCCGTTGAAGGGGGGGTATATCGACCCGGACACGCGCGGGTTGATCGGTGGGTGGGATATGTATccg             >  2:620315/1‑279 (MQ=14)
                                                                                                                                                                                                                                                                    agtaaaaGTGCTATTTTACCCCTGCCACCACCAGGGCCTGCCGGTGGCGATAAGCAGCGAAGACGGCAGAACGGCGTGGCGCGGAGAGTAGTATGAAATGGGCAACCAGCGTAATGGGGAAAACCCG‑CATCACCTGCACCCGCCGTACCGTCGACCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTTAGAGGTGGATGGAATATGTATCAGt            <  2:1004168/278‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gAAAAGTTCATTTTTACCGCTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTAt          <  1:691169/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cATTTTTCCCACTGCGACCACCGGGGCCTGCCTCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGTGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGaa  <  1:9886/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 ttttACCTCGGCGACCACCGGGGCCGGCCGCTGGCGCTCATCAGGGAAGACGGGAAGACGGGGTGGCGCGGGGAGAATGATGAATGGGGCAACCAGCTTAAGGAGGAGAACCCG‑CAGCTCCTGAACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGaa  <  2:528405/278‑1 (MQ=18)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                      
ACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCG‑CATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAA  >  NC_000913/1528907‑1529457

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