Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,550,6630AC25.6% 83.7 / 12.0 43N68T (AAC→ACC) fdnHformate dehydrogenase N subunit beta
Reads supporting (aligned to +/- strand):  ref base A (26/6);  new base C (0/11);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 2.15e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTG  >  NC_000913/1550532‑1550795
                                                                                                                                   |                                                                                                                                    
tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg                                                                                                                              >  2:162700/1‑140 (MQ=255)
         aGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAgcg                                                                                                                     >  2:905956/1‑140 (MQ=255)
           cAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGcc                                                                                                                   <  1:534732/140‑1 (MQ=255)
           cAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGcc                                                                                                                   <  1:566461/140‑1 (MQ=255)
           cAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGcc                                                                                                                   >  1:574285/1‑140 (MQ=255)
           aaGAAGTCGAAAAACTTATCGACTTTCCCCCCTGTATCGGCTGTAAAGCTTGTCGGGGGGCGTGTTGGGGGGGGACCGACATCGGTGAGGAAGGGGGGCACTGCGTGGGGTTTTACGATAACCCCGCCGATCTGAGCGcc                                                                                                                   <  2:574285/139‑1 (MQ=255)
            aGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCa                                                                                                                  >  1:1226391/1‑140 (MQ=255)
                gTCGAAAAACTTATCGACTTTTCCACCTGTAGCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGACCGACATCCGTGATGAGGGGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTc                                                                                                              <  2:1251394/140‑1 (MQ=37)
                   gCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTg                                                                                                           >  2:224739/1‑140 (MQ=255)
                   gCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTg                                                                                                           >  2:156806/1‑140 (MQ=255)
                           tATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGa                                                                                                   >  1:675738/1‑140 (MQ=255)
                            aTCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGAt                                                                                                  >  1:324171/1‑140 (MQ=255)
                             tCGACGTTTCCCGCTGTATCGGCTGTAAAGCCTGTCGGGGGGCGTGTTCGAAGTGGAACCACATCCGTGATGAGGTGGGGCACTGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATg                                                                                                 <  2:1226391/140‑1 (MQ=37)
                             tCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATg                                                                                                 <  1:843978/140‑1 (MQ=255)
                                          cTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCAGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGaaa                                                                                    <  2:1119182/140‑1 (MQ=255)
                                             tATGGGCTGAAAAGCTTGTCGGGGGGGGTTTTGGGGGGGAAACGACATCGGTGATGAGGGGGGGCACGGCGTGGGGTTTAACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCg                                                                                 <  2:1384650/140‑1 (MQ=255)
                                               tCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGACCGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAg                                                                                         <  1:1105612/130‑1 (MQ=255)
                                               tCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAg                                                                                         >  2:1105612/1‑130 (MQ=255)
                                                  aaggaaaaGCCTGTCAGGTGGCGTGTTGGGAGTGGAACGACACCCGTGATGAAGGGGGGCACTGCGTCGGGTTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAg                                                                            <  2:856459/135‑1 (MQ=255)
                                                   cTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGa                                                                           >  1:618418/1‑140 (MQ=255)
                                                             tGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCt                                                                 >  2:1051216/1‑140 (MQ=255)
                                                                cggggggCGTGTTGGGAGGGGAACGACATCCGTGATGAAGGGGGGCCGGGGGTGGGGTTTTACGAAACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa                                                              <  2:566053/135‑1 (MQ=255)
                                                                cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCtctg                                                              >  1:856459/1‑137 (MQ=255)
                                                                cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa                                                              >  1:9484/1‑140 (MQ=255)
                                                                cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa                                                              >  1:165901/1‑140 (MQ=255)
                                                                      gCGTGTTGGGAGTGGAACGACATCGGTGATGAAGGGGGGCACTGCGGGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCt                                                        <  2:165901/140‑1 (MQ=255)
                                                                      gCGTGTTCGGAGTGGAACGACATCCGTGATGAGGTGGGGCCCGGGGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCt                                                        <  1:905956/140‑1 (MQ=255)
                                                                         tgtTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGAt                                                     >  2:78782/1‑140 (MQ=255)
                                                                         tgtTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGAt                                                     >  2:4242/1‑140 (MQ=255)
                                                                           ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGC‑‑TAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGt                                                 >  1:369520/1‑140 (MQ=255)
                                                                                  tGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGa                                            <  1:1051216/140‑1 (MQ=255)
                                                                                  tGGAACCACATCCGGTATAGTGGGGGGCAGTGTGTGGGTGTTCACGAGACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGa                                            <  1:168138/140‑1 (MQ=255)
                                                                                      aCGACATCGGTGATGAGGGGGGGCACTGCGTGGGGTTTTAGGAAAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGc                                        <  2:324171/140‑1 (MQ=255)
                                                                                         aCATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg                                          <  1:1143335/135‑1 (MQ=255)
                                                                                         aCATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg                                          >  2:1143335/1‑135 (MQ=255)
                                                                                            tCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                  <  1:4242/140‑1 (MQ=255)
                                                                                            tCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                  <  1:78782/140‑1 (MQ=255)
                                                                                              ggTGATGAAGGGGGGCCCGGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                <  2:304572/139‑1 (MQ=255)
                                                                                              cGTGATGAGGTGGGGCACTGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                <  1:1117432/140‑1 (MQ=255)
                                                                                              cGTGATGAGGGGGGCCCCTGCGGGGGGGTTTACAATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                <  2:1266004/140‑1 (MQ=255)
                                                                                              cGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                <  1:703723/140‑1 (MQ=255)
                                                                                                 gatgaGGGGGGGCCCTGGGGGGGGTTTCACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGC‑‑TAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtgtg                           <  2:369520/140‑1 (MQ=37)
                                                                                                     aaGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGaa                         >  1:531710/1‑140 (MQ=255)
                                                                                                          gggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATcc                    >  1:649995/1‑140 (MQ=255)
                                                                                                           ggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATccc                   >  1:1382102/1‑140 (MQ=255)
                                                                                                                tGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTg              >  2:882508/1‑140 (MQ=255)
                                                                                                                 gCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGc             >  2:181792/1‑140 (MQ=255)
                                                                                                                  cGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGcc            >  2:414910/1‑140 (MQ=255)
                                                                                                                  cGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGcc            >  1:87477/1‑140 (MQ=255)
                                                                                                                       gggTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAg       <  2:9484/140‑1 (MQ=255)
                                                                                                                          tttACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCg    <  2:618418/140‑1 (MQ=255)
                                                                                                                            tACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTg  >  2:654896/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                    
TGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTG  >  NC_000913/1550532‑1550795

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