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

TGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGT  >  NC_000913/374885‑375462
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                      
tgaCTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGGGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAAGGATGGCGTGGGGATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACGGATGGCGCGAGAATTTGt                                                                                                                                                                                                                                                                                        >  2:369153/1‑300 (MQ=255)
   cTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAgg                                                                                                                                                                                                                                                                                     <  2:583864/300‑1 (MQ=255)
                        tATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                                                                                                                                                                                <  1:692537/300‑1 (MQ=255)
                           gCCCATCCCATCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGc                                                                                                                                                                                                                                                             >  1:401061/1‑300 (MQ=255)
                                   cAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCa                                                                                                                                                                                                                                                     >  1:325686/1‑300 (MQ=255)
                                          aCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATTTGCGCACGACTGCCGGATGCGGCGTGAACGTCTTATCTGGCCTACGGATAGCGTGGGAATTTGTCAGCCCGCTAAGACGCAGTGGCGTACCATCCAGCAGTC‑GCGCACAACTGCCGGATGCGGCGGGACCGACTTCCCCGACATCCGGATGGAGCAAGCATCTGTAGGCCTGATCGGACGCACTAGCGTCGCATCAGGCATCTgagc                                                                                                                                                                                                                                              >  2:387872/1‑297 (MQ=2)
                                             gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACg                                                                                                                                                                                                                                           >  1:340673/1‑300 (MQ=255)
                                               cGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAc                                                                                                                                                                                                                                         <  1:91064/300‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTCGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCCGCGCACGACTGCCGGCTGCGGCGt                                                                                                                                                                                                                          >  1:371046/1‑300 (MQ=255)
                                                                ccTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTaa                                                                                                                                                                                                                        >  1:4104/1‑300 (MQ=255)
                                                                        ggCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                                                                                                >  1:83263/1‑300 (MQ=255)
                                                                         gCCTACGGATGGCGTGGGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGTGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGCGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGACTGATAAGACGAGTTAGCGTCTCCTCGGGCATCTGCGCACGACTACAGGATATGGCGTAAACGCCTTAt                                                                                                                                                                                                               >  2:216298/1‑300 (MQ=14)
                                                                         gCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATgcgataaacgtcttgtcc                                                                                                                                                                                                               >  1:278828/1‑285 (MQ=9)
                                                                          ccTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                                                                                              <  1:235581/300‑1 (MQ=255)
                                                                              cGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCATCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTCGGCATGACAAGACGCGTTAGCATCGCATCCGGCA‑TCTGCGCACGACTGCCGGATGCGACGCGAACGCCTTATCGAGCCTACGGATGACGCGAGAATCTGTAGGCATGATAAGACGCGTTAGCGTCGAATACGGCATCTGCACACGGCTGCCGGATgcggctgaaataccttctccggc                                                                                                                                                                                                          >  2:408848/1‑282 (MQ=11)
                                                                               ggATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                                                                                         <  1:705944/300‑1 (MQ=255)
                                                                                             tttGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATACAGCGTAAACGCCTTATCCGGCCTACGGATGGCGCgg                                                                                                                                                                                           >  1:251202/1‑300 (MQ=255)
                                                                                                    gTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGTGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTCAGCGTCGAATCAGGCATCTGCGCCCGACTGACGGATGCGGCGTAAATGCATTATTCGGCCTACGGATGgcgtaggaacttg                                                                                                                                                                                    >  2:229336/1‑290 (MQ=14)
                                                                                                    gTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGGCGCGTTAGCGTCGTATCAGGCATCTGCGCATGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg                                                                                                                                                                                    >  2:152368/1‑300 (MQ=255)
                                                                                                          tAAGACGCTTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGcc                                                                                                                                                                              <  1:109797/300‑1 (MQ=255)
                                                                                                          tAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGcc                                                                                                                                                                              <  1:609132/300‑1 (MQ=255)
                                                                                                                       cGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAATGCCTTCCCCGGCCTACGGATAGCGCAAGAATTTGTAGGCCTGATAAGACGTGTTAGCGTCGCATCGAGCATCTGCGCACACCTGCCGGATGCGGCGCAAACGCCTTATATGGACTACGGATGGCGCGGGCATTTGTAGACATGATAAGAcgagc                                                                                                                                                                 >  2:85025/1‑297 (MQ=14)
                                                                                                                            cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATCAGACGCGTTAGCg                                                                                                                                                            >  1:323940/1‑300 (MQ=255)
                                                                                                                                        gcgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGGATGGCGCGGGGATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCCTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGACTACGGATGGCGCGGGACTTTGTAGGGCTGATAAGAAGCGTTAGCGTCGCATTAGGCa                                                                                                                                                >  2:241875/1‑300 (MQ=11)
                                                                                                                                        gcgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGGGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCNGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGACTTCTCCGGCCTACGGATGACGCGGGAATTTGTAGGCCTGATAAGACGCATTCGCGTCGCATCGGGCa                                                                                                                                                >  2:664350/1‑300 (MQ=14)
                                                                                                                                          gcgccgcttccGGACGCGGCGTGAACGCCTTATCCGGCCTGCGGATGGCGCGGGAATTTGTAGGCCGGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATc                                                                                                                                              <  2:82960/291‑1 (MQ=255)
                                                                                                                                                         gcgtcgTGAACGCCTTGTCCGGCCTACGTATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCg                                                                                                                               <  2:738473/296‑1 (MQ=255)
                                                                                                                                                              aggAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCATGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTACCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGc                                                                                                                          <  2:395661/298‑1 (MQ=255)
                                                                                                                                                                     ccTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGGCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAAc                                                                                                                   <  1:243578/300‑1 (MQ=255)
                                                                                                                                                                         aTCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCt                                                                                                               >  1:279648/1‑300 (MQ=255)
                                                                                                                                                                             ggCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                                                                           >  1:52996/1‑300 (MQ=255)
                                                                                                                                                                                       tGGCGCGGGAATTTTTATGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACTGATGGCACGAGAATTTGTATGCCTGATAAGACGCGTTAGCGTCACATCATGCATCTGCGCACGACGGCAGGAGGCGGCGTAAACGCCTTATCCTGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAt                                                                                                 <  1:333750/300‑1 (MQ=17)
                                                                                                                                                                                       tGGCGCGGGAATTTTTAGGCCTGATAAGACGCGTTAGCGTCGCATCATGCA‑TCGGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATTGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAt                                                                                                 <  1:165351/300‑1 (MQ=255)
                                                                                                                                                                                        tgcgcGGGAATCTGTATGTCTGACAAGACTCATCAGCGCCGCATCAGGCA‑TCGTCGCACGACTGCCGGATGTGGCGTGAACGCCTGATCCGACCTAAGGAGCGCGAGAGAATTTGTGGGACTGAGAAGGCGCGTTAGCGTCGCATCAGGCATCGGCGCACGACTGCTGGATGCGGCGGAAACGCCTTATCCGGCCTACGAATGGCGCTGGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAtt                                                                                                <  2:336458/299‑1 (MQ=14)
                                                                                                                                                                                                    tGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                                                    <  1:110816/300‑1 (MQ=255)
                                                                                                                                                                                                               gaaGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAAt                                                                         <  1:108539/299‑1 (MQ=255)
                                                                                                                                                                                                                                gCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTCCGCGGc                                                        >  1:266157/1‑300 (MQ=255)
                                                                                                                                                                                                                                gCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGACCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc                                                        >  1:269624/1‑300 (MQ=255)
                                                                                                                                                                                                                                 cATCAGGCA‑TCTGCGCATGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGAAGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGATGACGCGGGAATTTGTAGGCATGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGAATGCCGGATGCGAGAAACGTCTTGTCCAGCCTACATTTCGACCGTAGGCAGTAATTAAATAGTTCCGACTACGCGGCt                                                       >  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)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                      
TGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCA‑TCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTCGCAACGAGGTTCATCTTATGT  >  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: ‑