Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,176 A→G 8.3% intergenic (+295/‑283) 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,1760AG 8.3% 166.8 / 7.0 61intergenic (+295/‑283)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (22/31);  new base G (2/3);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

GCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTA  >  NC_000913/374890‑375459
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                              
gCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCgaga                                                                                                                                                                                                                                                                                              <  2:81992/289‑1 (MQ=255)
           gCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGc                                                                                                                                                                                                                                                                                    >  1:413467/1‑288 (MQ=255)
             aTTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCt                                                                                                                                                                                                                                                                                  <  1:479175/288‑1 (MQ=255)
                       cccATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGAcgcg                                                                                                                                                                                                                                                                       >  1:284637/1‑289 (MQ=255)
                                   cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAt                                                                                                                                                                                                                                                           <  2:413467/289‑1 (MQ=39)
                                     aCTGCCGGGTGCGGCGGGAAGGCCTTATCCGGCCTACGGAGGGCGGGAGAATTTGTAGGTCTGATACGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATc                                                                                                                                                                                                                                                          <  1:130853/288‑1 (MQ=255)
                                      cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCCTGAACGCCTTATCCGGCCTTCGGGTTGCGCGCGCATTTTTCGGCCTGATAAGACCCGCTAGCGTCGCATccg                                                                                                                                                                                                                                                        >  1:7233/1‑287 (MQ=11)
                                                gggcgTGACCGCCTTATCCGGCCTAGGGAGGGCGTGGGAATTTGTAGGTCTGATAAGACGCGGTAGCGTCGCATCAGGCATCGGCGCACGACTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTgcg                                                                                                                                                                                                                                               <  1:78413/287‑1 (MQ=255)
                                                   cgTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACg                                                                                                                                                                                                                                           >  2:288353/1‑289 (MQ=40)
                                                                              ggCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGCTACGACGCG‑TTAGCGTCGCATCAGGCAGATGCGCACGACTGCCGAGTACGGCCTGAACGCCTTATCGGGACCACGGGTGGAGAGCGAATTCGTAGGACTGATAAGACGCGTTCGAGACGCATCAGGACTGTGCGCACGACTGCCGGATGTGGCGTAACAGCCTTa                                                                                                                                                                                                                >  1:345108/1‑289 (MQ=14)
                                                                                cgcgagaATATAGAGGCCTGATGAGACGCGGTGCCGTCGCAGCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                               <  1:83535/285‑1 (MQ=14)
                                                                                                           gcgTTAGGGTCGCATACGGGATGTGCGAACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg                                                                                                                                                                                    <  1:288353/288‑1 (MQ=255)
                                                                                                              ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGCGCATGCCTTTTCCCGCCTACGGCTGGCGCGCGGACTTGTCGCCCTGTTCCGCCGCGTTAGCGTCCCATCTGGCCTCTGTCCAACCCTGCCGGACGCCGCCTAACCCCCTTACCCGGCCCTCCGCGGGCGTGGGCATTTGTat                                                                                                                                                                                 >  1:64412/1‑287 (MQ=11)
                                                                                                               tggagTCGCATCAGCGATCTGCGCTGGGCTGCCAGAGTCGGCGTGAAGGCCTTATCCGGCTTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGc                                                                                                                                                                               <  1:398636/285‑1 (MQ=14)
                                                                                                               tggCGTGGCACCAGCGATCGGCGGACGACTGCCGGATGCGGCGTGAACGCCTTATGCGGCCTAAGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGc                                                                                                                                                                               <  1:286318/287‑1 (MQ=255)
                                                                                                                    tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                                                                                                          >  2:487030/1‑289 (MQ=255)
                                                                                                                         tgatgctgcTGCGCAGGGCGGCCGGATGAGGCGTGAAGGCGGGAGCGGGCCTACGGCTGGCGCGGGAAGTGGTAGGCCGGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGa                                                                                                                                                                      <  1:106020/280‑1 (MQ=14)
                                                                                                                            ggCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATa                                                                                                                                                                                                                                                                              <  1:323187/181‑1 (MQ=39)
                                                                                                                                      cACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTa                                                                                                                                                                                                                                                                    >  2:323187/1‑181 (MQ=39)
                                                                                                                                            tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAgg                                                                                                                                                 <  2:45945/289‑1 (MQ=39)
                                                                                                                                              ccGGAGGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGc                                                                                                                                                <  2:269136/288‑1 (MQ=39)
                                                                                                                                               cggagtgcgcgTGAACGCCTGAGCCGGCGTACGGAGGGAGCGGGAAGTTGTAGGACTGCTAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGGACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAt                                                                                                                                              <  1:388931/281‑1 (MQ=11)
                                                                                                                                                     cggcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGcac                                                                                                                                        >  2:79602/1‑289 (MQ=40)
                                                                                                                                                            aaCGCCTTATCCGGCCTACACTTCGC‑CCGTAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACACTTCGC‑CCATACACCGCAGGCCTGCTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAGCGGTGCCGCCGGCATTTCTCGGCCTCCTAAGACCCG‑TTAGCGTCGCCTGCGCCATTTGCCCACAACTGc                                                                                                                                >  1:360942/1‑288 (MQ=9)
                                                                                                                                                                  ttATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGCAGGCCTGATAAGCCGCG‑TTAGCGTCGCATCAGACATCTTCACCAGACTGctc                                                                                                                              >  1:305479/1‑284 (MQ=255)
                                                                                                                                                                          ccTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGcaca                                                                                                                                       >  1:189102/1‑269 (MQ=255)
                                                                                                                                                                          ccTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGcaca                                                                                                                                       <  2:189102/269‑1 (MQ=255)
                                                                                                                                                                             aCGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACg                                                                                                               <  1:487030/288‑1 (MQ=255)
                                                                                                                                                                               ggATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATCGCGCGGGAATTTGTCGGCCTGATAAGACGCG‑TTAGCGTCGCATCCGGCATCTGCACCCGACTTCCGGAGGCG‑‑ATAAGCGTCt                                                                                                            >  1:269911/1‑289 (MQ=255)
                                                                                                                                                                                    gcgcGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCAGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGAATCAGGCATCTGCGCACGACCGCCGCATGCGGCGTAAACACCTTATCCGGCCTACGGATCGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGTCATCTGCACACAGCTGCCGGGTCTG‑‑CTAAAAGTCTTGTc                                                                                                        >  1:438283/1‑288 (MQ=14)
                                                                                                                                                                                               tGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACAttt                                                                                            <  2:4188/289‑1 (MQ=255)
                                                                                                                                                                                               tGTAGGCCTGATAAAACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACAttt                                                                                            <  2:46706/289‑1 (MQ=255)
                                                                                                                                                                                                      cTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAAACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGCATGCCGGAATCG‑‑ATAAACGGCCTGTTCGGCCTACACTTCCCCCGt                                                                                     >  1:46706/1‑289 (MQ=14)
                                                                                                                                                                                                        gATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGt                                                                                                                                                          <  2:120327/220‑1 (MQ=38)
                                                                                                                                                                                                          tAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                                                 <  2:284637/289‑1 (MQ=255)
                                                                                                                                                                                                           aaGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg                                                                                                                                                                                    <  1:216417/192‑1 (MQ=38)
                                                                                                                                                                                                               gccg‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                                                                            <  2:12937/287‑1 (MQ=255)
                                                                                                                                                                                                               cgcg‑TTGTCGTCGCATCAGGCAGCTGCCCACGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCGACGGGGGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATTCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                                                                            <  2:345108/289‑1 (MQ=255)
                                                                                                                                                                                                                  g‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGc                                                                                                                                                >  1:120327/1‑220 (MQ=38)
                                                                                                                                                                                                                     tAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                                                                       <  1:79602/289‑1 (MQ=255)
                                                                                                                                                                                                                      aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                                                                                                          >  2:216417/1‑192 (MQ=38)
                                                                                                                                                                                                                         ttCGCATCAGGGATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAg                                                                    <  1:194796/287‑1 (MQ=255)
                                                                                                                                                                                                                         gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGt                                                                   >  2:194796/1‑289 (MQ=255)
                                                                                                                                                                                                                          tCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACACTTCGC‑CCATAAACCGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATaa                                                                                                                                                                        >  1:313640/1‑188 (MQ=11)
                                                                                                                                                                                                                          tCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACACTTCGC‑CCATAAACCGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATaa                                                                                                                                                                        <  2:313640/188‑1 (MQ=11)
                                                                                                                                                                                                                           cGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGt                                                                   >  2:106262/1‑287 (MQ=255)
                                                                                                                                                                                                                            gctcaaGGAATCTGGGGACGACGGCGGGAGGCGGCGTGAACGCCTTGTCGGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACTTCGCCCGTAAGCCGtaggcctgataa                                                                   <  1:328391/283‑12 (MQ=11)
                                                                                                                                                                                                                             cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTg                                                               >  2:392905/1‑289 (MQ=255)
                                                                                                                                                                                                                                        tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcgcg                                                       >  2:423541/1‑286 (MQ=255)
                                                                                                                                                                                                                                           gcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGc                                                 <  1:106262/289‑1 (MQ=255)
                                                                                                                                                                                                                                           gcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTgg                                                  <  2:125788/288‑1 (MQ=255)
                                                                                                                                                                                                                                               gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCC‑‑GTAAACGTCTTGTCCTGCCTACATTTCGCCCGTAGGCAGCCCTTAAATAGTTCTGATTATCGGGCTGGCAATg                                             >  1:125788/1‑289 (MQ=17)
                                                                                                                                                                                                                                                                   aCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACAtccggcagtcg                                                                                    <  1:478153/230‑11 (MQ=255)
                                                                                                                                                                                                                                                                   aCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACAtccggcagtcg                                                                                    >  2:478153/1‑220 (MQ=255)
                                                                                                                                                                                                                                                                       cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACt                     <  1:392905/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ttATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGt                                                                                                              >  2:418940/1‑199 (MQ=255)
                                                                                                                                                                                                                                                                        ttATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGt                                                                                                              <  1:418940/199‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  tACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGt          <  1:423541/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  caCGGGTGGCGCGAGACTTTTTAGGCCTGTTAAGACGCGTTAGCGTCGCATCAGGAATCGGGGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGt          <  1:219575/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          gcgcGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACTTCGC‑CCGTAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACACTTCGCCCATaaaccgtaggcctgataa                                                                   >  1:7495/1‑208 (MQ=14)
                                                                                                                                                                                                                                                                                          gcgcGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACTTCGC‑CCGTAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACACTTCGCCCATaaaccgtaggcctgataa                                                                   <  2:7495/225‑18 (MQ=14)
                                                                                                                                                                                                                                                                                          gcgcGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTa  >  2:329189/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                              
GCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCG‑‑ATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTA  >  NC_000913/374890‑375459

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: ‑