Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U000961,708,7510AG23.2% 81.6 / 19.8 57E611E (GAA→GAG) rsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base A (16/27);  new base G (12/1);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 9.41e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAT  >  U00096/1708608‑1708894
                                                                                                                                               |                                                                                                                                               
tccgtcgcCAATGCTCGCAAGCGTTAACAGCAACTGTCTAATGCGTATCTCTTACATCAGGGCGATCCGGGAAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCaga                                                                                                                                            <  2:4025257/144‑1 (MQ=9)
     cgcCAAAGCGCGCAAGCTGGAACAGCAACAGCCAAATGCGGAACCAGCACACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAc                                                                                                                                       <  2:1639199/149‑1 (MQ=14)
      gcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAGGCAGCCGACGAGCCCGATATTGCCAGTGACTAAGCGCGCAAGCTGGCCCAGCATCAGGCTAAGGCGCAACCAGACGAGCa                                                                                                                                      >  2:4036157/1‑149 (MQ=14)
         aaaGCGCTCTAGCTTGATCAGCAACAGGCTAATTCGGAACCAGAACAACAGGCCGATCCTCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                                                                                   <  2:1639200/149‑1 (MQ=14)
         aaaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTAGATCCGCGAAAGGCAGCCGTCGACGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                                                                   >  2:1639295/1‑149 (MQ=11)
          acgctcgcAAGCTGGCACAGCAACAGGCAACTGCGGACCCAGCACCCCAGGCCGATCCTCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                                                                  <  2:1639203/144‑1 (MQ=14)
          aatctcgcAAGCTGGAACATCATCAGCATAATTCGGAACCAGAACAACAGGCCGATCCTCGCAAATCTGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                                                                  <  2:1639205/144‑1 (MQ=14)
               cggaaGCTGGAAGACCAACAGCAAAATCAGGACGCAGAACCAGAGGACGATCCGCGCAAATCCGCCGTCGAAGCCGCTATGGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATcc                                                                                                                             <  2:1639215/146‑1 (MQ=14)
                    ccTGGAACACCTTCAGGAAAATCCGGACGCAGAACACGAGGCCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCa                                                                                                                        <  2:1639218/148‑1 (MQ=14)
                      tGGAACCGCAACATGCTACTGCGCCACCAGACGACCAGGCCGATCCTCGCACAGCCGCCGCCGCCGCCGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCaaa                                                                                                                      <  2:4685744/149‑1 (MQ=11)
                        gAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGGAGAACAGGTCGATCCGCGCAAAgc                                                                                                                    >  1:3745731/1‑149 (MQ=255)
                         aaCATCAACAGGCTAATGCGGAACCAGAACCACAGGCCGATCCCCGCAAACCCGCCGTCGAAGCCGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                                                                   <  2:1639209/149‑1 (MQ=14)
                         aaCACCAACAGCCTACTCCGGACCCAGAACAACAGGCCGATCCGCGCAAATCAGCCGTCGAAGCCGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                                                                   <  2:1639208/149‑1 (MQ=14)
                                  aGGCTAATGCGCAACCAGAACACCCGGTCGATCCGCGCAAAGCCGCCGGCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGaa                                                                                                          <  2:2814873/149‑1 (MQ=255)
                                                       aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                                                                     <  2:1639212/149‑1 (MQ=255)
                                                          aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                                                                   >  1:1639213/1‑100 (MQ=37)
                                                             ccGATCCCCGCAAACCTGCCGTCGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGc                                                                               <  2:1639214/148‑1 (MQ=11)
                                                                 tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGATGGGCAGGTCGATCCGCGCTACGCCGCGGTCGGAGCCGCTATTGCCCGTGCCAACGCGCGCAAGc                                                                           >  2:1639221/1‑149 (MQ=14)
                                                                         aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGtt                                                                                                                                  >  2:1639332/1‑85 (MQ=255)
                                                                              gccgTTGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                                                                    <  1:1639224/79‑1 (MQ=255)
                                                                              gccgTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                                                                    <  1:1639227/79‑1 (MQ=255)
                                                                              gccgTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                                                                  <  1:1639225/81‑1 (MQ=255)
                                                                                 gTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGc                                                                                         <  1:1639234/119‑1 (MQ=255)
                                                                                 gTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGc                                                                                         >  2:1639234/1‑119 (MQ=255)
                                                                                  tggatgCCGCTATTGCCCGCGCCATAGCGCGCAAGCTGGAACATCATCAGGATAATCCCGAACCAGAACAACAGGGCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCt                                                          <  1:1639244/144‑1 (MQ=11)
                                                                                  tCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACa                                                                                                                                      >  1:1639341/1‑73 (MQ=255)
                                                                                    gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcgcg                                                                                                                          <  2:1639421/83‑1 (MQ=255)
                                                                                    gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCcgcg                                                                                                                          >  1:1639421/1‑83 (MQ=255)
                                                                                    gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCg                                                                                                                                 <  1:1639240/76‑1 (MQ=255)
                                                                                     aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCa                                                                                                                        >  2:1639343/1‑84 (MQ=255)
                                                                                     aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                                                                                   <  1:1639242/73‑1 (MQ=255)
                                                                                       tttgCTATAGCCCGTGCCAAAGCGCGAAAGCTGGAACAGCAACAGGCTAATTCGGAACCAGAAGAACAGTTCTATCCGCTCAAAGCCTCCGCCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGc                                                     <  2:1639246/146‑1 (MQ=11)
                                                                                       gCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAgc                                                                                                                    >  1:1639433/1‑86 (MQ=255)
                                                                                            tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                                                                   >  1:1639354/1‑66 (MQ=255)
                                                                                            tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                                                                   <  2:1639354/66‑1 (MQ=255)
                                                                                             atcgCCTGTGCCAAAGCGCTCATCCAGGACCACCAACACCCCAATGCGGACCCATTCGCACCTCTCGTTCCGCGCCCAGCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGcc                                               <  2:4713772/146‑1 (MQ=9)
                                                                                             aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCaga                                                                                                                                            >  1:1639356/1‑56 (MQ=255)
                                                                                             aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCAGCCGt                                                                                                              >  2:1639355/1‑86 (MQ=255)
                                                                                                ggccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTACTGCGGAACCAGAAGACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                                            <  1:4539726/147‑1 (MQ=255)
                                                                                                   cGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCaaa                                                                                                                      >  1:1639286/1‑72 (MQ=39)
                                                                                                     tGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAgg                                                                                                                                    >  1:1639316/1‑56 (MQ=255)
                                                                                                     tGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCTGCCGt                                                                                                              >  2:1639362/1‑78 (MQ=11)
                                                                                                        cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGAt                                                                                                                               >  2:1639365/1‑58 (MQ=11)
                                                                                                          aaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTc                                  <  2:1639249/149‑1 (MQ=255)
                                                                                                            gcgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACa                                                                                                                                      >  1:1639368/1‑47 (MQ=255)
                                                                                                              gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                                                                                            >  2:1639369/1‑55 (MQ=11)
                                                                                                                     cTGCAACAGCAACATTCTAATGGCTAACGAGAAGACCTTGTCGTTCCGCGCAAACCCCCCGGCCTATCTTCTTTTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaa                       <  2:1639286/149‑1 (MQ=9)
                                                                                                                      tGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                                                                    >  1:1639372/1‑39 (MQ=40)
                                                                                                                      tGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGt                                                                                                              <  1:2362326/61‑1 (MQ=35)
                                                                                                                      tGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcgagcc                      >  2:3753700/1‑145 (MQ=255)
                                                                                                                       ggAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                                           <  1:1639251/79‑1 (MQ=34)
                                                                                                                           cAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAt                                                                                                  >  1:1639373/1‑68 (MQ=255)
                                                                                                                            aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGcac                                                                                             >  1:1639374/1‑70 (MQ=255)
                                                                                                                                   ggCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAg                                                                                                                     >  1:1639503/1‑41 (MQ=38)
                                                                                                                                     cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgc                                                                                                                    >  1:1639387/1‑40 (MQ=40)
                                                                                                                                     cTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCaa                                                                                                                       >  2:1639385/1‑37 (MQ=255)
                                                                                                                                          gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCGGCTATTGCCCGAGCCAACGCACGCAAGCTGGAAAAGCAACAGGCAAATGCGGCGCCAGAACTACAGGTGGATCCGCGCACAGCGGCCGTTGAAGCCGCTAt  >  1:1639399/1‑149 (MQ=11)
                                                                                                                                               aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGcc                                               >  1:4495635/1‑99 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAT  >  U00096/1708608‑1708894

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