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

CATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAA  >  NC_000913/372501‑372774
                                                                                                                                           |                                                                                                                                      
cATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCGTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGt                                                                                                                                        >  1:639811/1‑140 (MQ=255)
cATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCGTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGt                                                                                                                                        <  2:639811/140‑1 (MQ=255)
cATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGt                                                                                                                                        <  1:685913/140‑1 (MQ=255)
 aTTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGtt                                                                                                                                       <  1:366609/140‑1 (MQ=255)
  ttGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTa                                                                                                                                      >  1:63034/1‑140 (MQ=255)
          gTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGaa                                                                                                                              >  2:610026/1‑140 (MQ=255)
             cTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGaaaaa                                                                                                                           <  2:925580/140‑1 (MQ=255)
              ttCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAAc                                                                                                                          >  2:471012/1‑140 (MQ=255)
               tCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACt                                                                                                                         <  1:610026/140‑1 (MQ=255)
                ccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTg                                                                                                                        <  2:917491/140‑1 (MQ=255)
                ccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTg                                                                                                                        <  2:108453/140‑1 (MQ=255)
                ccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTg                                                                                                                        <  2:777880/140‑1 (MQ=255)
                   gcgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTa                                                                                                                                      <  1:530627/123‑1 (MQ=255)
                   gcgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTa                                                                                                                                      >  2:530627/1‑123 (MQ=255)
                     gcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCa                                                                                                                   <  2:939134/140‑1 (MQ=255)
                           acaggtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGa                                                                                                             <  1:118487/137‑1 (MQ=255)
                              ggtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGa                                                                                                          <  1:471012/140‑1 (MQ=255)
                                  gggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGAcgcg                                                                                                      <  2:90609/140‑1 (MQ=255)
                                  gggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGAcgcg                                                                                                      <  1:677222/140‑1 (MQ=255)
                                   agacaggCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGt                                                                                                     <  1:24240/134‑1 (MQ=255)
                                          cATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGt                                                                                                     <  1:587974/133‑1 (MQ=255)
                                          cATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGt                                                                                                     >  2:587974/1‑133 (MQ=255)
                                          cATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGGCCCGCGGGGTTAGAACCGAAAAAATGCGCCATGAAAATGACGCGtaaactgt                                                                                              >  2:118487/1‑135 (MQ=255)
                                            ttCGCGACTGCTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGaa                                                                                            <  2:882701/140‑1 (MQ=255)
                                                 gACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGt                                                                                       >  2:296544/1‑140 (MQ=255)
                                                    tGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTTGACCTGAAAAAATGCGCCATGAAGAAGACGCGTAATAAagaagactgtct                                                                                    >  2:24240/1‑134 (MQ=255)
                                                      gTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTa                                                                                  >  1:801893/1‑140 (MQ=255)
                                                      gTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGGCCTTGAAAACTGCGCCATGAAGATGACGCGGAATAACGAAGAGGGTTTTa                                                                                  >  2:833371/1‑140 (MQ=255)
                                                      gTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGGTTCtt                                                                                  >  2:890410/1‑139 (MQ=255)
                                                        cGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGc                                                                                <  1:307602/140‑1 (MQ=255)
                                                                     cGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCg                                                                   <  2:801893/140‑1 (MQ=255)
                                                                         aCCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAATGCGCCATGAAGATGACCGGGAAAAAAGAAGAGGGTTTTAGCGGGGGCGGCAGCCaaat                                                               >  2:940351/1‑138 (MQ=255)
                                                                              ggCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTgg                                                          <  1:833371/140‑1 (MQ=255)
                                                                                  gATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACa                                                      >  1:797925/1‑140 (MQ=255)
                                                                                         cctcctGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCtt                                               >  2:854327/1‑140 (MQ=255)
                                                                                          ctcctGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAATGCGCCCTGAAAATGACGCGGAAAAAAGAAGAGGGTTTTAGCGGGCGCGGCAGCGAAAGCCTGGGGCAACCGCCTa                                              >  2:766754/1‑140 (MQ=255)
                                                                                             gagacaggTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATg                                           <  1:484550/134‑1 (MQ=255)
                                                                                              tgtgGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGc                                          >  1:392416/1‑140 (MQ=255)
                                                                                                     tgtATGTCATCGGCGGTCCGGCGCAACGTCCCGCGGGGGTAGACCTGTAAAACCCCGCCCTGAAGAAGACGCGTAATAACGAAGAGGTTTCTATCCGGCGCCGCAGCGAATGCCCTGGACATCCGCTTTATTCGGCCgcc                                   >  2:548807/1‑138 (MQ=255)
                                                                                                       tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctg                                 >  2:599064/1‑140 (MQ=255)
                                                                                                       tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctg                                 >  1:811583/1‑140 (MQ=255)
                                                                                                       tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGACAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctg                                 >  2:135428/1‑140 (MQ=255)
                                                                                                       tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGGAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCCTAATGCGGCCGTctg                                 >  2:211689/1‑140 (MQ=255)
                                                                                                       tATGTCATCGGCGGTCCGGCGCAACGGCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAAAGAAGAGGGTTCTAGCGGGCGCGGCAGCGAAAGCCTGGGACATCCGCTTAATGCGGCCGTctt                                 >  2:599550/1‑139 (MQ=255)
                                                                                                       tATGTCATCGGCGGTCCGGCGAAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctg                                 >  2:699932/1‑140 (MQ=255)
                                                                                                        aTGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctgg                                <  2:811583/140‑1 (MQ=255)
                                                                                                        aTGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctgg                                >  1:546488/1‑140 (MQ=255)
                                                                                                            cATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctggctgg                            >  1:373327/1‑140 (MQ=255)
                                                                                                            cATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctggctgg                            >  1:289537/1‑140 (MQ=255)
                                                                                                             aTCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCtggctggc                           >  1:711487/1‑140 (MQ=255)
                                                                                                              tCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCCTGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCa                          >  2:695468/1‑140 (MQ=255)
                                                                                                               cggcggTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCCTGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGGGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTTGCTGGCAc                         >  2:981443/1‑140 (MQ=255)
                                                                                                               cggcggTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGGAATAAAGAAGAGGGTTCTAGCGGGGGCGGGAGCGAAAGCCTGGGGCATCCGCCTTATGCGGCCGTctggctgtctc                         >  2:484550/1‑136 (MQ=255)
                                                                                                               cggcggTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAATGCGCCCTGAAGATGACGCGTAATAACGAAAAGGGTTTTTGCGGGCGCGGCAGCGAAAGCCCGGGGCATCCCCTTAATTCGGCCGGCTGGCTGGCAc                         >  2:495205/1‑140 (MQ=255)
                                                                                                                     tCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCaaaa                   >  1:301156/1‑140 (MQ=255)
                                                                                                                       cGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATg                 <  1:548807/140‑1 (MQ=255)
                                                                                                                         gcgcAACGTCCGGCGGGGTTAGACCTGAAAGACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc               <  2:711487/140‑1 (MQ=255)
                                                                                                                         gcgcAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc               <  1:495205/140‑1 (MQ=255)
                                                                                                                                      cggGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGaa  >  2:834753/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                      
CATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAA  >  NC_000913/372501‑372774

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: