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

GGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTC  >  NC_000913/1550526‑1550801
                                                                                                                                         |                                                                                                                                          
ggTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAAcc                                                                                                                                          >  2:462326/1‑140 (MQ=255)
   gCGTGATTACAAAGCAGAAGTCGAAAAACTTACCGTCTTTTCCCCCTGTATCGGCTGTAAAGCCTGTCGGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAGGTGGGGCACTGCGTCGGGTTTTACGATACCCccg                                                                                                                                       <  2:8971/140‑1 (MQ=37)
    cGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgc                                                                                                                                      >  2:733266/1‑140 (MQ=255)
      tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg                                                                                                                                    >  2:65919/1‑140 (MQ=255)
      tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg                                                                                                                                    >  2:842206/1‑140 (MQ=255)
       gATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGa                                                                                                                                   >  2:664388/1‑140 (MQ=255)
       gATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGa                                                                                                                                   >  2:142784/1‑140 (MQ=255)
                 aaGAAGTCGAAAAACTTATCGACGTTTCCACCTGTATCGGCTGAAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGACCGACATCCGTGATGAAGGGGGGCACTGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGcc                                                                                                                         <  2:102780/139‑1 (MQ=25)
                   aaaGTGGAAAAACTTATCGACTTTCCCCCCTGTGTCGGCTGTAAGGCCTGTCGGGTGGCTTGTTGGGGGTGCAACGCCACCCGTGATGAGGGGGGGCCCGGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCaa                                                                                                                       <  2:868134/139‑1 (MQ=255)
                      gccGCAAAACTTACCGACGTTCCCACCTGTATCGGTTGTAAAGCCTGTCAGGTGGGGTGTTCGGAGTGGAACGACATCCGTGAGGAAGGGGGGCACGGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTc                                                                                                                    <  2:296230/138‑1 (MQ=255)
                      gTCGAAAAATTTATCGACGTTCCCCCCTGTAGCGGTTGAAAGGCTTGACAGGGGGCGTGTTCGGGGTGGACCGACATCCGTGATGAAGGGGGGCACTGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTc                                                                                                                    <  2:520682/140‑1 (MQ=255)
                           aaaaCTTATCGACGTTTCCCCCTGTATCGGCTGTAAAGCCTGTCGGGTGGCGTGTTGGGAGTGGAACGACATCGGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGa                                                                                                               <  2:518050/140‑1 (MQ=37)
                           aaaaCTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGa                                                                                                               <  2:311608/140‑1 (MQ=255)
                                      ccGTTCCCCCCTGTATGGGTTGTAAAGCTTGTCAGGTGGCGTGTTCGGAGTGGACCGACATCCGTGATGAAGTGGGGCACTGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATgcgc                                                                                                    <  2:882900/139‑1 (MQ=255)
                                         tttCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCttt                                                                                                 >  1:882900/1‑140 (MQ=255)
                                                     tCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGaa                                                                                     >  1:491073/1‑140 (MQ=255)
                                                                  cTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGc                                                                        <  1:57160/140‑1 (MQ=255)
                                                                   tGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCt                                                                       >  2:573823/1‑140 (MQ=255)
                                                                     tgaGGGGGCGTGTTGGGAGTGGAACGACATCCGTGATGAAGGGGGGCCCTGCGGGGGGTTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTgg                                                                     <  2:934191/138‑1 (MQ=255)
                                                                     tcgggTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGGGGGGCACTGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAACCCGAACAGAACGGCAAGCTgg                                                                     <  2:348527/137‑1 (MQ=255)
                                                                      cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa                                                                    >  2:371202/1‑140 (MQ=255)
                                                                                      aGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTaaa                                                    >  1:232500/1‑140 (MQ=255)
                                                                                        gggAACAATACCCGGGATGAGGTGGGCCCCGGCGTCGGGTTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGa                                                  <  2:386212/139‑1 (MQ=255)
                                                                                         ggAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGAc                                                 >  2:664521/1‑140 (MQ=255)
                                                                                          gAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCATAAAGACg                                                <  1:65919/140‑1 (MQ=255)
                                                                                                 aTCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTAt                                         <  1:371202/140‑1 (MQ=255)
                                                                                                  ccgttaatgaGGTGGGCCACTGGGTGGGGTTTTGCGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                        <  2:232500/134‑1 (MQ=255)
                                                                                                    cGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                      <  1:664388/140‑1 (MQ=255)
                                                                                                    cGTGATGAAGGGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                      <  1:352588/140‑1 (MQ=255)
                                                                                                     gTGATGAAGTGGGGCCCTGCGTGGGGTTTTAGAATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCAc                                     <  2:48912/140‑1 (MQ=38)
                                                                                                     gTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCAc                                     <  1:733266/140‑1 (MQ=255)
                                                                                                     gTGATGAAGTGGGGCACGGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCAc                                     <  2:723448/140‑1 (MQ=255)
                                                                                                        atgaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtgt                                  >  2:447231/1‑140 (MQ=255)
                                                                                                          gaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGa                                >  2:567836/1‑140 (MQ=255)
                                                                                                               ggggCACTGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATc                           <  1:567836/140‑1 (MQ=255)
                                                                                                                    ccTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGc                      <  1:462326/139‑1 (MQ=255)
                                                                                                                     cTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCt                     >  1:657507/1‑140 (MQ=255)
                                                                                                                     cTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCt                     >  2:416422/1‑140 (MQ=255)
                                                                                                                            ggggTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGaa              >  1:722025/1‑140 (MQ=255)
                                                                                                                                        aaCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTc  >  1:440044/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                          
GGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTC  >  NC_000913/1550526‑1550801

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