Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 375,168 G→A intergenic (+287/‑291) mhpE → / → mhpT 4‑hydroxy‑2‑oxovalerate aldolase/3‑hydroxyphenylpropionate/3‑ hydroxycinnamate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913375,1680GA98.2% 158.0 / ‑4.0 56intergenic (+287/‑291)mhpE/mhpT4‑hydroxy‑2‑oxovalerate aldolase/3‑hydroxyphenylpropionate/3‑ hydroxycinnamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/24);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

TGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGT  >  NC_000913/374885‑375462
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                      
tgaCTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGGGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAAGGATGGCGTGGGGATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACGGATGGCGCGAGAATTTGt                                                                                                                                                                                                                                                                                        >  2:369153/1‑300 (MQ=255)
   cTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAgg                                                                                                                                                                                                                                                                                     <  2:583864/300‑1 (MQ=255)
                        tATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                                                                                                                                                                                <  1:692537/300‑1 (MQ=255)
                           gCCCATCCCATCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGc                                                                                                                                                                                                                                                             >  1:401061/1‑300 (MQ=255)
                                   cAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCa                                                                                                                                                                                                                                                     >  1:325686/1‑300 (MQ=255)
                                          aCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATTTGCGCACGACTGCCGGATGCGGCGTGAACGTCTTATCTGGCCTACGGATAGCGTGGGAATTTGTCAGCCCGCTAAGACGCAGTGGCGTACCATCCAGCAGTCGCGCACAACTGCCGGATGCGGCGGGACCGACTTCCCCGACATCCGGATGGAGCAAGCATCTGTAGGCCTGATCGGACGCACTAGCGTCGCATCAGGCATCTgagc                                                                                                                                                                                                                                              >  2:387872/1‑297 (MQ=2)
                                             gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACg                                                                                                                                                                                                                                           >  1:340673/1‑300 (MQ=255)
                                               cGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAc                                                                                                                                                                                                                                         <  1:91064/300‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTCGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCCGCGCACGACTGCCGGCTGCGGCGt                                                                                                                                                                                                                          >  1:371046/1‑300 (MQ=255)
                                                                ccTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTaa                                                                                                                                                                                                                        >  1:4104/1‑300 (MQ=255)
                                                                        ggCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                                                                                                >  1:83263/1‑300 (MQ=255)
                                                                         gCCTACGGATGGCGTGGGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGTGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGCGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGACTGATAAGACGAGTTAGCGTCTCCTCGGGCATCTGCGCACGACTACAGGATATGGCGTAAACGCCTTAt                                                                                                                                                                                                               >  2:216298/1‑300 (MQ=14)
                                                                         gCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATgcgataaacgtcttgtcc                                                                                                                                                                                                               >  1:278828/1‑285 (MQ=9)
                                                                          ccTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                                                                                              <  1:235581/300‑1 (MQ=255)
                                                                              cGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCATCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTCGGCATGACAAGACGCGTTAGCATCGCATCCGGCATCTGCGCACGACTGCCGGATGCGACGCGAACGCCTTATCGAGCCTACGGATGACGCGAGAATCTGTAGGCATGATAAGACGCGTTAGCGTCGAATACGGCATCTGCACACGGCTGCCGGATgcggctgaaataccttctccggc                                                                                                                                                                                                          >  2:408848/1‑282 (MQ=11)
                                                                               ggATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                                                                                         <  1:705944/300‑1 (MQ=255)
                                                                                             tttGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATACAGCGTAAACGCCTTATCCGGCCTACGGATGGCGCgg                                                                                                                                                                                           >  1:251202/1‑300 (MQ=255)
                                                                                                    gTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGTGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTCAGCGTCGAATCAGGCATCTGCGCCCGACTGACGGATGCGGCGTAAATGCATTATTCGGCCTACGGATGgcgtaggaacttg                                                                                                                                                                                    >  2:229336/1‑290 (MQ=14)
                                                                                                    gTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGGCGCGTTAGCGTCGTATCAGGCATCTGCGCATGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg                                                                                                                                                                                    >  2:152368/1‑300 (MQ=255)
                                                                                                          tAAGACGCTTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGcc                                                                                                                                                                              <  1:109797/300‑1 (MQ=255)
                                                                                                          tAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGcc                                                                                                                                                                              <  1:609132/300‑1 (MQ=255)
                                                                                                                       cGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAATGCCTTCCCCGGCCTACGGATAGCGCAAGAATTTGTAGGCCTGATAAGACGTGTTAGCGTCGCATCGAGCATCTGCGCACACCTGCCGGATGCGGCGCAAACGCCTTATATGGACTACGGATGGCGCGGGCATTTGTAGACATGATAAGAcgagc                                                                                                                                                                 >  2:85025/1‑297 (MQ=14)
                                                                                                                            cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATCAGACGCGTTAGCg                                                                                                                                                            >  1:323940/1‑300 (MQ=255)
                                                                                                                                        gcgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGGATGGCGCGGGGATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCCTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGACTACGGATGGCGCGGGACTTTGTAGGGCTGATAAGAAGCGTTAGCGTCGCATTAGGCa                                                                                                                                                >  2:241875/1‑300 (MQ=11)
                                                                                                                                        gcgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGGGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCNGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGACTTCTCCGGCCTACGGATGACGCGGGAATTTGTAGGCCTGATAAGACGCATTCGCGTCGCATCGGGCa                                                                                                                                                >  2:664350/1‑300 (MQ=14)
                                                                                                                                          gcgccgcttccGGACGCGGCGTGAACGCCTTATCCGGCCTGCGGATGGCGCGGGAATTTGTAGGCCGGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATc                                                                                                                                              <  2:82960/291‑1 (MQ=255)
                                                                                                                                                         gcgtcgTGAACGCCTTGTCCGGCCTACGTATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCg                                                                                                                               <  2:738473/296‑1 (MQ=255)
                                                                                                                                                              aggAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCATGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTACCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGc                                                                                                                          <  2:395661/298‑1 (MQ=255)
                                                                                                                                                                     ccTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGGCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAAc                                                                                                                   <  1:243578/300‑1 (MQ=255)
                                                                                                                                                                         aTCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCt                                                                                                               >  1:279648/1‑300 (MQ=255)
                                                                                                                                                                             ggCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                                                                           >  1:52996/1‑300 (MQ=255)
                                                                                                                                                                                       tGGCGCGGGAATTTTTATGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACTGATGGCACGAGAATTTGTATGCCTGATAAGACGCGTTAGCGTCACATCATGCATCTGCGCACGACGGCAGGAGGCGGCGTAAACGCCTTATCCTGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAt                                                                                                 <  1:333750/300‑1 (MQ=17)
                                                                                                                                                                                       tGGCGCGGGAATTTTTAGGCCTGATAAGACGCGTTAGCGTCGCATCATGCATCGGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATTGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAt                                                                                                 <  1:165351/300‑1 (MQ=255)
                                                                                                                                                                                        tgcgcGGGAATCTGTATGTCTGACAAGACTCATCAGCGCCGCATCAGGCATCGTCGCACGACTGCCGGATGTGGCGTGAACGCCTGATCCGACCTAAGGAGCGCGAGAGAATTTGTGGGACTGAGAAGGCGCGTTAGCGTCGCATCAGGCATCGGCGCACGACTGCTGGATGCGGCGGAAACGCCTTATCCGGCCTACGAATGGCGCTGGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAtt                                                                                                <  2:336458/299‑1 (MQ=14)
                                                                                                                                                                                                    tGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                                                    <  1:110816/300‑1 (MQ=255)
                                                                                                                                                                                                               gaaGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAAt                                                                         <  1:108539/299‑1 (MQ=255)
                                                                                                                                                                                                                                gCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTCCGCGGc                                                        >  1:266157/1‑300 (MQ=255)
                                                                                                                                                                                                                                gCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGACCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc                                                        >  1:269624/1‑300 (MQ=255)
                                                                                                                                                                                                                                 cATCAGGCATCTGCGCATGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGAAGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGATGACGCGGGAATTTGTAGGCATGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGAATGCCGGATGCGAGAAACGTCTTGTCCAGCCTACATTTCGACCGTAGGCAGTAATTAAATAGTTCCGACTACGCGGCt                                                       >  2:318495/1‑300 (MQ=17)
                                                                                                                                                                                                                                           tttgagcgcgattGCCGGATGCTGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGGAGTCCTGATAAGACGCGTTAGCGTCGCATCACGCATCGGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCGACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCTTCGCATCAGGCATCTGCACACGACTGCCTGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATgga                                              <  2:33643/288‑3 (MQ=18)
                                                                                                                                                                                                                                            cTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGGATGGCGCAAGAATTTGTAGGCCTGATAAGACTCTTTATCGTTGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCTTTATCCGTCCTACGGATTGCGCGGGAATTTGTGGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTAt                                             <  2:358406/300‑1 (MQ=18)
                                                                                                                                                                                                                                             tacgcACGACGGCCGGGTGCAGCTGGGACGCTTTATCCGACCTACTGATGGCGCGAGCATTTGTAGGTCGGATAAGACGCGTTATCGTCGGGTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCTTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATc                                            <  2:387800/298‑1 (MQ=18)
                                                                                                                                                                                                                                                   cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGACCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTcgc                                      >  2:310597/1‑300 (MQ=255)
                                                                                                                                                                                                                                                   cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCAGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTcgc                                      >  2:210116/1‑300 (MQ=255)
                                                                                                                                                                                                                                                     aCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGa                                    >  1:675151/1‑300 (MQ=255)
                                                                                                                                                                                                                                                             aTGCGGTGTGGACGCCTTGTCTGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGAGAAGACGCGTTATCGTCGTATCAGGCATCTGCGCACGACTGCCGGATGCGGCGGAAACGCCTTATTCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCAATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTg                            <  2:26173/300‑1 (MQ=21)
                                                                                                                                                                                                                                                             aTGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTg                            <  1:581906/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                cggcgTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCTTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCAc                         <  1:14222/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gcgTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTc                       >  1:567136/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                         cGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGNTGCGATAAACGTCTTGTCCGGCATACATTTCGCACGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCCACGa                >  1:541905/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                           ccTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAgg              >  1:683792/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                tCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGGAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCAt         <  1:576840/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   tgCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGTGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCNGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTATCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCtt      <  2:178967/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   ggcgtACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACATGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTCGGCAGTCGTTCAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACGCGCAACGAAGTTCATCtt      >  2:322145/5‑300 (MQ=21)
                                                                                                                                                                                                                                                                                    gCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTa     >  1:575409/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                       gaCGGATGGCGCGAGAATTTTGGGGTCTGATAAGACGCGGTATCTTCGCACCAGNCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGt  <  2:573558/299‑1 (MQ=25)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                      
TGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGT  >  NC_000913/374885‑375462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: