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

CGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTC  >  NC_000913/1550530‑1550801
                                                                                                                                      |                                                                                                                                          
cGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgc                                                                                                                                      >  2:694466/1‑140 (MQ=255)
  tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg                                                                                                                                    >  1:974344/1‑140 (MQ=255)
    aTTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGAt                                                                                                                                  <  1:213822/140‑1 (MQ=255)
       aCAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTg                                                                                                                               >  1:216583/1‑140 (MQ=255)
           aGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAgcg                                                                                                                           >  1:820629/1‑140 (MQ=255)
             aaGAAGTCCCAAAGTTTATCGCGTTTCCCACGTGGTTGGGCTGAAAAGCTTGGCGGGTGGCGTGTTGGGAGGGAAACGACATCCGGGATGAGGGGGGGACCGGCGGCGGGGTTTAGGATACCCCCGCCGATCTGAGCGcc                                                                                                                         <  2:93550/139‑1 (MQ=255)
                            ttATCGCCTTTTCCCCCTGTATCGGCTGTAAAGCCTGTAAGGGGGCGTGTGGGGGGGGGAACGACATCCGGGATGAAGTGGGGCACTGCGGGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTg                                                                                                          <  2:150231/140‑1 (MQ=255)
                            ttATCGACTTTTCCACCTGTATCGGCTGTAAAGCCTGTCGGGTGGCGTGTTGGGAGTGGAACGACATCGGTGATGAGGGGGGGCACTGCGTGGGGTTTTACAATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTg                                                                                                          <  2:943939/140‑1 (MQ=25)
                              aTCGACGTTTCCCCCTGTATCGGCTGTAAAGCCTGTCGGGGGGCGTGTTCGGAGTGGAACGACATCGGTGATGAGGTGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGAt                                                                                                        <  2:452150/140‑1 (MQ=37)
                                     tttCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCttt                                                                                                 >  1:1224216/1‑140 (MQ=255)
                                      ttCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGTGCTTTa                                                                                                >  1:335428/1‑140 (MQ=255)
                                                tATGGGCTGTAAAGCCTGTCGGGTGGCGTGTTGGGAGTGGAACGACATCCGTGATGAGGGGGGGCACTGGGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGTGCTTTAGCGAAACCg                                                                                       <  2:335428/140‑1 (MQ=255)
                                                 aTCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGa                                                                                      >  1:372940/1‑140 (MQ=255)
                                                 aTCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGa                                                                                      >  1:707203/1‑140 (MQ=255)
                                                  tCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGaa                                                                                     >  2:57803/1‑140 (MQ=255)
                                                                tGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCt                                                                       >  2:1208093/1‑140 (MQ=255)
                                                                tGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCt                                                                       >  1:993815/1‑140 (MQ=255)
                                                                   cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa                                                                    >  1:960017/1‑140 (MQ=255)
                                                                   cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa                                                                    >  1:1184611/1‑140 (MQ=255)
                                                                   cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa                                                                    >  2:881767/1‑140 (MQ=255)
                                                                        ggCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGc                                                               >  1:16865/1‑140 (MQ=255)
                                                                         gCGTGTTGGGAGTGGAACGACATCCGTGATGAAGTGGGGCCCTGGGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCt                                                              <  2:1184611/140‑1 (MQ=255)
                                                                              ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                         >  1:975723/1‑140 (MQ=255)
                                                                              ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                         >  1:1388628/1‑140 (MQ=255)
                                                                                acaggTGGAACGACATCGGTGATGAGGGGGGGCCCGGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGt                                                       <  2:1028218/136‑1 (MQ=255)
                                                                                   aGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTaaa                                                    >  2:770044/1‑140 (MQ=255)
                                                                                     gggAACGACACCCGTGATGAAGGGGGCCACTGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGa                                                  <  2:372940/139‑1 (MQ=255)
                                                                                      ggAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGAc                                                 >  1:692124/1‑140 (MQ=255)
                                                                                       gAACGACATCCGTGATGAGGTGGGGCACTGGGGCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg                                                <  2:974344/140‑1 (MQ=255)
                                                                                       gAACGACATCCGTGATGAAGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg                                                <  2:820629/140‑1 (MQ=255)
                                                                                             cATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTa                                          >  1:960051/1‑140 (MQ=255)
                                                                                             cATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTa                                          >  2:1259423/1‑140 (MQ=255)
                                                                                               tCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                        <  1:770044/140‑1 (MQ=255)
                                                                                               tCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTgctgt                                        >  1:1028218/1‑136 (MQ=255)
                                                                                               ccggTGATGAAGTGGGCCACGGCGTGGGGTTTTACAATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                        <  2:975723/137‑1 (MQ=255)
                                                                                                 cGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg                                                <  1:928030/130‑1 (MQ=255)
                                                                                                 cGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg                                                >  2:928030/1‑130 (MQ=255)
                                                                                                  gTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCAc                                     >  1:402666/1‑140 (MQ=255)
                                                                                                  gTGATGAAGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCAc                                     <  2:16865/140‑1 (MQ=255)
                                                                                                    gatgaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtg                                   >  2:600643/1‑140 (MQ=255)
                                                                                                                cACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                       <  1:57803/140‑1 (MQ=255)
                                                                                                                cACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                       <  1:881767/140‑1 (MQ=255)
                                                                                                                cACTGCGGGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                       <  2:216583/140‑1 (MQ=255)
                                                                                                                  cTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCt                     >  2:138217/1‑140 (MQ=255)
                                                                                                                      gTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCt                 <  1:694466/140‑1 (MQ=255)
                                                                                                                            gTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGc           >  2:1145334/1‑140 (MQ=255)
                                                                                                                             tttACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCg          >  2:1377694/1‑140 (MQ=255)
                                                                                                                             tttACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCg          >  1:204799/1‑140 (MQ=255)
                                                                                                                               tACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTg        >  1:1275573/1‑140 (MQ=255)
                                                                                                                                   aTAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCg    <  1:600643/140‑1 (MQ=255)
                                                                                                                                     aaCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTc  >  2:1353753/1‑140 (MQ=255)
                                                                                                                                     aaCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTc  >  2:1003027/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                          
CGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTC  >  NC_000913/1550530‑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: