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

CTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGTCGCTC  >  NC_000913/374960‑375531
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                            
ttACGGAGGGCGTGAGAATTTGTAGGTCGGGTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGGCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTa                                                                                                                                                                                                                                                                                              <  1:114669/287‑1 (MQ=255)
         gCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGGCTGCTCAGACGCGTTAGCGTCGCATCGGGCATCTGCGCACGACTGCCGGCTGCGGCGTCAACGCCTTaa                                                                                                                                                                                                                                                                                    >  2:622248/1‑288 (MQ=255)
                aaTTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGCTGGCGCGGGAATTTGTCGGCCTGCTAAGACACGTTAGCGTGCCATCAGGCATCTGCGCACGCCTGCCCGCTGCGGCGTaagctacgtattcggct                                                                                                                                                                                                                                                                              >  1:202813/1‑273 (MQ=11)
                          tcagctAAGGCGCGTGAGCGTCGCAGCAGGCGTTTGCGCACGACTGCCGGAGGCGGCGGGAGGGTCTTCTCCGGCATACGGATGGCGCGAGAATTTTTGGGCCTGATAAGAGGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGGGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                                                                                                                                                                   <  1:809867/284‑1 (MQ=14)
                               taggACGCGTTAGAGACGCACAAGGCATGTGGGCACGACTGCCGGACGCGGCGTGACCGGCTTATCCGGCGTACGGATGGCGCGGGAAATTGTAGGACTGATAAGAGGCGTTAGCGTCGCATCAGGCATCTGCGGACGACTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACGGGTGGCGGGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              <  1:147681/286‑1 (MQ=14)
                                          aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGGGCGGGCATTTGTAGGc                                                                                                                                                                                                                                                    >  2:594785/1‑288 (MQ=255)
                                                        cATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                                                                                                                                                                       <  2:333590/287‑1 (MQ=255)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGATGGACCGGGAATTTGTAGGCCTGATaacatacgtacgcg                                                                                                                                                                                                                                 >  1:84440/1‑277 (MQ=255)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              <  1:889671/260‑1 (MQ=255)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              >  2:889671/1‑260 (MQ=255)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                                                                                                                                                                                 >  1:627605/1‑289 (MQ=40)
                                                            tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGCTGGCGCGAGAATTTGTAGGCCTGATAAGAAGCGTTAGCGTCGCATCATGCATCTGCGAACGACTGCCGGACGCGGCGTACACGCTCTATCCGGCCTACGGACCGCCCGGCCATTTACAGGCATGAtcatactacgttgtc                                                                                                                                                                                                                                 >  1:126074/1‑275 (MQ=11)
                                                                 aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCa                                                                                                                                                                                                                            >  2:864652/1‑289 (MQ=39)
                                                                  cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACTTCGCCCGTAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCATGGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTAcactt                                                                                                                                                                                                                                                                       >  2:78999/1‑240 (MQ=11)
                                                                  cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACTTCGCCCGTAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCATGGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTAcactt                                                                                                                                                                                                                                                                       <  1:78999/244‑5 (MQ=11)
                                                                       gcgggtgggggcgGGAAGGCCTTGTGCGGCGTACGGATGGCGGGGGAATTTGTAGGCCGGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAgg                                                                                                                                                                                                                       <  1:197416/279‑1 (MQ=11)
                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                                                                                                                                                                                                      <  2:627605/289‑1 (MQ=39)
                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTCGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                                                                                                                                                                                                      <  2:408581/288‑1 (MQ=21)
                                                                               cggcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAACTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGATTGCGCGGGAATTTGTAGGCCTCATAAGACGGGTTAGCGTCCCATCAGGCATCTGca                                                                                                                                                                                                               >  1:626808/1‑288 (MQ=255)
                                                                                              aTCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                                                                                                                                                <  2:316130/256‑1 (MQ=40)
                                                                                              aTCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                                                                                                                                                >  1:316130/1‑256 (MQ=40)
                                                                                                   gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              <  1:801647/221‑1 (MQ=255)
                                                                                                   gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                                                                                                              >  2:801647/1‑221 (MQ=255)
                                                                                                   gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCCGCCCACGACTGCCGGGTGCGGTaa                                                                                                                                                                                          >  1:425315/1‑289 (MQ=255)
                                                                                                     cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAAc                                                                                                                                                                                        <  2:780792/289‑1 (MQ=255)
                                                                                                       aCGGATGGCGCGGGAATTTGTAGGCCTGAGAAGACGCGTTGGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                                                                                                                                                                                      <  2:352049/289‑1 (MQ=255)
                                                                                                            tGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 <  2:16053/289‑1 (MQ=255)
                                                                                                              tggcggacaTTTGTAGGCTGTGTGAGACGCGTTGCAGACGCATCGTGCACCGGCGCACGTCTGCCGGCTGCGGGTGGACGGCCTTTGCCGGCCTAAGGTTGGCGGGAGAATTTGTAGGCCTGATAAGACGCGTTGGCGTCGCATCGGGCAGCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                                                                                                                                                <  1:622248/280‑1 (MQ=14)
                                                                                                              gggcgGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATa                                                                                                                                                                                                                                              <  1:570506/224‑1 (MQ=255)
                                                                                                              gcgcGGTAATGGGGAGGCGGGGTAGGACGGGGTAGGGGCGGATGAGGCATCGGCGCACGACTGCCGGGTGGGGCGGGAGCGCCTTATCCGGCCGACGGGTGGCGCGAGGATTTGTAGGGCTGAGAAGAGGCGTTAGCGTCGCATCGGGCGTCTGAGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                                                                                                                                                <  1:805332/288‑1 (MQ=14)
                                                                                                              gcgcGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCCTCAGGCATCTGCGCACGACTGCCGGATGCCGCGGAAAATCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGGGTTGACGTCGCACGCGGCCTCTGCACCCGGCCCCCGGATCCGATAAACggcatttc                                                                                                                                                                                >  1:352049/1‑281 (MQ=14)
                                                                                                              gcgcGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATa                                                                                                                                                                                                                                              >  2:570506/1‑226 (MQ=255)
                                                                                                                                 tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAACCGCCTTATCCGGTCTACGGATGGCGCGGGAATTTGTAGGCCTGTTAAGACGCGTTCGCGTCCCATCAGGCATCTGTACACCGCTTCCGTCTGCGATAAGCGTCTTGTGCGGCCtcactttccctcgt                                                                                                                                                             >  1:378983/1‑275 (MQ=14)
                                                                                                                                  gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTG‑‑AAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 >  1:488315/1‑265 (MQ=255)
                                                                                                                                  gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTG‑‑AAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 <  2:488315/265‑1 (MQ=255)
                                                                                                                                    tAAGAGGCGTTAGCGTGGCATCAGGCAGCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                                                                                                                         <  1:733655/289‑1 (MQ=255)
                                                                                                                                    tAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                                                                                                                         <  2:84440/289‑1 (MQ=255)
                                                                                                                                         cgcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                                                                                                                                                    <  2:707577/289‑1 (MQ=255)
                                                                                                                                          gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATACGACGCGTTAGCGTCCCCTCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCCGTcct                                                                                                                                                   >  2:383352/1‑287 (MQ=255)
                                                                                                                                          aacgtgcattttgtatgcccatcccagcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATACGACGCGTTAGCGTCGCATCAGGCCTCTGCGCACGACTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGATGGCGCGGGCATTTTTAGGCCTGCTACGACGCGTTACCGTCCCATCCTGCACATGCCAACGACTTCCGGATGCGATAACCTGCCTGTCCGGCCTACATTTCTCCCCTAGGCAGTcct                                                                                                                                                   >  1:793117/27‑287 (MQ=11)
                                                                                                                                            gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTa                                                                                                                                                                                                                                                                                              <  2:836658/148‑1 (MQ=35)
                                                                                                                                              tAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                                                                                                                                               <  2:416394/289‑1 (MQ=255)
                                                                                                                                                     gCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCt                                                                                                                                        >  2:400260/1‑289 (MQ=255)
                                                                                                                                                      cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                                                                                                    >  1:836658/1‑148 (MQ=37)
                                                                                                                                                             cagcTGCGAACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCATACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcg                                                                                                                                 <  1:923602/285‑1 (MQ=255)
                                                                                                                                                                 ttaccACGACGGCCGTATGCGGCGTGCACGCCTTGTCCGGCCTACGAGTGGCGCGAGAATTTGTAGGCCTGATAGGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc                                                                                                                             <  2:505115/284‑1 (MQ=255)
                                                                                                                                                                 cgggcggaccgtgcgtaTGCGGCGCGACCGCCTGATCCGGCATACGGTGGGAGGGATAATGTGTAGGCCTGGTAAGACGCGTTAGCGGCGCATCATGCATCGGCGCGCGACGGCCGGATGCGGCGTAATCGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGCACACGACGGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCt                                                                                                                            <  1:512868/273‑1 (MQ=255)
                                                                                                                                                                   cgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCCCGGGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGCCATCTGCACACGACTGCCGGCTGCGATCACCGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCCTTAACTTGTTCTTATTTcgcc                                                                                                                               >  1:751004/1‑283 (MQ=17)
                                                                                                                                                                    gcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTgg                                                                                                                          <  2:425315/288‑1 (MQ=255)
                                                                                                                                                                    acACGACTGCCGGATGGGGCGTGAACGCCTTATCCGGCCTACGGGGGGCGCGAGAAGTTGTAGGCCTGATAAGACGCGTTAGCGTGGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTgg                                                                                                                          <  1:886655/287‑1 (MQ=255)
                                                                                                                                                                      aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGCTAAGACGCCTTAGCGGCCCATCAGGCATCTGCCCCCGACTCCCGCATGCACTAAACGGCTTGTTCGGCCTACCTTTCCCCCGTCGGCAGTCATCAATCAGCTCTTCTTACCCGTCTGCCac                                                                                                                       >  1:665419/1‑288 (MQ=14)
                                                                                                                                                                      aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCCCCCGTAGGCCGCCATTAAATAGTTCTTATTACGCGCCGGGcca                                                                                                                       >  1:409093/1‑287 (MQ=255)
                                                                                                                                                                       cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACACTTCGCCCATAAACCGTAGGCCTGATCAGACGCGTTAGCGTCGCATCAGGCATCTCCGCACGAATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCCTACGCAAACGACGCACAGATGCGATAAACGTCCTGTCAGGCCGAGAATCCCCCCGGCGGCAGCCACGAAACACTTCTGCTTATGCTGCTGGCaact                                                                                                                     >  1:666201/1‑287 (MQ=9)
                                                                                                                                                                          cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCAGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTa                                                                                                                   >  1:641026/1‑289 (MQ=255)
                                                                                                                                                                           tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTAt                                                                                                                  <  2:871040/289‑1 (MQ=255)
                                                                                                                                                                             ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                                                                                                                                                       <  1:871932/250‑1 (MQ=255)
                                                                                                                                                                             ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                                                                                                                                                       >  2:871932/1‑250 (MQ=255)
                                                                                                                                                                                  tGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 >  2:488981/1‑219 (MQ=255)
                                                                                                                                                                                  tGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                                                                                                                 <  1:488981/219‑1 (MQ=255)
                                                                                                                                                                                                gcTTATCCGGCCGACGGGTGGCGCTAGATTTGGTAGGCCAGGTCATAGATGTTAGCGTCGCATCAGGCATCTGCGCACGACGGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGCGCGGTAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCAc                                                                                              <  1:715221/288‑1 (MQ=21)
                                                                                                                                                                                                 cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCAc                                                                                              <  2:66480/288‑1 (MQ=255)
                                                                                                                                                                                                      ccGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCaa                                                                                        <  2:175199/289‑1 (MQ=255)
                                                                                                                                                                                                         gCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTGGCGCGGGAATTTGTAGGCCTGATTAGACGGGTTATCGCCCCATCCGGCATCTGCACACGACTGCCGGATGCGATTAACGTCTTGTCCGGCCTACATTTCGCCCTGCGGCAGTCATTTCTCATTTCTGATTTCGCGGTTGGCTTTTTTCCATTTGCTATTCACCTTCACTCGCAACgt                                                                                     >  1:30527/1‑288 (MQ=17)
                                                                                                                                                                                                           cTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAgg                                                                                   >  2:896475/1‑289 (MQ=255)
                                                                                                                                                                                                             aCGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGtt                                                                                 <  2:901586/289‑1 (MQ=255)
                                                                                                                                                                                                             aCGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGt                                                                                  <  1:864652/288‑1 (MQ=255)
                                                                                                                                                                                                              cGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTGAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGGCTTGTCCGGCCTACATTTCGCCCGTAGGCAGGCACTAAATAGCTCTGATTACGCGGCTGGCAATGTACCCGTCCCCACTCACCTGCACTCGCaagggcggta                                                                                >  1:464972/1‑281 (MQ=18)
                                                                                                                                                                                                                  tGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCAGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATc                                                                             <  2:641026/288‑1 (MQ=255)
                                                                                                                                                                                                                    gcgtGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTa                                                                          >  2:380714/1‑289 (MQ=255)
                                                                                                                                                                                                                        gagaATTTGTGGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGAATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGt                                                                       <  2:378983/288‑1 (MQ=255)
                                                                                                                                                                                                                         gtgtataagagaccggagaaggctcgTTAGCCTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTACCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAGGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTc                                                                      <  2:404441/264‑1 (MQ=25)
                                                                                                                                                                                                                         agaATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACCTTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCACTGTCTCAGTCCCGATTCACCTGCACTCTCAACGCGGTTACTCTTATGTc                                                                      >  1:494070/1‑288 (MQ=255)
                                                                                                                                                                                                                                     ccTGATAAGACGCTTTAGCGTCGCAGCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTAcccc                                                          <  1:830940/288‑1 (MQ=255)
                                                                                                                                                                                                                                     ccTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTAcccc                                                          <  2:591458/288‑1 (MQ=255)
                                                                                                                                                                                                                                         aTAAGACGCGTTAGCGTCGCAGCAGGCATCGGCGCACGACTGCCGGATGCGGCGTAAAAGCCTTAGCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTtcat                                                     <  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)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                            
CTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGTCGACTCGTACCCCTTCATCATCTTCATCCCGCCTGATGCTGACCATCGGGCTTTGTTTTTTGGTCGCTC  >  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: