Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,8440CG20.4% 85.7 / 14.6 56A642A (GCC→GCGrsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base C (27/16);  new base G (5/6);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATT  >  NC_000913/1708719‑1708991
                                                                                                                             |                                                                                                                                                   
cgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAAACGCAATCGAATCCTCTATTGCAAGAGACAAAGAAACCAAGCTCGCCCAGAAAAATACTAATGCAGAGCCAGAAGAACAGGTGGATcc                                                                                                                              >  2:1493609/1‑149 (MQ=2)
                    ggCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCCGTAGCAGCGGCTATTGCCCGAGAGCGGGAGAGCGACCTGGAGCGGAAACATGGTAATGACGGTCCAGAACATCAGGTTGATCAGCGCAAAGCGGCCGtggccg                                                                                                          >  2:1493708/1‑144 (MQ=2)
                         aTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg                                                                                                                                                  >  1:1493504/1‑104 (MQ=18)
                         aTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg                                                                                                                                                  <  2:1493504/104‑1 (MQ=18)
                                  cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTTCTGAAGAGCCACGAGAACGGGTAGATCCTAGCAATGCGGCCGTCGCGGCTGCTAtagaa                                                                                               >  2:1493630/1‑142 (MQ=11)
                                            ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCACCAGGCAAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGGCCAAGCCGCTATAGCCCGAGCCAAAgcgc                                                                                  >  2:3723309/1‑149 (MQ=14)
                                                gATCCTCGCAAAGCTGCCGTCGAAGCAGCTATAGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                                                                                                                                <  1:1493517/99‑1 (MQ=14)
                                                gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaa                                                                              <  2:1493522/149‑1 (MQ=255)
                                                gATCCCCGCAAACCAGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCa                                                                                                                                               <  1:1493518/84‑1 (MQ=14)
                                                 acccccgcAACCCCGCCGCCGAAGCCGCTATTGATTGTTCTAAAGCTCGCAAGCCGTAACAGAAACAGTCTATTTCCGAACAAGAAGATCTTGCCTATTCTCGCAAAGCTGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCaaa                                                                             <  2:1493577/144‑1 (MQ=9)
                                                           cccccccgTCGACGCCGCTACCGCCCCCGCCAAAGCTCTCAAGCTTTAACAGCATCAGCCTACTGCGTACACATCTCAACCGGTCGACGCGCGCAAAGCGTCCTTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc                                                                   <  2:1493589/144‑1 (MQ=255)
                                                              cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGCACAGCacc                                                                >  1:1493659/1‑147 (MQ=255)
                                                               gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                                                                                                                              <  1:1493530/86‑1 (MQ=255)
                                                               gccgTCGAAGCAGCTATTGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                                                                                                                                 <  1:1493525/83‑1 (MQ=11)
                                                                   tCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAgc                                                                                                                     >  1:1493646/1‑91 (MQ=255)
                                                                      tattcGCTATTGCCCGTGCCAAAGCACGCAAGCTTGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAAt                                                        <  2:1493534/145‑1 (MQ=255)
                                                                      aaGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAGCAGGTCGATCCGCGCAAAGCAGCAGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGAACAGCAGCCGGCTAAt                                                        >  2:1493594/1‑149 (MQ=14)
                                                                       aGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                                               >  1:1493658/1‑109 (MQ=255)
                                                                           gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTGGATCCGCGCAAAgc                                                                                                                     >  1:1493662/1‑83 (MQ=255)
                                                                             tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAgc                                                                                                                     >  1:1493700/1‑81 (MQ=255)
                                                                              aTTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAGAACGAGGTCGAGACGAGGAAAGCAGCCGTCGAGGACTCTATTGCCCGCGGCAAAGAAAGAAAACGGGAACAGCAACCGGCTAGTGCGGGGcc                                                >  2:1493665/1‑149 (MQ=9)
                                                                                 gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTAGCAGCGGCTATTGCCCGCGCAAAAGCCCGAAAACCCGCACAGCAACAGGCTGATGCCGAGcc                                                >  2:1493702/1‑146 (MQ=11)
                                                                                  cccccgCCAAAGCCCGCAAGCTTGAACACCAACAGCCTAATGCCTACCCAGAAGAACAGCGCGACCCGCGCTAAGACGCGTGCGAATCCGATATCGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCaga                                             <  1:1493599/143‑1 (MQ=16)
                                                                                  cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCAGTCGCCGCCGCTATTGCCCGCGCACAGGCCCGCAAACCCGCCCAGCAGCAGGCTAATGCCGAGCCcgaa                                            >  2:1493668/1‑149 (MQ=11)
                                                                                          aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAgcgg                                                                                                                   >  2:1493679/1‑68 (MQ=255)
                                                                                          aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAg                                                                                                                      >  1:1493680/1‑67 (MQ=255)
                                                                                           aaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                                                                                                                                            >  1:1493544/1‑44 (MQ=255)
                                                                                           aaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                                                                                                                                            <  2:1493544/44‑1 (MQ=255)
                                                                                            accACGCAAGCTGGAACAGCAACAGGCTAATGCCGATCCATAACACCAGGGCGTTCCTCGCATAGCCTCAGTCTAAGCCTCTATTGCCCGATTCATAGCTCGCAATCTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTg                                  <  2:1493545/147‑1 (MQ=14)
                                                                                                 gCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATcc                                                                                                                              >  1:1493686/1‑52 (MQ=255)
                                                                                                    aGCTGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg                                                                                                             <  1:1493547/66‑1 (MQ=25)
                                                                                                     gCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCAGTCGAAGCGGCTATTGCCCTCGACAAAGCGCGCAATCGTGAACAGCAACCGGCTAGTGAGGATACAGAAGAACAGGTTGATCCgcgcg                         >  2:1493692/1‑148 (MQ=14)
                                                                                                       tGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                      <  2:1493550/102‑1 (MQ=255)
                                                                                                       tGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                      >  1:1493550/1‑102 (MQ=12)
                                                                                                       tGGAACAGCAACAGGCTAATGCCGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                                                     <  1:1493548/71‑1 (MQ=21)
                                                                                                        ggAACAGCAACAGGCTAATGCCGAGCCAGAACAACAg                                                                                                                                      <  1:1493552/37‑1 (MQ=255)
                                                                                                        ggAACAGCAACAGGCTAATGCCGAGCCAGAACAACAg                                                                                                                                      >  2:1493552/1‑37 (MQ=255)
                                                                                                        ggAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATc                                                                                                                               >  1:1493613/1‑44 (MQ=255)
                                                                                                        ggAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCgcgc                                                                                          >  1:1493651/1‑81 (MQ=255)
                                                                                                          aaCAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc                                                                   <  1:1493554/102‑1 (MQ=255)
                                                                                                          aaCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc                                                                   >  2:1493554/1‑102 (MQ=12)
                                                                                                            cAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc                                                                                                                     >  1:1493698/1‑50 (MQ=255)
                                                                                                               cAACAGGCTAATGCCGAGCCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg                                                                                                 <  1:1493556/67‑1 (MQ=255)
                                                                                                                aaCAGGCTAATTCCGAGCCAGAACAACAGGTCGATCcgc                                                                                                                            <  1:3655954/39‑1 (MQ=255)
                                                                                                                  cAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGAGCGCAAACTGTAACAGAAAACTGATAATCAGGAGGAAGAACAAAAGCTTCATACACACAAAGCTGCAGTGGa            >  2:1493706/1‑149 (MQ=9)
                                                                                                                   aGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAc                                                                >  1:1493558/1‑96 (MQ=16)
                                                                                                                   aGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCTATCGGGAACAGCAACCGGCTAATGCGGGGCCAGATGAACAGGTTGATCCGCGCACAGCTGCCGTCgca           >  2:1493557/1‑147 (MQ=14)
                                                                                                                    ggCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGTGCCAGAAGAACAGGTTGATc                               >  1:1493622/1‑128 (MQ=255)
                                                                                                                     gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                                                         <  2:1493559/85‑1 (MQ=11)
                                                                                                                     gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                                                         >  1:1493559/1‑85 (MQ=11)
                                                                                                                      cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                                             <  2:1493558/96‑1 (MQ=255)
                                                                                                                       tAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa                                                                       <  1:1493560/85‑1 (MQ=255)
                                                                                                                       tAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa                                                                       >  2:1493560/1‑85 (MQ=255)
                                                                                                                       tAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc                                                                                                                     >  1:1493676/1‑39 (MQ=255)
                                                                                                                        aaTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGTGCAAGAAGAACAGGTTGATCAGCGAAAAGCTGCAGTAGAAGCgcc      >  2:1493736/1‑147 (MQ=14)
                                                                                                                           gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCg                                                                                             >  1:1493726/1‑59 (MQ=255)
                                                                                                                            ccGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAtt  >  1:1493694/1‑149 (MQ=255)
                                                                                                                             |                                                                                                                                                   
CGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATT  >  NC_000913/1708719‑1708991

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: