Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U000961,708,7510AG32.4% 34.8 / 21.8 37E611E (GAA→GAG) rsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base A (5/18);  new base G (9/2);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 1.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCC  >  U00096/1708609‑1708867
                                                                                                                                              |                                                                                                                    
cccTCGCCAAAGCGAGCAAGCTGGAACAGCATCAGGCTAATTCGGAAGCAGAACAACAGGTCGATCCTCGCAAGTCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaa                                                                                                                <  2:1050323/149‑1 (MQ=14)
         atgctcgcAAGCTGGAAGAGCAACAGGATAATGAGGAACCAGAACACCTGGTCGATCCTCGTAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                                       <  2:1050328/144‑1 (MQ=14)
          aacccgcAAGCAGGAACATCAACAGGAAAATTAGGAACCAGCACAAGATGTCGATCCTCGCAAACCAGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCg                                                                                                      <  2:1050331/145‑1 (MQ=14)
              cgcAATCTGGAAGATCATCAGGATAATTCGGAACCAGAAAAACATGCCGATCCTCGCAAAGCCGCCGTCGAAGCAGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa                                                                                                     <  2:1050325/146‑1 (MQ=14)
                cAAGCTGGAACAGCATCAGGCTAATTAGGAACCAGAACACCATGCCGATCCGCGGAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc                                                                                                <  1:1050330/149‑1 (MQ=14)
                           aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGCGAGCAAGATGGAACCGCAGCAGGAGACTGAGGAGCCAGAAGCAAAGATAGATCCGCGCAAAgcggcc                                                                                     >  2:1050391/1‑149 (MQ=11)
                                         gccgAACCAGAACCCCTGGCTGATCCTCGAAAAACAGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt                                                                       <  2:1050344/146‑1 (MQ=11)
                                           ggAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGCCaga                                                                                                                 >  2:1050393/1‑105 (MQ=255)
                                                 aGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGAGCGCAACCTGGAACAGAAACAGGCAAATGCCGATCCAGAAGATCAGGTAGATCAGCTCAAAGCGGCCGTCGAAGCGGCTATTGCCCGTg                                                               >  2:1050444/1‑149 (MQ=11)
                                                    acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCa                                                            <  2:1050332/149‑1 (MQ=255)
                                                      aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGAAACAGGCTAATGCGGTACCAGCAGAACAGCTCGATCCGCGCAAAGCCGCCGTCGCAGCAGCGATTGCCCGTGCCaaa                                                          >  1:2600368/1‑149 (MQ=14)
                                                       acaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACa                                                                                                           >  1:1050446/1‑99 (MQ=255)
                                                              gATCCTCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTGATGCGGAACCAGAAGAACa                                                                                                           <  1:1050338/92‑1 (MQ=255)
                                                              gATCCGCGCAAAGCAGCCGTCGAAGCCGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc                                                                                                <  1:1050336/103‑1 (MQ=255)
                                                                     gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGcc                                                                                                                    >  1:1050410/1‑76 (MQ=255)
                                                                     gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCaga                                                                                                                 >  1:1050409/1‑79 (MQ=21)
                                                                         aGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGTAGACAGAAGAACCGGGGGAGCCGCGCAAAGCCGCCGTCTCAGCAGCTCTTTCCAGTGAGATAGCACTCAAGATCGAAAAGc                                       >  1:1286488/1‑149 (MQ=2)
                                                                             gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                                                     <  1:2881422/67‑1 (MQ=255)
                                                                             gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt                                                                                                    <  1:1050346/84‑1 (MQ=255)
                                                                             gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt                                                                                                    <  1:1050347/84‑1 (MQ=255)
                                                                             gccgTCGAAGCCGCTATCGCCCGTGACAAAGCGCGCAAGCTGGAACAGCAACATTCTAATGCGGAACCCCATCAACAGTTCGTTCCTCGTAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACa                                   <  2:1050348/149‑1 (MQ=11)
                                                                             gccgTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa                                                                                                     <  1:1050345/83‑1 (MQ=255)
                                                                                    aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCaga                                                                                                                 >  1:1050418/1‑64 (MQ=255)
                                                                                         gCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                                        <  1:1050351/84‑1 (MQ=18)
                                                                                         gCTATAGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGccgccg                                                                                    <  1:1050352/88‑1 (MQ=255)
                                                                                                cccGTGCCAAAGCGCTCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATc                                                                                                   <  2:1050426/66‑1 (MQ=14)
                                                                                                cccGTGCCAAAGCGCACAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATc                                                                                                   >  1:1050426/1‑66 (MQ=14)
                                                                                                 ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAg                                                                                                          >  1:1050429/1‑58 (MQ=255)
                                                                                                  cGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGa                                                                                <  1:1050428/83‑1 (MQ=11)
                                                                                                  cGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCTGCCGTCGa                                                                                >  2:1050428/1‑83 (MQ=11)
                                                                                                       cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                             <  1:1050354/97‑1 (MQ=18)
                                                                                                            cgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAAc                                    <  1:1050355/117‑1 (MQ=35)
                                                                                                              cgcAAGCTGGAACAGCAACATGATAATGATTATCCTGAACAAGAGTTCTATCCTCGTAAAGACGACGTCGATGCAGCTATTTCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc  <  2:1050376/149‑1 (MQ=11)
                                                                                                              cgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgctgc                                                                                      >  1:1050498/1‑65 (MQ=25)
                                                                                                                           aGCAACAGGCTAATGCGGATCCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTc                                                                                  <  1:1050436/56‑1 (MQ=11)
                                                                                                                                     tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGCCGCGGCTATTGCCCg                                                                 >  1:1050514/1‑63 (MQ=11)
                                                                                                                                            gAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGaa                                                                               >  1:1050442/1‑42 (MQ=255)
                                                                                                                                              |                                                                                                                    
CCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCC  >  U00096/1708609‑1708867

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: ‑