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

AGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGC  >  U00096/1708562‑1708848
                                                                                                                                                   |                                                                                                                                           
agcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                                                                                             <  2:1162891/148‑1 (MQ=255)
  caACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGc                                                                                                                                           <  2:1162892/148‑1 (MQ=255)
     tcGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaa                                                                                                                                        <  1:1307994/147‑1 (MQ=255)
      cGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                                                                                                                       <  2:1162894/148‑1 (MQ=255)
      cGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                                                                                                                       <  2:1162893/148‑1 (MQ=255)
       gTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAg                                                                                                                                      <  1:1162896/148‑1 (MQ=255)
         tGCTGATCCACGTAAAACTGCCGGTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcg                                                                                                                                    <  2:1162899/148‑1 (MQ=255)
           cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                                                                                  <  2:1162901/148‑1 (MQ=255)
           cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                                                                                  <  2:1162903/148‑1 (MQ=255)
           cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                                                                                  <  2:1162902/148‑1 (MQ=255)
           cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                                                                                  <  2:1162900/148‑1 (MQ=255)
            tGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcgc                                                                                                                                 <  1:1162906/148‑1 (MQ=255)
            tGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcgc                                                                                                                                 <  2:1162907/148‑1 (MQ=255)
               tCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                                                                                              <  2:1162908/148‑1 (MQ=255)
                  aCGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACGGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                           <  2:2072124/148‑1 (MQ=255)
                  aCGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                           <  2:1928183/148‑1 (MQ=255)
                  aCGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                           <  2:1928182/148‑1 (MQ=255)
                   cgccaaccTGCCTTTTAAGAAGCTATGGCCCGCGCCAACGCGAGCAAGCTGGAACAGCAACAGGCTAATCCGTACCCAGAATAAGAGGCAGATCCCCGCAAATCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                                                                                                                          <  2:1162946/141‑1 (MQ=2)
                     tAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGATCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGaa                                                                                                                        <  2:1162911/148‑1 (MQ=255)
                     tAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGaa                                                                                                                        <  1:1162909/148‑1 (MQ=255)
                     tAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGaa                                                                                                                        <  1:1162910/148‑1 (MQ=255)
                      aaatCTGCCGTTTAAGCATCTTTCGCCCGCGCCAAAGCTCGCAATCTGGAACAGCTACAGGCTAATGCGGAAGCAGAACACCTTGCCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAAc                                                                                                                       <  2:1162912/148‑1 (MQ=14)
                            gCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCACAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAAc                                                                                                                 <  2:1162916/148‑1 (MQ=255)
                                  gAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAACAGAACAACATGTCGATCCTCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTa                                                                                                           <  1:1162917/148‑1 (MQ=255)
                                   acgaagcTATCTCACGCGCCTACTCTCGCAATCTGGAACAGCAACAGGCTAATTCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                                                                                          <  2:1162925/144‑1 (MQ=11)
                                     gcagcTAGCTCCGGCGCCTAAGCGCGCAAGCTGGAACAGCATCAGGATAATGAGGAAGCAGAACCCCAGGGCGATCCTCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATg                                                                                                        <  2:1162921/148‑1 (MQ=14)
                                       agcTATCGCCCGCGCCACACCTCGCAAGCTGGAACAGCATCAGGATAATGCGGAACCAGCACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                                                                                                      <  1:1162922/148‑1 (MQ=14)
                                                     ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCa                                                                                                                                         <  1:1162974/99‑1 (MQ=18)
                                                     ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCa                                                                                                                                         >  2:1162974/1‑99 (MQ=18)
                                                     aatcctggcgcAACCAGGAACACCATCAGGATAATGCGGAACCAGAAACAGAGGCTGATCCTCGCAAATCAGCCGTCGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                                                                        <  2:1162914/141‑1 (MQ=11)
                                                       tcacctcgcAATCTGGCAGATCATCAGGAAAATCCCGAAGCAGAACAAGAGGTAGATCCGCGTAAACCCGCCGTCGAAGCAGCTATCGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                      <  2:1162933/142‑1 (MQ=9)
                                                        taGCGCGTAAGCTTGAAGATCTTCAGGATAATTAGGAAGCAGCAACAGATGTTGATCCTCGCAAATCTGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                                     <  2:1162943/147‑1 (MQ=9)
                                                        aatctcgcAAGCTGGAACATCATCAGGATAATTAGGAACCAGAACAACAGGTCGATCCTCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                                     <  2:1162926/143‑1 (MQ=14)
                                                                      gAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAg                                                                       <  1:1162929/148‑1 (MQ=255)
                                                                              aCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaacg                                                                                                                            >  2:1162982/1‑84 (MQ=0)
                                                                              aCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaacg                                                                                                                            <  1:1162982/87‑4 (MQ=0)
                                                                                 ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaa                                                                                                                               <  1:1162986/81‑1 (MQ=18)
                                                                                 ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaa                                                                                                                               >  2:1162986/1‑81 (MQ=18)
                                                                                  gCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                                                                                              <  1:1162987/81‑2 (MQ=17)
                                                                                  gCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                                                                                              >  2:1162987/1‑80 (MQ=17)
                                                                                                             gATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAg                                                                                                               <  1:1162937/69‑1 (MQ=25)
                                                                                                             gATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAgaagaa                                                                                          <  1:1162938/90‑1 (MQ=255)
                                                                                                                     cAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTACAGCGGAACCCCATCAACAGGTCGATCCGCGCAAAGCTGCCGTCGAGGCCGCTAATGAAAGAGACAAAGCACGCGCGCTggga                        >  2:1163058/1‑146 (MQ=11)
                                                                                                                            gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                <  1:1162940/85‑1 (MQ=255)
                                                                                                                            gccgTCGAAGCAGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa                                                                                  <  1:1162942/83‑1 (MQ=255)
                                                                                                                            gccgTAGAAGCCGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                <  1:1162941/85‑1 (MQ=11)
                                                                                                                                 cGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                                                                       >  1:1163052/1‑73 (MQ=255)
                                                                                                                                   aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                         >  1:2041445/1‑101 (MQ=255)
                                                                                                                                    aGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAAc                                                                                         >  1:1163005/1‑68 (MQ=255)
                                                                                                                                    aGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAAc                                                                                         <  2:1163005/68‑1 (MQ=255)
                                                                                                                                    aGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                      <  1:1162945/71‑1 (MQ=255)
                                                                                                                                        gCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTa                                                                                                           <  1:2339953/46‑1 (MQ=255)
                                                                                                                                           aTTGCCCGTGCCAAAGCGCGCAAGCTTGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGc  <  1:1686069/148‑1 (MQ=255)
                                                                                                                                               cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCg                                              >  1:1162988/1‑100 (MQ=37)
                                                                                                                                                 cGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt                                                    >  2:1162972/1‑92 (MQ=37)
                                                                                                                                                   |                                                                                                                                           
AGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGC  >  U00096/1708562‑1708848

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