Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,603,740 A→G 5.3% E241G (GAG→GGG)  lsrC → autoinducer 2 import system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,603,7400AG 5.3% 417.6 / 3.4 145E241G (GAG→GGG) lsrCautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base A (53/77);  new base G (4/5);  total (62/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 150 total aligned reads displayed.

CGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGTTAAACCCGCTTCGTCAGGTA  >  NC_000913/1603456‑1604027
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                  
cgAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATg                                                                                                                                                                                                                                                                                                <  2:635550/289‑1 (MQ=255)
cgAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATg                                                                                                                                                                                                                                                                                                <  1:896477/289‑1 (MQ=255)
      gAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGGGATGAAAGCa                                                                                                                                                                                                                                                                                          >  2:383245/1‑289 (MQ=255)
             ctctCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCa                                                                                                                                                                                                                                                                                    >  2:1060295/1‑288 (MQ=255)
             ctctCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCATGGCTGCTGGCAAAGACGGCGTTGGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACAGGGCGGGTGATGAAAGCAATTTCa                                                                                                                                                                                                                                                                                    >  2:1100858/1‑288 (MQ=255)
             ctcgCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAAGTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGCACGCAGTTTGTATGCCACGGGCGCTAATTTACAGGGCGCTCGTCAACTGTGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCa                                                                                                                                                                                                                                                                                    <  2:618094/288‑1 (MQ=255)
               ctcCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGGTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGc                                                                                                                                                                                                                                                                                  >  1:402158/1‑288 (MQ=255)
                          tgctgctTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTggg                                                                                                                                                                                                                                                                      <  1:1415527/289‑1 (MQ=255)
                                ttGGCGTTTCAGCAATTTGTTGGTGGACGAGAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTGGGACGAAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGAGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTGATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCa                                                                                                                                                                                                                                                                 <  1:117945/288‑1 (MQ=255)
                                       ttCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTg                                                                                                                                                                                                                                                         <  2:402158/289‑1 (MQ=255)
                                        tCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGc                                                                                                                                                                                                                                                        <  1:168413/289‑1 (MQ=255)
                                            cAAGTGGTTGGTTGACGATAATTCTGGTGGCATTTATGTCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTggg                                                                                                                                                                                                                                                                      <  2:306137/271‑1 (MQ=255)
                                             aaTTGGTTGGTTGACGATAATTCTGGAGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTGTTTATGCCACGGGCGATAATTTACAGTGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCgg                                                                                                                                                                                                                                                    <  2:898582/288‑1 (MQ=255)
                                               ttggttggTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGGGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtgg                                                                                                                                                                                                                                                 <  1:383245/289‑1 (MQ=255)
                                               ttggtgggTTGACGATAATTCTGGTGTCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGAAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtgg                                                                                                                                                                                                                                                 <  2:1563950/289‑1 (MQ=255)
                                                 ggttggttGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGGTGGAGATGAAAGCAATTGCAGCCCGCGTGCTGGGGGGCATTAGTTTGCTCggtggt                                                                                                                                                                                                                                                >  2:81771/1‑288 (MQ=255)
                                                   ttggttgACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTc                                                                                                                                                                                                                                              <  2:1008960/288‑1 (MQ=255)
                                                     ggttgCCGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCgg                                                                                                                                                                                                                                           <  2:1179171/289‑1 (MQ=255)
                                                      gttgACGATAATTCTGGTGGCTTTTATGGCCTGGCTGCTGGCTAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCTATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGCTCTTGGGTTTGCTTTGCCGCTTGGTTTTATTGCCAACCAGAACGGTACCGCGCTGGGACTGAAAACAATTTCATCATGCGTGTTTGGCGGCCTTCGtgctgtccctttttctcat                                                                                                                                                                                                                                          >  1:306137/1‑271 (MQ=255)
                                                            gATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGGGGCACTGGCGGGAATTGGGTTTGCTTCGCAGATTGGGTTTATCGCCTACCAGGCCTGGACCGGGCTGGGGATGGAAGCAATTGTAGGCTGCGTGCTGGGCGGCAATAGTTTGCTCGGTGTTTCTggggagtgt                                                                                                                                                                                                                                    >  1:1522492/1‑282 (MQ=255)
                                                             aTAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCCAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCCGACCGGTCCCGGGCTGGCGGTGAACGCACTTGCAGCCTGCGTGCTGGGCGGCATTCGTTTGCTCGGTGGTTTCGGTGCgccca                                                                                                                                                                                                                                   >  2:742272/1‑285 (MQ=255)
                                                             aTAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCa                                                                                                                                                                                                                                   >  2:1139048/1‑289 (MQ=255)
                                                                 tgcTGGTGTCAGTTATGGCCTGGCTGGTGGAAAAGACGGCGGTGGGGCGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCGCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAGCGGCTGCATGGCGGCACTGGCGGGAGTTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGTGATGAAAGCAAGTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTg                                                                                                                                                                                                                                <  1:48222/286‑1 (MQ=255)
                                                                          agagacaagCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACt                                                                                                                                                                                                                     <  2:1066710/281‑1 (MQ=255)
                                                                           aTTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTCGCATTTTCGTTGAACGGCTGCATGGCGGCCCTGGCGGGTCTTGTGTTTGCTTCGCAGATTGGTTTTATCCTCAACCCGACCGGTACCGGGCTGGCGATGAAAGCAATTGCCGCCTGCGTGCTCGGCGGCCTTCGTTTGCTCGGAGGTTTGGGTGCGATCGTTGGTTCGGTACTc                                                                                                                                                                                                                    >  1:187771/1‑289 (MQ=255)
                                                                            tttATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGGTCCGGTGCGAGCATTTGGGCGGGACTCg                                                                                                                                                                                                                   >  1:178449/1‑289 (MQ=255)
                                                                                           tgggCAAAGTCGGCGTTTGGACGCAGGTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGa                                                                                                                                                                                                    <  1:81771/287‑1 (MQ=255)
                                                                                             ggCAAAGACGGCGTTTGGACGCGGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                                                                                                                                                                                  <  2:1494108/289‑1 (MQ=255)
                                                                                             ggCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGAc                                                                                                                                                                                                   >  1:698514/1‑288 (MQ=255)
                                                                                             ggCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                                                                                                                                                                                  >  2:124885/1‑289 (MQ=255)
                                                                                             ggCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGGGTTTGGTTCGCCGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCCGGGGCTGAAAGCAATTGCACCCTGTCGGCTGGGCGGCGTTTGTTTGGTTGGTGGGTCTGGTGCGATGTTTGGTGCCGGTCCTCgggcaggggttggtgct                                                                                                                                                                                                   >  1:1066710/1‑274 (MQ=255)
                                                                                              gCAAAGGCGGCGTTTGGACGCAGTTTTTATGCCACGGGCGAGAGTTTACAGGGCCCTCGCCAAATGTGCGTTCGTACAGAATCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                                                                                                                                                                                  <  2:1510272/288‑1 (MQ=255)
                                                                                              gCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCACTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                                                                                                                                                                                  <  1:1535034/288‑1 (MQ=255)
                                                                                                  aGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGAt                                                                                                                                                                                             <  1:491771/289‑1 (MQ=255)
                                                                                                     cGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGAATGGTTCCTGACGCAGATCGa                                                                                                                                                                                          <  1:1139048/289‑1 (MQ=255)
                                                                                                      ggCGTTTGGACGCGGTGTTTATGCCTGTGGCAATGATTTGCATGGCGGTCCTGAACGTGGAGTTCTGACGGAAGCCATTCGCATTGTGGGATTTTCGTTGGGCGTCGCCGTGGAGGCACTGGCTGGAATGGTGTTTGCTTCGCAGAGTGGTTTTTTCCTCAACCAGACCGGTACCGGGATGGAGATGAAAGCAATTGGAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTAGGCGCATGGTTCCTGCCGCAGATCGa                                                                                                                                                                                          <  1:1202382/288‑1 (MQ=255)
                                                                                                        cGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGGGCGGTACTCGGCGCATGGTTCCTGa                                                                                                                                                                                                    >  1:1051437/1‑276 (MQ=255)
                                                                                                           ttGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCg                                                                                                                                                                                     <  2:337597/288‑1 (MQ=255)
                                                                                                                  cAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTg                                                                                                                                                                             >  2:1277060/1‑289 (MQ=255)
                                                                                                                   aGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTg                                                                                                                                                                             <  2:698514/288‑1 (MQ=255)
                                                                                                                   aGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTg                                                                                                                                                                             <  2:662258/288‑1 (MQ=255)
                                                                                                                        ttATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCCGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGTTGCGGCACTCGGGCCCTGGTTCCTGACGCCGCTCTATTGCGTGCAGGTGCTGtt                                                                                                                                                                        >  1:652306/1‑288 (MQ=255)
                                                                                                                                     gATAATTTACCGGGCGTGCGTCAACTGGCCGTTGGTACGGAAGCCATTGGCATTGTGGCATTTTCGTTGAACGGCGACATGGCGGCACTGGCGGAAATTGTGTGTGCTTTGCAGATTGGTTTTATCCCCAACCTGAGCGGTACCGTGCTGGAGATGAAAGTAATTGCAGCCCGCGTGCTGGGCGGCATTAGTTTGCTCTGTGGGTCCGGTGCGATCATTGGTGCGGTACTTGTCGCATTGTTCCTGACGTAGATTCATAGCGTACTGGTGCTGTTGCGCATTCCGGCa                                                                                                                                                          <  1:1100858/288‑1 (MQ=255)
                                                                                                                                     gATAATTTACAGTGCGCTCTTCAACTGCGCGTTCGTACTGAAGCCAGTCGCAGTGTTGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGTTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCa                                                                                                                                                          <  1:1149083/288‑1 (MQ=255)
                                                                                                                                     gATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCa                                                                                                                                                          <  1:1277060/288‑1 (MQ=255)
                                                                                                                                                 ggCGCTCGTCAACTGGGCGTTCGTACTGAAGCGATTCGCATTGTGGCATTTTCGTTGAACGGCGGCATGGCGGCACTGGCGGGAATTGTGTTTGATTCGCAGATTGGTTTTATCACCAACCAGACCGGTACCGGGCTGGAGATGAATGGAGTTGAATGCTGCGTTCCGGGCGGCAATAGGTTTCGCGGTGGGCTCGGGTCGATCATTGATGAGGTCGATGGTGCAGGGTTTCTGCAGCACAACGATAGCGTGCTGCTGCTGTTGCGCATTCCGGCATGGAGGAAAGAtt                                                                                                                                             >  1:1303005/1‑289 (MQ=255)
                                                                                                                                                 ggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGAt                                                                                                                                                                                                                                     >  2:852193/1‑203 (MQ=255)
                                                                                                                                                 ggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGAt                                                                                                                                                                                             <  1:649520/242‑1 (MQ=255)
                                                                                                                                                   cgcTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCATCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCCGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGCATGAtttt                                                                                                                                           >  2:559341/1‑289 (MQ=255)
                                                                                                                                                   cgcTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtttt                                                                                                                                           >  2:1426966/1‑289 (MQ=255)
                                                                                                                                                        gTCAACTGGGGGTTCGTACTGTAGCCAGTTGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGAGTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAAATGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGCTCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATcgc                                                                                                                                      <  1:1006562/289‑1 (MQ=255)
                                                                                                                                                        gTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATcgc                                                                                                                                      >  2:229959/1‑289 (MQ=255)
                                                                                                                                                        gTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATGGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATcgc                                                                                                                                      <  2:461456/289‑1 (MQ=255)
                                                                                                                                                               gggCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCCGGGCGGCCTTAGTTTGCTCGGTGGTTCCGGTGCCATCATTGGTGAGGTACGAGGGGCATGGTTGCTGACGCCGTTCGATAGCGTACTGGGTGTGTTGCGCTTTGCGGGATGGTGGAATGGTGTTACGCCGGGTCTg                                                                                                                              >  1:29157/1‑289 (MQ=255)
                                                                                                                                                                ggCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGGATTGTGTTTGCTTCGCAGATTGGTTTTCTCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGCTCCGGTGCGATCATTGGTCCGGTACTCGGCGCATGGTTCCTGACGCAGATGGATAGCGTACTGGTGATGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGTGGGTCTgg                                                                                                                             >  2:1186356/1‑289 (MQ=255)
                                                                                                                                                                  cGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGtt                                                                                                                           >  2:960375/1‑289 (MQ=255)
                                                                                                                                                                   gTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTAACGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTCGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCCTGGTTCCTGACCCCGATCGATCGCGTACTGGTGCTGTTGCGCATTCCGGCATGGCTGAATGATTTTATCGCggggttggg                                                                                                                             >  1:265668/1‑281 (MQ=255)
                                                                                                                                                                         agtGAAGCCATTCGCATTGTGGGATTTTCGTGTAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                                                                                                                    <  1:994081/287‑1 (MQ=255)
                                                                                                                                                                         aCTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                                                                                                                    <  2:895194/289‑1 (MQ=255)
                                                                                                                                                                              aGCCATTCGCATTGTTGGTTTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCGGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAACGCAATGGCAGCCTGCGTGCTTGGCGGCATTAGTTTGCTCGGTGGTGCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTg                                                                                                               <  1:346925/289‑1 (MQ=255)
                                                                                                                                                                                ccATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGaa                                                                                                                                                  <  2:723324/253‑1 (MQ=255)
                                                                                                                                                                                  aTTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCCACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCCGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGATGCAGATTCATCGCGTACTGGTGCTGTTGGCGCTTCCGGCATGGGGGAATGGTTTTCTTGCGGGTCCTGGtttggttggggggctgg                                                                                                            >  1:1371130/1‑274 (MQ=255)
                                                                                                                                                                                   ttCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                                                                                                                                                                                  >  1:123470/1‑203 (MQ=255)
                                                                                                                                                                                   ttCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                                                                                                                                                                                  <  2:123470/203‑1 (MQ=255)
                                                                                                                                                                                          tgtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                                                                                                                                                                                          <  2:1084717/157‑1 (MQ=255)
                                                                                                                                                                                           gtgGCATTTTGGTGGGACGGCTGGATGGCGGCACTGGCGGGAGTTGTGTTTGCGTGGGAGATTGGTTGTCTCCCCAACCAGACCGGTACCGGGCTGGGGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCGCGGTGGTTCCGGTGCGATCATTGGTGTGGTTCTCGGCGCATGGTTCCGGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATg                                                                                                  <  1:1075085/289‑1 (MQ=255)
                                                                                                                                                                                              gCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGa                                                                                                <  2:1031069/288‑1 (MQ=255)
                                                                                                                                                                                                aTTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGc                                                                                             <  2:301091/289‑1 (MQ=255)
                                                                                                                                                                                                  tttCGTTGACCGGCTGCATGGGGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCt                                                                                           <  2:1108995/289‑1 (MQ=255)
                                                                                                                                                                                                  tttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                                                                                                                                                                                          <  1:1468691/149‑1 (MQ=255)
                                                                                                                                                                                                  tttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGAAGCAGATCGATAGCGTACTGGTGCTGTTGCGAATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCt                                                                                           >  2:1517049/1‑289 (MQ=255)
                                                                                                                                                                                                      gTTGAACGGCTGCATGGCGGCACTGGCGGGAACTGTGTTTGCCTCGCAGATTGGTTTTATCCCCAACCCGACCGGTACCGGGCTGGAGAAGAAAGCAATTGCAGCCTGCGTGCTCGGCGGCATTAGTTAGCTCGGTGGTTGCGGTGCGATCATTGGTGAGGTACTCAGCGCAGGGATCCCGACGATAATCGAAAGAGGACCGGAGCTGTTGCGAAATCCGGCATGGTGGAACGATTTTATTCAGGGTCTTGTTCTGGTTGCGGTGCCggcggtttggggcagcacgcgt                                                                                       >  1:837240/1‑269 (MQ=255)
                                                                                                                                                                                                                    tGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACg                                                                         >  1:389506/1‑289 (MQ=255)
                                                                                                                                                                                                                            cTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTAc                                                                 >  2:392866/1‑289 (MQ=255)
                                                                                                                                                                                                                                   gAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACggcggc                                                           <  2:1562223/288‑1 (MQ=255)
                                                                                                                                                                                                                                   gAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACggcggc                                                           <  2:1374283/288‑1 (MQ=255)
                                                                                                                                                                                                                                         tgtTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTtgtg                                                                                   <  1:938128/258‑1 (MQ=255)
                                                                                                                                                                                                                                         tgtTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACGGGTGCTGTTGCGCATTCCGGCATGGTGGAATGGTTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGTTGGACGCCTGCGTTGTGCGCTGGAACGCAATCTCAGGCGGCAAAAata                                                    >  1:1238875/1‑289 (MQ=255)
                                                                                                                                                                                                                                           tttGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGc                                                                                             <  2:954969/246‑1 (MQ=255)
                                                                                                                                                                                                                                           tttGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGc                                                                                             >  1:954969/1‑246 (MQ=255)
                                                                                                                                                                                                                                             tGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGGTCGCGGGAACGTAATCCTGCGCGGCAAAACTTTGc                                                 >  1:1485772/1‑288 (MQ=255)
                                                                                                                                                                                                                                             gggtgtgaaGATGGGTTGTAACCCCAACCTGAACGGTGCCGGGCTGGAGATGGAGGAAATTGGAGAATGCGGGTTGGGCGGCATGAGTTGGCGCGGTGGGTCCGGGGCGACCATTGGGGTGGTACTCGGCGGATGGTTTCTGACGTAGATGGATAGCGGACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATCATTGTATCGCGGGTTTGGTTCTGCTGTCGGTGCTGTTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATAGGcc                                                <  1:742272/281‑1 (MQ=255)
                                                                                                                                                                                                                                              gCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                                                                                                                    <  2:807446/220‑1 (MQ=255)
                                                                                                                                                                                                                                                              ttATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTTTTTTGTGGAAGCCTGCGTTTTGCGCTGGAACGTAATCTAAGGCGGCAAAACTATGCCCGCTTTATGAAGCCAcc                               >  1:194606/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                 tcccaggccATACCGGTACCGGGGTGGAGATGAAAGCAATTGCAGCCTGCGGGCTGGGCGGCATTATTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCTCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGcc                            <  2:793624/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTgcg                                                                                 >  1:1189829/1‑236 (MQ=255)
                                                                                                                                                                                                                                                                   cccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCAt                          <  2:237013/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATcg                                                                                                                                       >  1:857861/1‑179 (MQ=255)
                                                                                                                                                                                                                                                                      aaCCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGAGCGATCATTGGTGCGGTACTCTGCGCAGGGTACCTGACGCAGATCGACAGCGTACAGGTGCAGGTGGGCATTCCGGCATGGTAGAATGAATTTCTCGCGGGTCGGGTTCTGGTGGCGGTGCTGTGGTTTGAGGGATGACTGCGTTGTGCGCGTGGACGTACTCTAGGGCGGCGAAAATAGGCCCGTCTTAttgtgtcccccacttcc                        >  1:1101309/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                            aCCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCTTTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGa                                                                                                                                                                                          <  1:1214261/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aCCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCTTTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGa                                                                                                                                                                                          >  2:1214261/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                             ccGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGTTaaa                  <  2:1485772/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCCTGGTTCCTGACGCAGATCGATAGCGTACTGGTTCTGCTGCGCATTCCGGTCTGGTGGAATGGTTTTATCGCTGGTCTGGTTTTGCTGGCGGCGCTGGTGTTTGATGGACGCCGGCGCTGTGTGCTGGACTGTAATCCACGACCGCAAAAATATGCACGACTTATGACGCCACCGCCATCCGTTAAACACg              >  1:499348/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                 tACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGAt                                                                                                                                                                                             <  1:32870/114‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     gggCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGTTAAACCCGCTTCGt        <  2:1123995/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      ggCTGGGGATGAAGGAAATTGCCGCCTGCGGGCTGAGTGGCCTTAGTTTGCTTGCTGGTTCGGGTGCGATCATGGGGGCGGTACTCGGCTCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGTCATGGTGGAATGATTTTATCGCGTGTCTGGTTCTGCTGGCGGTGCTTGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGTTAAACCCGCTTCGt        <  1:477920/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      ggCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGTTAAACCCGCTTCGTc       >  2:505787/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                           gagaTGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTAt                                        <  2:485396/251‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           gagaTGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGTTAAACCCGCTTCGTCAGGTa  >  2:1052063/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                           gagaTGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCCATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTCCTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGCTTTTATCGCGGGTCTGGGTATGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCAGCGTTATGACGCCATCGCGATCCGTTAAACCCGGTACGCCAggca  >  1:96493/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                  
CGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGTTAAACCCGCTTCGTCAGGTA  >  NC_000913/1603456‑1604027

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