Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,247 A→G 4.4% intergenic (+366/‑212) mhpE → / → mhpT 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913375,2470AG 4.4% 290.6 / 3.0 92intergenic (+366/‑212)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (39/47);  new base G (2/2);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-01

CTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGTCGCTC  >  NC_000913/374960‑375531
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                            
ttACGGAGGGCGTGAGAATTTGTAGGTCGGGTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGGCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTa                                                                                                                                                                                                                                                                                              <  1:114669/287‑1 (MQ=255)
         gCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGGCTGCTCAGACGCG‑TTAGCGTCGCATCGGGCATCTGCGCACGACTGCCGGCTGCGGCGTCAACGCCTTaa                                                                                                                                                                                                                                                                                    >  2:622248/1‑288 (MQ=255)
                aaTTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGCTGGCGCGGGAATTTGTCGGCCTGCTAAGACACG‑TTAGCGTGCCATCAGGCATCTGCGCACGCCTGCCCGCTGCGGCGTaagctacgtattcggct                                                                                                                                                                                                                                                                              >  1:202813/1‑273 (MQ=11)
                          tcagctAAGGCGCGTGAGCGTCGCAGCAGGCGTTTGCGCACGACTGCCGGAGGCGGCGGGAGGGTCTTCTCCGGCATACGGATGGCGCGAGAATTTTTGGGCCTGATAAGAGGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGGGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                                                                                                                                                                   <  1:809867/284‑1 (MQ=14)
                               taggACGCGTTAGAGACGCACAAGGCATGTGGGCACGACTGCCGGACGCGGCGTGACCGGCTTATCCGGCGTACGGATGGCGCGGGAAATTGTAGGACTGATAAGAGGCGTTAGCGTCGCATCAGGCATCTGCGGACGACTGCCGGATGCGGCGGGAAC‑GCCTTATCCGGCCTACGGGTGGCGGGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              <  1:147681/286‑1 (MQ=14)
                                          aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGGGCGGGCATTTGTAGGc                                                                                                                                                                                                                                                    >  2:594785/1‑288 (MQ=255)
                                                        cATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                                                                                                                                                                       <  2:333590/287‑1 (MQ=255)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGATGGACCGGGAATTTGTAGGCCTGATaacatacgtacgcg                                                                                                                                                                                                                                 >  1:84440/1‑277 (MQ=255)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              <  1:889671/260‑1 (MQ=255)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              >  2:889671/1‑260 (MQ=255)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                                                                                                                                                                                 >  1:627605/1‑289 (MQ=40)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGCTGGCGCGAGAATTTGTAGGCCTGATAAGAAGCG‑TTAGCGTCGCATCATGCATCTGCGAACGACTGCCGGACGCGGCGTACACGCTCTATCCGGCCTACGGACCGCCCGGCCATTTACAGGCATGAtcatactacgttgtc                                                                                                                                                                                                                                 >  1:126074/1‑275 (MQ=11)
                                                                 aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCa                                                                                                                                                                                                                            >  2:864652/1‑289 (MQ=39)
                                                                  cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACACTTCGCCCG‑TAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTAcactt                                                                                                                                                                                                                                                                       >  2:78999/1‑240 (MQ=11)
                                                                  cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACACTTCGCCCG‑TAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTAcactt                                                                                                                                                                                                                                                                       <  1:78999/244‑5 (MQ=11)
                                                                       gcgggtgggggcgGGAAGGCCTTGTGCGGCGTACGGATGGCGGGGGAATTTGTAGGCCGGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAgg                                                                                                                                                                                                                       <  1:197416/279‑1 (MQ=11)
                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                                                                                                                                                                                                      <  2:627605/289‑1 (MQ=39)
                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTCGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                                                                                                                                                                                                      <  2:408581/288‑1 (MQ=21)
                                                                               cggcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAACTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGATTGCGCGGGAATTTGTAGGCCTCATAAGACGGGTTAGCGTCCCATCAGGCATCTGca                                                                                                                                                                                                               >  1:626808/1‑288 (MQ=255)
                                                                                              aTCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                                                                                                                                                <  2:316130/256‑1 (MQ=40)
                                                                                              aTCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                                                                                                                                                >  1:316130/1‑256 (MQ=40)
                                                                                                   gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              <  1:801647/221‑1 (MQ=255)
                                                                                                   gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              >  2:801647/1‑221 (MQ=255)
                                                                                                   gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCCGCCCACGACTGCCGGGTGCGGTaa                                                                                                                                                                                          >  1:425315/1‑289 (MQ=255)
                                                                                                     cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAAc                                                                                                                                                                                        <  2:780792/289‑1 (MQ=255)
                                                                                                       aCGGATGGCGCGGGAATTTGTAGGCCTGAGAAGACGCGTTGGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                                                                                                                                                                                      <  2:352049/289‑1 (MQ=255)
                                                                                                            tGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 <  2:16053/289‑1 (MQ=255)
                                                                                                              tggcggacaTTTGTAGGCTGTGTGAGACGCGTTGCAGACGCATCGTGCACCGGCGCACGTCTGCCGGCTGCGG‑GTGGACGGCCTTTGCCGGCCTAAGGTTGGCGGGAGAATTTGTAGGCCTGATAAGACGCG‑TTGGCGTCGCATCGGGCAGCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                                                                                                                                                <  1:622248/280‑1 (MQ=14)
                                                                                                              gggcgGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATa                                                                                                                                                                                                                                              <  1:570506/224‑1 (MQ=255)
                                                                                                              gcgcGGTAATGGGGAGGCGGGGTAGGACGGGGTAGGGGCGGATGAGGCATCGGCGCACGACTGCCGGGTGGGGCGGGAGC‑GCCTTATCCGGCCGACGGGTGGCGCGAGGATTTGTAGGGCTGAGAAGAGGCG‑TTAGCGTCGCATCGGGCGTCTGAGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                                                                                                                                                <  1:805332/288‑1 (MQ=14)
                                                                                                              gcgcGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCCTCAGGCATCTGCGCACGACTGCCGGATGCCGCGGAAAATCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGGGTTGACGTCGCACGCGGCCTCTGCACCCGGCCCCCGGATCCGATAAACggcatttc                                                                                                                                                                                >  1:352049/1‑281 (MQ=14)
                                                                                                              gcgcGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATa                                                                                                                                                                                                                                              >  2:570506/1‑226 (MQ=255)
                                                                                                                                 tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAACCGCCTTATCCGGTCTACGGATGGCGCGGGAATTTGTAGGCCTGTTAAGACGCGTTCGCGTCCCATCAGGCATCTGTACACCGCTTCCGTCTGCGATAAGCGTCTTGTGCGGCCtcactttccctcgt                                                                                                                                                             >  1:378983/1‑275 (MQ=14)
                                                                                                                                  gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTG‑‑AAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 >  1:488315/1‑265 (MQ=255)
                                                                                                                                  gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTG‑‑AAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 <  2:488315/265‑1 (MQ=255)
                                                                                                                                    tAAGAGGCGTTAGCGTGGCATCAGGCAGCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                                                                                                                         <  1:733655/289‑1 (MQ=255)
                                                                                                                                    tAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                                                                                                                         <  2:84440/289‑1 (MQ=255)
                                                                                                                                         cgcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                                                                                                                                                    <  2:707577/289‑1 (MQ=255)
                                                                                                                                          gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATACGACGCGTTAGCGTCCCCTCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCCGTcct                                                                                                                                                   >  2:383352/1‑287 (MQ=255)
                                                                                                                                          aacgtgcattttgtatgcccatcccagcACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATACGACGCG‑TTAGCGTCGCATCAGGCCTCTGCGCACGACTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGATGGCGCGGGCATTTTTAGGCCTGCTACGACGCGTTACCGTCCCATCCTGCACATGCCAACGACTTCCGGATGCGATAACCTGCCTGTCCGGCCTACATTTCTCCCCTAGGCAGTcct                                                                                                                                                   >  1:793117/27‑287 (MQ=11)
                                                                                                                                            gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTa                                                                                                                                                                                                                                                                                              <  2:836658/148‑1 (MQ=35)
                                                                                                                                              tAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                                                                                                                                               <  2:416394/289‑1 (MQ=255)
                                                                                                                                                     gCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCt                                                                                                                                        >  2:400260/1‑289 (MQ=255)
                                                                                                                                                      cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                                                                                                    >  1:836658/1‑148 (MQ=37)
                                                                                                                                                             cagcTGCGAACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCATACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcg                                                                                                                                 <  1:923602/285‑1 (MQ=255)
                                                                                                                                                                 ttaccACGACGGCCGTATGCGGCGTGCAC‑GCCTTGTCCGGCCTACGAGTGGCGCGAGAATTTGTAGGCCTGATAGGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc                                                                                                                             <  2:505115/284‑1 (MQ=255)
                                                                                                                                                                 cgggcggaccgtgcgtaTGCGGCGCGACC‑GCCTGATCCGGCATACGGTGGGAGGGATAATGTGTAGGCCTGGTAAGACGCG‑TTAGCGGCGCATCATGCATCGGCGCGCGACGGCCGGATGCGGCGTAATCGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGCACACGACGGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCt                                                                                                                            <  1:512868/273‑1 (MQ=255)
                                                                                                                                                                   cgcACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCCCGGGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGCCATCTGCACACGACTGCCGGCTGCGATCACCGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCCTTAACTTGTTCTTATTTcgcc                                                                                                                               >  1:751004/1‑283 (MQ=17)
                                                                                                                                                                    gcACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTgg                                                                                                                          <  2:425315/288‑1 (MQ=255)
                                                                                                                                                                    acACGACTGCCGGATGGGGCGTGAAC‑GCCTTATCCGGCCTACGGGGGGCGCGAGAAGTTGTAGGCCTGATAAGACGCG‑TTAGCGTGGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTgg                                                                                                                          <  1:886655/287‑1 (MQ=255)
                                                                                                                                                                      aCGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGCTAAGACGCCTTAGCGGCCCATCAGGCATCTGCCCCCGACTCCCGCATGCACTAAACGGCTTGTTCGGCCTACCTTTCCCCCGTCGGCAGTCATCAATCAGCTCTTCTTACCCGTCTGCCac                                                                                                                       >  1:665419/1‑288 (MQ=14)
                                                                                                                                                                      aCGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCCCCCGTAGGCCGCCATTAAATAGTTCTTATTACGCGCCGGGcca                                                                                                                       >  1:409093/1‑287 (MQ=255)
                                                                                                                                                                       cGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGTCTACACTTCGC‑CCATAAACCGTAGGCCTGATCAGACGCG‑TTAGCGTCGCATCAGGCATCTCCGCACGAATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCCTACGCAAACGACGCACAGATGCGATAAACGTCCTGTCAGGCCGAGAATCCCCCCGGCGGCAGCCACGAAACACTTCTGCTTATGCTGCTGGCaact                                                                                                                     >  1:666201/1‑287 (MQ=9)
                                                                                                                                                                          cTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCAGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTa                                                                                                                   >  1:641026/1‑289 (MQ=255)
                                                                                                                                                                           tGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTAt                                                                                                                  <  2:871040/289‑1 (MQ=255)
                                                                                                                                                                             ccGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                                                                                                                                                       <  1:871932/250‑1 (MQ=255)
                                                                                                                                                                             ccGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                                                                                                                                                       >  2:871932/1‑250 (MQ=255)
                                                                                                                                                                                  tGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 >  2:488981/1‑219 (MQ=255)
                                                                                                                                                                                  tGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 <  1:488981/219‑1 (MQ=255)
                                                                                                                                                                                                gcTTATCCGGCCGACGGGTGGCGCTAGATTTGGTAGGCCAGGTCATAGATG‑TTAGCGTCGCATCAGGCATCTGCGCACGACGGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGCGCGGTAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCAc                                                                                              <  1:715221/288‑1 (MQ=21)
                                                                                                                                                                                                 cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCAc                                                                                              <  2:66480/288‑1 (MQ=255)
                                                                                                                                                                                                      ccGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCaa                                                                                        <  2:175199/289‑1 (MQ=255)
                                                                                                                                                                                                         gCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTGGCGCGGGAATTTGTAGGCCTGATTAGACGGGTTATCGCCCCATCCGGCATCTGCACACGACTGCCGGATGCGATTAACGTCTTGTCCGGCCTACATTTCGCCCTGCGGCAGTCATTTCTCATTTCTGATTTCGCGGTTGGCTTTTTTCCATTTGCTATTCACCTTCACTCGCAACgt                                                                                     >  1:30527/1‑288 (MQ=17)
                                                                                                                                                                                                           cTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAgg                                                                                   >  2:896475/1‑289 (MQ=255)
                                                                                                                                                                                                             aCGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGtt                                                                                 <  2:901586/289‑1 (MQ=255)
                                                                                                                                                                                                             aCGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGt                                                                                  <  1:864652/288‑1 (MQ=255)
                                                                                                                                                                                                              cGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TGAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGGCTTGTCCGGCCTACATTTCGCCCGTAGGCAGGCACTAAATAGCTCTGATTACGCGGCTGGCAATGTACCCGTCCCCACTCACCTGCACTCGCaagggcggta                                                                                >  1:464972/1‑281 (MQ=18)
                                                                                                                                                                                                                  tGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCAGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATc                                                                             <  2:641026/288‑1 (MQ=255)
                                                                                                                                                                                                                    gcgtGAGAATTTGTAGGTCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTa                                                                          >  2:380714/1‑289 (MQ=255)
                                                                                                                                                                                                                        gagaATTTGTGGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGAATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGt                                                                       <  2:378983/288‑1 (MQ=255)
                                                                                                                                                                                                                         gtgtataagagaccggagaaggctcg‑TTAGCCTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTACCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAGGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTc                                                                      <  2:404441/264‑1 (MQ=25)
                                                                                                                                                                                                                         agaATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACCTTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCACTGTCTCAGTCCCGATTCACCTGCACTCTCAACGCGGTTACTCTTATGTc                                                                      >  1:494070/1‑288 (MQ=255)
                                                                                                                                                                                                                                     ccTGATAAGACGCT‑TTAGCGTCGCAGCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTAcccc                                                          <  1:830940/288‑1 (MQ=255)
                                                                                                                                                                                                                                     ccTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTAcccc                                                          <  2:591458/288‑1 (MQ=255)
                                                                                                                                                                                                                                         aTAAGACGCG‑TTAGCGTCGCAGCAGGCATCGGCGCACGACTGCCGGATGCGGCGTAAAAGCCTTAGCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTtcat                                                     <  2:126074/289‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCTGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATCAACGCCATGTCCGGCCTCCTTTTCGCTAGCTGGCTGTCATTAACTACTTCTGATTACGCGGCTGGACATGTATCAGTCGCGATTCACCGGCACTCGCACCGAGGATCAACTTATCTCCAATCggctgcctgcatattatcat                                             >  1:404441/1‑268 (MQ=18)
                                                                                                                                                                                                                                                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCTGCCTACATTTCGCCCGTAGGCAGTCATTAAAAAGTTCTTATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCCCTCGCAACCAGGTTCATCTTATGTCGACTCGTTCCCCTTCATCCTCTTcct                                            >  1:739305/1‑286 (MQ=255)
                                                                                                                                                                                                                                                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCAt                                            >  2:185573/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                cAGGCATCTGCGAACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCt                               <  2:409093/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ggCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTg                              >  2:900305/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                       cTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCAt                         >  1:497889/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                        tGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCg                       >  2:621873/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                            cACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCa                                                                                                                                                        >  1:388634/1‑156 (MQ=255)
                                                                                                                                                                                                                                                                            cACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCa                                                                                                                                                        <  2:388634/156‑1 (MQ=255)
                                                                                                                                                                                                                                                                             aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCtt                  >  2:12369/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                             aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCtt                  >  1:675276/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                              cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                                                                                                                         >  1:136999/1‑153 (MQ=255)
                                                                                                                                                                                                                                                                              cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                                                                                                                         <  2:136999/153‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        aTGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTgg        <  2:751004/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGt       <  1:900305/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGt       <  2:344904/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGt       <  2:359733/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGTc      <  2:739305/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              cgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGTCGctc  <  1:380714/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              cgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGTCGc    >  2:104968/1‑286 (MQ=255)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                            
CTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAC‑GCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGTCGCTC  >  NC_000913/374960‑375531

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