Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913372,6400TG30.0% 71.7 / 11.1 40L176V (TTA→GTA) mhpD2‑hydroxypentadienoate hydratase
Reads supporting (aligned to +/- strand):  ref base T (6/22);  new base G (12/0);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 3.32e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAA  >  NC_000913/372516‑372774
                                                                                                                            |                                                                                                                                      
tCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACt                                                                                                                         >  2:382051/1‑140 (MQ=255)
 ccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTg                                                                                                                        <  1:618193/140‑1 (MQ=255)
  cGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGTCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTgc                                                                                                                       >  2:1084393/1‑140 (MQ=255)
    gcgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTgcgc                                                                                                                     <  1:116403/140‑1 (MQ=255)
              tggtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATg                                                                                                           <  1:806158/140‑1 (MQ=255)
              tggtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATg                                                                                                           <  1:605570/140‑1 (MQ=255)
                gtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGAc                                                                                                         >  2:48514/1‑140 (MQ=255)
                gtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGAc                                                                                                         >  1:771830/1‑140 (MQ=255)
                 tggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGAcg                                                                                                        >  1:1093345/1‑140 (MQ=255)
                    ggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGt                                                                                                     >  2:379379/1‑140 (MQ=255)
                     gAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCCACGCTCCGCGGGGGTTGACCCGAAAAAAACCGCCCTGTAGATGACGCGtg                                                                                                    >  2:605570/1‑139 (MQ=255)
                     gAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGta                                                                                                    <  2:1093345/140‑1 (MQ=255)
                        ccGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGtaata                                                                                                 >  1:967811/1‑140 (MQ=255)
                             ttCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGaa                                                                                            <  1:379379/140‑1 (MQ=255)
                                 cgACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAgg                                                                                        <  1:48514/140‑1 (MQ=255)
                                  gACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGt                                                                                       >  2:619794/1‑140 (MQ=255)
                                                 tttGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCgg                                                                        >  2:1109104/1‑140 (MQ=255)
                                                 tttGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACCGCGCCATGAAGAAGACGCGTAATAAAGAAGAGGTTTCTAGGGGGCGCgg                                                                        >  2:21457/1‑140 (MQ=255)
                                                      cGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCg                                                                   <  2:380162/140‑1 (MQ=255)
                                                      cGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCg                                                                   <  1:1109104/140‑1 (MQ=255)
                                                      cGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCCTGAAGATGACGCGTAATAACGAAGAGGGTTCTAACGGGCGCGGCAGCg                                                                   >  2:728559/1‑140 (MQ=255)
                                                       gATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGa                                                                  <  1:1018617/140‑1 (MQ=255)
                                                       gATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGa                                                                  <  1:579674/140‑1 (MQ=255)
                                                           ccGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCCGAAAAACTGGGCCCTGAGGATCACGCGGAATAAAAAAGAGGGTTCTAACGGGCGCGGCAGCGAATgg                                                              >  2:1018617/1‑139 (MQ=255)
                                                               ggCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTgg                                                          >  1:758734/1‑140 (MQ=255)
                                                               ggCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAATGCGCCCTGAAGATGACGCGGAATAACGAAGAGGGTTTTAGCGGGCGCGGGAGCGAAAGCCTgg                                                          >  2:579674/1‑140 (MQ=255)
                                                                gCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCg                                                                   >  1:811469/1‑130 (MQ=255)
                                                                gCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCg                                                                   <  2:811469/130‑1 (MQ=255)
                                                                   gATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACa                                                      <  1:382051/140‑1 (MQ=255)
                                                                   gATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACa                                                      <  1:21457/140‑1 (MQ=255)
                                                                                gtgGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCg                                         <  1:912032/140‑1 (MQ=255)
                                                                                gtgGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCg                                         <  2:758734/140‑1 (MQ=255)
                                                                                        tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGGAGCGAAAGCCTGGGGCATCCGCTTAATGCGGGCGGctg                                 >  2:912032/1‑140 (MQ=255)
                                                                                        tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctg                                 >  1:851949/1‑140 (MQ=255)
                                                                                           gTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctggct                              <  1:579500/140‑1 (MQ=255)
                                                                                           gTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctggct                              >  1:97495/1‑140 (MQ=255)
                                                                                              aTCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCtggctggc                           >  2:651121/1‑140 (MQ=255)
                                                                                               tCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCAGTCTGGCTGGCa                          <  1:1084393/140‑1 (MQ=255)
                                                                                                cggcggTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCAc                         >  1:132533/1‑140 (MQ=255)
                                                                                                        cGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATg                 <  2:132533/140‑1 (MQ=255)
                                                                                                              aaCGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGt           <  2:790346/140‑1 (MQ=255)
                                                                                                                cGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCt         >  2:626554/1‑140 (MQ=255)
                                                                                                                  tCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTgg       >  2:865450/1‑140 (MQ=255)
                                                                                                                     ggcggGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTg    <  2:109739/140‑1 (MQ=255)
                                                                                                                       cggGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGaa  <  2:204678/140‑1 (MQ=255)
                                                                                                                            |                                                                                                                                      
TCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAA  >  NC_000913/372516‑372774

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