Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913372,6400TG29.8% 72.3 / 16.6 47L176V (TTA→GTA) mhpD2‑hydroxypentadienoate hydratase
Reads supporting (aligned to +/- strand):  ref base T (7/26);  new base G (14/0);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 3.40e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCT  >  NC_000913/372503‑372779
                                                                                                                                         |                                                                                                                                           
ttGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTa                                                                                                                                           >  1:917159/1‑140 (MQ=255)
ttGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTa                                                                                                                                           >  2:821553/1‑140 (MQ=255)
 tGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAg                                                                                                                                          <  2:492133/140‑1 (MQ=255)
            ttCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAc                                                                                                                               >  2:796875/1‑140 (MQ=255)
             tCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACt                                                                                                                              >  1:129745/1‑140 (MQ=255)
             tCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACt                                                                                                                              >  1:415298/1‑140 (MQ=255)
             tCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACt                                                                                                                              >  2:224584/1‑140 (MQ=255)
              ccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGTGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTg                                                                                                                             <  1:502125/140‑1 (MQ=255)
              ccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTg                                                                                                                             <  2:937950/140‑1 (MQ=255)
                 gcgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTg                                                                                                                             <  2:1147820/137‑1 (MQ=255)
                 gcgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTg                                                                                                                             >  1:1147820/1‑137 (MQ=255)
                  cgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAACGCGcc                                                                                                                         >  2:1283455/1‑140 (MQ=255)
                   gcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCa                                                                                                                        <  2:415298/140‑1 (MQ=255)
                     tGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAATGCGCCCTg                                                                                                                      >  2:195337/1‑140 (MQ=255)
                           tggtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATg                                                                                                                <  2:798686/140‑1 (MQ=255)
                           tggtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATg                                                                                                                <  1:224584/140‑1 (MQ=255)
                              tggGGAGCCGCGTTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGAcg                                                                                                             <  2:1282275/140‑1 (MQ=255)
                               ggggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGAcgc                                                                                                            >  1:822973/1‑140 (MQ=255)
                                gggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGAcgcg                                                                                                           <  1:1283455/140‑1 (MQ=255)
                                   aGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGtaa                                                                                                        <  2:161142/140‑1 (MQ=255)
                                   aGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGtaa                                                                                                        >  1:254111/1‑140 (MQ=255)
                                   aGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGtaa                                                                                                        >  1:161142/1‑140 (MQ=255)
                                     ccGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGgaata                                                                                                      >  1:980446/1‑140 (MQ=255)
                                          ttCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGaa                                                                                                 <  2:129745/140‑1 (MQ=255)
                                          ttCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGaa                                                                                                 <  2:953071/140‑1 (MQ=255)
                                          ttCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGaa                                                                                                 <  1:821553/140‑1 (MQ=255)
                                          ttCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGaa                                                                                                 <  1:715351/140‑1 (MQ=255)
                                                                   cGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCg                                                                        <  2:917159/140‑1 (MQ=255)
                                                                   cGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCg                                                                        <  2:566588/140‑1 (MQ=255)
                                                                   cGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCg                                                                        >  2:479404/1‑140 (MQ=255)
                                                                    gATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGa                                                                       <  1:796875/140‑1 (MQ=255)
                                                                    gATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGa                                                                       <  2:86215/140‑1 (MQ=255)
                                                                     aTACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGaa                                                                      <  1:625287/140‑1 (MQ=255)
                                                                      tACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCgg                                                                             >  1:506814/1‑132 (MQ=255)
                                                                      tACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCgg                                                                             <  2:506814/132‑1 (MQ=255)
                                                                        ccGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGc                                                                   >  2:707146/1‑140 (MQ=255)
                                                                                        ctcctGTGGGGTGTATGTCATCGGCGGTCCGGGGCAAAGTCCGGCGGGGGTAGACCTGAAAAAATGCGCCATGGAGATGAAGCGTAATAACGGAGAGGGTTTTAGCGGGCGCGGCAGCGAAAACCTGGGGCATCCGCTTa                                                   >  2:242856/1‑140 (MQ=255)
                                                                                            tgtgGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGc                                               >  2:1218062/1‑140 (MQ=255)
                                                                                             gtgGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCg                                              <  1:479404/140‑1 (MQ=255)
                                                                                               ggggTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGc                                            <  1:242856/140‑1 (MQ=255)
                                                                                                   tgtATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTc                                        >  1:403102/1‑140 (MQ=255)
                                                                                                     tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCCGCCGTctg                                      >  1:14525/1‑140 (MQ=255)
                                                                                                          cATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCCTGAAGATGACGCGTAATAACGTAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctggctgg                                 >  2:1359108/1‑140 (MQ=255)
                                                                                                           aTCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCtggctggc                                >  1:148127/1‑140 (MQ=255)
                                                                                                                 ggTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCaa                          >  1:325354/1‑140 (MQ=255)
                                                                                                                  gTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCaaa                         <  1:195337/140‑1 (MQ=255)
                                                                                                                   tCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCaaaa                        >  2:369354/1‑140 (MQ=255)
                                                                                                                   tCCGGCGCAACGTCCCGCGGGGGTAGATCTGAAAAAAAGCCCCCTGAAAAAGACGCGGAAAAACGAAAAGGTTTCTAGCGGGCGCGGCAGGGAAAACCCGGGGCAACCGCGTTAAGCGGCCGTCTCGCTGCCACACacaa                        >  2:656425/1‑140 (MQ=255)
                                                                                                                       gcgcAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc                    <  1:707146/140‑1 (MQ=255)
                                                                                                                       gcgcAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc                    <  2:1291603/140‑1 (MQ=255)
                                                                                                                       gcgcAACGTCCGGCGGGCTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc                    <  2:1205629/140‑1 (MQ=255)
                                                                                                                                         ttAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCt  <  2:27431/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
TTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCT  >  NC_000913/372503‑372779

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