Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U000961,708,7090TA27.5% 36.8 / 23.6 57R597R (CGT→CGA) rsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base T (8/29);  new base A (6/8);  total (14/43)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAA  >  U00096/1708567‑1708840
                                                                                                                                              |                                                                                                                                   
ccGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                                                                                                               <  2:1067306/149‑1 (MQ=255)
  gTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgc                                                                                                                             <  1:1067307/149‑1 (MQ=255)
   ttGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTCGAACAGAAACAGAATAATGAAGAACAAGAAGAACAGCTAGAACATAAAAAAGAAGCCGTCGAAGCCGCTATTGCCCGTGCCTGAgcg                                                                                                                            >  2:1067373/1‑149 (MQ=12)
      cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcgc                                                                                                                         <  2:1067308/149‑1 (MQ=255)
       tGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCa                                                                                                                        <  2:1067309/149‑1 (MQ=255)
        gATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCaa                                                                                                                       <  1:1067312/149‑1 (MQ=255)
           ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGc                                                                                                                     <  1:1067314/148‑1 (MQ=255)
            cACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGCAACCAGAACAACAGGTCGATCCGCGCAAAGACGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGAGCGCAAGCTg                                                                                                                   >  1:1067363/1‑149 (MQ=255)
               gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGaa                                                                                                                <  2:1067315/149‑1 (MQ=255)
                tAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                                                                 <  1:1067316/147‑1 (MQ=255)
                 aaaaCTGCCGCTGAAGCAGCTAGCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAACCACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACa                                                                                                              <  2:1067317/149‑1 (MQ=255)
                        ccGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAg                                                                                                       <  2:1067322/149‑1 (MQ=255)
                            tGACGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTa                                                                                                   <  2:1067323/149‑1 (MQ=255)
                             gAATCAGCTATCGCCTGCGCTAAATCTCGGAATCTGGAAGATCAACAGGCTAATGCGGAAGCAGAAAAAGAGGCCGATCCTCGCAAATCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                                                                                  <  2:1067324/149‑1 (MQ=14)
                                gcagcTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGTAACCAGAACAACATGTCTATCCGCGCATATCCGCAGCCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGc                                                                                               <  2:1067326/149‑1 (MQ=11)
                                  agcTATCTCCCGCGCCTAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCgg                                                                                             <  2:1067327/149‑1 (MQ=255)
                                     tATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGc                                                                                                                                   <  2:1067328/108‑1 (MQ=255)
                                     tATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGc                                                                                                                                   >  1:1067328/1‑108 (MQ=255)
                                            cgcgCCAACTCCCGGAACCAGGAACAGCATCAGCATACTGCCGAAGCAGCACAACAGGCCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATAGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaaga                                                                                   <  2:1067337/149‑1 (MQ=11)
                                               gacAAATCGCGGAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                                                                <  2:1067329/147‑1 (MQ=255)
                                                 atcctcccgcAAGCTGGAACAGCTTCAGCATAATTCCGAACCAGCACAAGTGGCCGATCCTCGCAAACCAGCCGTCGAAGCAGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                              <  2:1067339/142‑1 (MQ=14)
                                                     cgtcggaaGCTGGAACACCTTCAGCATAATGCGGACCCAGCACAAGAGGGCGATCCGCGGAAATCTGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa                                                                          <  2:1067335/143‑1 (MQ=14)
                                                        cggaaGCAGGAAGAGCAACAGGATACTCCGGAAGCAGCACCCCTTGGCGATCCTCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATcc                                                                       <  2:1067342/146‑1 (MQ=14)
                                                        cgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGaa                                                                                                                <  1:1067331/108‑1 (MQ=35)
                                                                ggAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCa                                                                                                                        >  1:1067392/1‑92 (MQ=18)
                                                                   aCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                                                                                                               <  1:1067393/82‑1 (MQ=18)
                                                                   aCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                                                                                                               >  2:1067393/1‑82 (MQ=18)
                                                                          cAGGCAACTCCGGAACCAGCACCACCGGGCGATCCCCGCAAAGCCGCCGTCGAAGCCGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGa                                                     <  1:1067340/149‑1 (MQ=14)
                                                                            ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                                                                                                               <  1:1067397/73‑1 (MQ=21)
                                                                            ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                                                                                                               >  2:1067397/1‑73 (MQ=21)
                                                                              cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCa                                                                                                                        >  1:1067399/1‑78 (MQ=21)
                                                                              cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCa                                                                                                                        <  2:1067399/78‑1 (MQ=21)
                                                                                 ctcacgAACCAGAAAAACATGCCGATCCTAGAAAATCCGCCTTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATTCGGAACCAGAACAACAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCt                                              <  2:1067352/144‑1 (MQ=11)
                                                                                       aaCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAAc                                                                                 >  1:1067402/1‑108 (MQ=255)
                                                                                                   aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCa                            <  2:1067338/149‑1 (MQ=255)
                                                                                                      ccGACCCCCGCAAACCCGCCGTTTAAGCACCTCTTCCCCGGGCCACAGCGCGCAAGCTGGAACACCAACACGACTATCCTGATACAGCACAAGAGGCCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGc                         <  2:1067367/148‑1 (MQ=2)
                                                                                                       cGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg                                                                               <  1:1067344/94‑1 (MQ=255)
                                                                                                        gATCCTCGCAAATCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCg                                                      <  1:1067347/118‑1 (MQ=255)
                                                                                                        gATCCTCGCAAAGCCGCCGTCGAAGCAGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                             <  1:1067349/95‑1 (MQ=255)
                                                                                                        gATCCTCGCAAACCCGCCGTCGAAGCCGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                             <  1:1067346/95‑1 (MQ=14)
                                                                                                        gATCCGCGCAAACCCGCCGTCGAAGCCGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                                                                                  <  1:1067348/74‑1 (MQ=11)
                                                                                                           ccGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCt                    <  1:1067350/149‑1 (MQ=255)
                                                                                                              cgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAAt                                                                                                 >  2:1067419/1‑69 (MQ=25)
                                                                                                                  aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGtag                                                                           >  1:1067422/1‑85 (MQ=255)
                                                                                                                       gccgTTGAAGCCGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccgc                                                           <  1:1067361/98‑1 (MQ=16)
                                                                                                                       gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAgg                                                                              >  1:1067423/1‑79 (MQ=255)
                                                                                                                       gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                          <  1:1067353/67‑1 (MQ=255)
                                                                                                                       gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                                                                                        <  1:1067359/69‑1 (MQ=255)
                                                                                                                       gccgTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                                                                <  1:1067356/77‑1 (MQ=255)
                                                                                                                       gccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                            <  1:1067357/81‑1 (MQ=255)
                                                                                                                       gccgTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                             <  1:1067360/80‑1 (MQ=255)
                                                                                                                       gccgTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAc                                                                                 <  1:1067354/76‑1 (MQ=255)
                                                                                                                          gTCGAAGCAGCTATTGCCCGTGCCATAGCGCGCAAGCTGGAACAGCATGAGGCTAATGCTGAAGCATAAGATCTGGTCGATCCGCGTAAATCAGCCGTCGAAGCCGCTATTGCTCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGc     <  2:1067362/149‑1 (MQ=14)
                                                                                                                           tCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCgg                                                                                             >  1:1067473/1‑60 (MQ=255)
                                                                                                                            cGAAGCCGCTATTGCCCTTGCCAAAGCGCGCAAGCTTGAACAGCAACAGGTTTATGCGGAACCTTATCAACAGGTCGCTCCTCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTa   <  2:1067363/149‑1 (MQ=11)
                                                                                                                             gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTCGAAGACGCTATAGCCAGGGCCAAAGAAAGCAAGCTGGAAACGCAACAGGCGaa  >  2:1067444/1‑149 (MQ=14)
                                                                                                                                         gCCCGTGCCAAAGCGCGCAAGCGGGAACAGCAAcc                                                                                                        >  1:1067432/1‑34 (MQ=255)
                                                                                                                                              |                                                                                                                                   
CCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAA  >  U00096/1708567‑1708840

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