Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,550,6630AC26.1% 92.6 / 10.4 46N68T (AAC→ACC) fdnHformate dehydrogenase N subunit beta
Reads supporting (aligned to +/- strand):  ref base A (30/4);  new base C (0/12);  total (30/16)
Fisher's exact test for biased strand distribution p-value = 4.68e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTC  >  NC_000913/1550528‑1550801
                                                                                                                                       |                                                                                                                                          
tGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAAcccc                                                                                                                                        >  2:345501/1‑140 (MQ=255)
  cGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgc                                                                                                                                      >  1:134104/1‑140 (MQ=255)
    tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg                                                                                                                                    >  1:959748/1‑140 (MQ=255)
    tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg                                                                                                                                    >  2:570434/1‑140 (MQ=255)
            aaGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAgc                                                                                                                            >  1:1246058/1‑140 (MQ=255)
               aaGAAGCCGCAAAACTTATCGACGTTTCCACTTGTTTCGGCTGAAATGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGAAACCGGTGAGGAGGGGGGGCTCTGCGTGGGGTTTTACGATAACCCCGCCGATCTGAGCGcc                                                                                                                         <  2:318697/139‑1 (MQ=21)
                aGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCa                                                                                                                        >  1:1013330/1‑140 (MQ=255)
                aGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCa                                                                                                                        >  2:897782/1‑140 (MQ=255)
                 gAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCaa                                                                                                                       >  1:995560/1‑140 (MQ=255)
                    gTCGCAAATCTTACCGACTTTTCCACCTGTACCGGCTGTAAAGCCTGTCGGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTGGGGGTTTACAAAACCCCCGCCGATCTGAGCGCCAAGTc                                                                                                                    <  2:134104/140‑1 (MQ=255)
                                           cgccTGTATCGGCTGAAAAGCCTGTCGGGTGGCGTGTTGGGGGTGGACCGACATCGTTGAGGAGGGGGGGCACTGCGGCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCg                                                                                             <  2:1311022/138‑1 (MQ=255)
                                                 tATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCg                                                                                       >  2:717850/1‑140 (MQ=255)
                                                 tATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCg                                                                                       >  1:273705/1‑140 (MQ=255)
                                                 tATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCg                                                                                       >  2:1290313/1‑140 (MQ=255)
                                                   tCGGCTGTAAAGCCTGTCAGGTGGGGTGTTCGGAGTGGAACGACATCCGTGATGAAGGGGGGCACGGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGaa                                                                                     <  2:1380442/140‑1 (MQ=38)
                                                        tGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGaa                                                                                <  1:482797/140‑1 (MQ=255)
                                                           aaaGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACgg                                                                             >  2:521972/1‑140 (MQ=255)
                                                           aaaGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACgg                                                                             <  1:1229620/140‑1 (MQ=255)
                                                                 tGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCt                                                                       >  2:425337/1‑140 (MQ=255)
                                                                 tGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGGTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCt                                                                       >  2:1059633/1‑140 (MQ=255)
                                                                   tCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTgg                                                                     >  2:263359/1‑140 (MQ=255)
                                                                        tGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTgg                                                                <  1:897782/140‑1 (MQ=255)
                                                                          gCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCt                                                              <  2:521192/140‑1 (MQ=255)
                                                                           ggtgtTGGGAGTGGAACGACATCCGTGATGAGGGGGGGCCCTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTg                                                             <  1:345501/139‑1 (MQ=255)
                                                                               ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                         >  1:480299/1‑140 (MQ=255)
                                                                               ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                         <  1:263359/140‑1 (MQ=255)
                                                                               ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                         >  1:826303/1‑140 (MQ=255)
                                                                                 cGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGt                                                       >  1:391875/1‑140 (MQ=255)
                                                                                      tGGAACGACATCCGTGATGAAGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGa                                                  <  1:1480149/140‑1 (MQ=255)
                                                                                      gggACCGACATCCGTGAGGAGGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGa                                                  <  2:391875/139‑1 (MQ=255)
                                                                                      gggAACGACATCCGTGATGAAGGGGGGCACTGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGa                                                  <  2:826303/139‑1 (MQ=255)
                                                                                      gggAACGACACCCGTGAGGAGGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGa                                                  <  2:1013330/139‑1 (MQ=255)
                                                                                          aCGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGc                                              <  1:521972/140‑1 (MQ=255)
                                                                                             aCATCCGTGATGAAGTGGGGCACTGCGTCGGGGGTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGt                                           <  1:1059633/140‑1 (MQ=255)
                                                                                                tCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                        >  2:638007/1‑140 (MQ=255)
                                                                                                tCCGTGATGAAGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                        <  2:480299/140‑1 (MQ=255)
                                                                                                cccGTGAGGAAGGGGGGCACGGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                        <  2:463212/139‑1 (MQ=255)
                                                                                                   aagatgaGGGGGGCCCCTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCAc                                     <  2:785460/138‑1 (MQ=37)
                                                                                                              gggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATcc                          >  2:1367718/1‑140 (MQ=255)
                                                                                                              gggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATcc                          >  1:990933/1‑140 (MQ=255)
                                                                                                              gggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATcc                          >  1:548448/1‑140 (MQ=255)
                                                                                                               ggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                      <  1:998714/127‑1 (MQ=255)
                                                                                                               ggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                      >  2:998714/1‑127 (MQ=255)
                                                                                                                 cACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                       <  2:959748/140‑1 (MQ=255)
                                                                                                                 cACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATccctg                       >  1:785460/1‑138 (MQ=255)
                                                                                                                 cACGGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                       <  2:911575/140‑1 (MQ=255)
                                                                                                                   cTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCt                     >  1:523315/1‑140 (MQ=255)
                                                                                                                     gCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGc                   >  1:3674/1‑140 (MQ=255)
                                                                                                                      cGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGcc                  >  1:12443/1‑140 (MQ=255)
                                                                                                                        tCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTg                >  1:1068005/1‑140 (MQ=255)
                                                                                                                              tttACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCg          >  2:401391/1‑140 (MQ=255)
                                                                                                                                tACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTg        >  1:1070818/1‑140 (MQ=255)
                                                                                                                                      aaCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTc  >  2:548219/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                          
TGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTC  >  NC_000913/1550528‑1550801

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