Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913370,8880AC48.6% 41.2 / 10.4 36L204F (TTA→TTC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (15/3);  new base C (0/17);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 2.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGAT  >  NC_000913/370761‑371038
                                                                                                                                   |                                                                                                                                                      
ccacccctctCAAAAACGGCGTGTTGTTCTTGGTTTCGGGTGGGTTTTCCCACCACCCGCCGGTC‑CCCGACCGGGCAAAG‑G‑CCGATCCCC‑ATGGGCGCGCCCTTTTGTTGGGGGGGGGAAAAATTTTCCCCGCCAGTGAGCGCGAATTGc                                                                                                                                  <  1:125840/144‑1 (MQ=255)
            aTAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGt                                                                                                                     >  2:704641/1‑151 (MQ=255)
                aCGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAtt                                                                                                                 >  1:540545/1‑151 (MQ=255)
                   cgTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAgc                                                                                                              >  2:400463/1‑151 (MQ=255)
                          ttctggggTCCGGGTGGTTTTCCCCACCGCCCGCGGGGC‑CCCGACCGGGAAAAA‑G‑CGGAGGCCA‑ATAGGCGCGCCCGTGTGTGGGGGGGGGGAAAATATTGCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCt                                                                                                        <  2:250307/145‑1 (MQ=255)
                              tGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTgaga                                                                                                    >  1:799415/1‑150 (MQ=255)
                                ggttggggggggCTTCCCCATCGCCCGCGGGGG‑CCGGAACGGGGGAAA‑G‑CCGTCCCCC‑ATATGCGCGCCCGTTTGTGGGGGAGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAg                                                                                                  <  1:35482/141‑1 (MQ=255)
                                 gTTCCGGAGGGCTTTCCCATCCGCGGGCGGGC‑GCCAAAGGGGGGAAACG‑TCGCTCCCT‑AAATGCGCACCCGTTTGTGGGGGGGCGGAAATGTTTGCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTg                                                                                                       <  1:501338/145‑1 (MQ=255)
                                    ccGGTGGGCTTTCCCATCAGCCGCCGGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgt                                                                                             >  2:543191/1‑151 (MQ=255)
                                      ggggggCTTCCCCACCCGCCGCGGGTG‑CCGGAACGGGGAAAA‑G‑CCGATGCCC‑TTATGCGCGCCCGTTTGTGGGGGGGGGGGAAAGATTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttg                                                                                              <  1:6430/145‑1 (MQ=255)
                                      ggTGGGCTTCCCCACCAGCCGCCGGGC‑CCGGACCTGGAAAAA‑G‑CCGATGCCC‑TTATGCGCGCCCGTCTGTGGGGGGGCGGAAAATTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGtt                                                                                                <  2:11423/146‑1 (MQ=255)
                                           gctccccctccggccgctGGTC‑CCCGATCGGGTGAAA‑G‑CGGACCCCA‑ATGTGCGCGCCCGTCTGTGGGGGGGGGGAAAATATTGTCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAg                                                                                         <  2:349826/136‑1 (MQ=255)
                                           gCTTTCCCATCAGCCGCCGGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGCCAGCAACGGGTGATTAGCGCCGCTGAGAAGCTTGTTGAGGAt                                                                                      >  1:72116/1‑151 (MQ=255)
                                            cTTTCCCATCAGCCGCCGGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATc                                                                                     >  2:257413/1‑151 (MQ=255)
                                                cctccccgccgcgGGTC‑CCCGAGC‑GGGGAAA‑GCCCGAGGCCA‑ATACGCGCGCCTGTTGGGGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGAt                                                                                      <  2:833589/139‑1 (MQ=255)
                                                  cagcggcggcgggTG‑CCCG‑ACGGGGAAAA‑G‑CCCGTGCCCAATATGCGCGGCGTTGTGTTGGGGGGGGGGAAATTTTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGa                                                                                       <  1:744650/132‑1 (MQ=255)
                                                      cggcggcgGTGCCCCGACCGGGCAAAA‑G‑CTGACGCCC‑ATAGGCGCGCCCGTTTGTGGGGGGGGGGAAAAGATTACCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGaa                                                                               <  2:615794/142‑1 (MQ=255)
                                                            cgGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATcc                                                                     >  1:433401/1‑151 (MQ=255)
                                                             ggTC‑CCCGACGTGGAGAAC‑G‑CCGATGCCC‑ATATGCGCGCCCGTCTGTGGGGGGGGGGAAAATATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCt                                                                           <  2:660925/144‑1 (MQ=255)
                                                                   ccGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCa                                                                >  2:325503/1‑150 (MQ=255)
                                                                      aaCTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAccc                                                            >  1:414501/1‑151 (MQ=255)
                                                                          ggCGAAG‑G‑CCGCTCCCC‑ATAGGCGCGACGGTTTGGGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGa                                                          <  1:715228/149‑1 (MQ=255)
                                                                                          ccc‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGtt                                          >  2:156968/1‑151 (MQ=255)
                                                                                                    gcgcccttgtgggggggggggGAAAAGATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCAt                                       <  1:714118/132‑1 (MQ=255)
                                                                                                    gagctcGTCTGTGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACt                                   <  1:551041/144‑1 (MQ=255)
                                                                                                        ccttttgtgGGGGGGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATg                                      <  2:560540/137‑1 (MQ=255)
                                                                                                          tgttgggggggggggggAAAGATTACCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTg                             <  1:14664/142‑1 (MQ=255)
                                                                                                          gttggtgggggggggggAAATTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTgat                              <  1:321955/141‑3 (MQ=255)
                                                                                                          gTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAGTTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTgg                            >  2:173530/1‑150 (MQ=255)
                                                                                                           tgtgtgGGGGGGGGGAAAAGTTTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACttt                                 <  1:277976/142‑1 (MQ=255)
                                                                                                                 ggggggggGAAAAGTTTCACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCt                              <  2:460995/141‑1 (MQ=255)
                                                                                                                   gggggggaaaaGATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACg                  <  1:345873/143‑1 (MQ=255)
                                                                                                                   ggAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACg                  >  2:162402/1‑151 (MQ=255)
                                                                                                                    gggggggAAAGTTTGCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACg                  <  1:369911/146‑1 (MQ=255)
                                                                                                                        gggAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATAc             >  2:517568/1‑151 (MQ=255)
                                                                                                                              gATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAAc       >  1:367602/1‑151 (MQ=255)
                                                                                                                                   aCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGAt  >  2:312587/1‑151 (MQ=255)
                                                                                                                                   aCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGAt  >  1:120123/1‑151 (MQ=255)
                                                                                                                                   |                                                                                                                                                      
CCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTG‑CCCGAACTGGCGAAA‑G‑CCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGAT  >  NC_000913/370761‑371038

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