Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,816,5790GA41.5% 40.0 / 91.1 65intergenic (‑30/‑97)rph/yicCtruncated RNase PH/UPF0701 family protein YicC
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base A (13/14);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGGCA  >  NC_000913/3816532‑3816627
                                               |                                                
aCGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCt                                                >  1:1644746/1‑50 (MQ=255)
 cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTg                                               >  1:1317543/1‑50 (MQ=255)
 cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                                <  1:3322861/49‑1 (MQ=255)
  ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                              <  1:4947512/50‑1 (MQ=255)
    ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTgcg                                             >  1:340429/1‑49 (MQ=255)
    ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                             >  1:3188675/1‑49 (MQ=255)
      tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCAt                                          >  1:3689746/1‑50 (MQ=255)
      tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCAt                                          >  1:2086203/1‑50 (MQ=255)
        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCAt                                          <  1:1523636/48‑1 (MQ=255)
        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                          <  1:1279982/48‑1 (MQ=255)
        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                        <  1:4933239/50‑1 (MQ=255)
           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATtat                                       <  1:3619955/48‑1 (MQ=255)
           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATAc                                     >  1:4563396/1‑50 (MQ=255)
           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATAc                                     >  1:653583/1‑50 (MQ=255)
           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                     <  1:2240549/50‑1 (MQ=255)
            aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACg                                    <  1:580633/50‑1 (MQ=255)
            aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                    <  1:571628/50‑1 (MQ=255)
             cGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACgg                                   <  1:4350924/50‑1 (MQ=255)
             cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACgg                                   <  1:622853/50‑1 (MQ=255)
              gCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGa                                  >  1:2060058/1‑50 (MQ=255)
              gCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGa                                  <  1:1307037/50‑1 (MQ=255)
              gCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGa                                  <  1:3508960/50‑1 (MQ=255)
              gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                  >  1:462795/1‑50 (MQ=255)
              gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                  >  1:3407102/1‑50 (MQ=255)
               cATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGAc                                 >  1:542416/1‑50 (MQ=255)
               cATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGAc                                 <  1:1410817/50‑1 (MQ=255)
               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                                 >  1:4039646/1‑50 (MQ=255)
                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                                <  1:179509/50‑1 (MQ=255)
                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                                <  1:3977885/50‑1 (MQ=255)
                 tatTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACtt                               <  1:2165455/50‑1 (MQ=255)
                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                               <  1:5170343/50‑1 (MQ=255)
                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                               <  1:2623119/50‑1 (MQ=255)
                   ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                             <  1:645684/50‑1 (MQ=255)
                      aaaTCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCcggc                          <  1:2066243/50‑1 (MQ=255)
                      aaaTCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCcggc                          >  1:4836613/1‑50 (MQ=255)
                      aaaTCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCcggc                          <  1:582648/50‑1 (MQ=255)
                       aaTCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCcggcg                         <  1:3984862/50‑1 (MQ=255)
                       aaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcg                         >  1:5585833/1‑50 (MQ=255)
                        aTCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCcggcgg                        <  1:3893766/50‑1 (MQ=255)
                         tctcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                       >  1:3010226/1‑50 (MQ=255)
                          ctcCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGt                       >  1:3548549/1‑49 (MQ=255)
                           tcCGGCTTGAAACAGATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                     >  1:2608557/1‑50 (MQ=255)
                            ccGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTAt                    >  1:4806324/1‑50 (MQ=255)
                            ccGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAt                    >  1:2577059/1‑50 (MQ=255)
                            ccGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAt                    >  1:3470159/1‑50 (MQ=255)
                             cGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTa                     <  1:3648159/48‑1 (MQ=255)
                             cGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTAtt                   >  1:132897/1‑50 (MQ=255)
                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAt                    <  1:659679/49‑1 (MQ=255)
                              ggCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTc                  >  1:3112189/1‑50 (MQ=255)
                               gCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTcc                 >  1:247127/1‑50 (MQ=255)
                               gCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc                 >  1:2569298/1‑50 (MQ=255)
                               gCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc                 >  1:1028641/1‑50 (MQ=255)
                                cTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTcc                 >  1:390315/1‑49 (MQ=255)
                                cTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTcc                 >  1:3835787/1‑49 (MQ=255)
                                  tGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt              >  1:3054132/1‑50 (MQ=255)
                                   gAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt              <  1:1875294/49‑1 (MQ=255)
                                    aaaCAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATc             <  1:578669/49‑1 (MQ=255)
                                     aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt           <  1:357794/50‑1 (MQ=255)
                                     aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt           >  1:2658205/1‑50 (MQ=255)
                                      aCAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTg          <  1:784736/50‑1 (MQ=255)
                                        aaaTGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa         >  1:807173/1‑49 (MQ=255)
                                         aaTGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACa       <  1:3879966/50‑1 (MQ=255)
                                          aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc        >  1:5369526/1‑48 (MQ=255)
                                           tgtgGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg     >  1:5551173/1‑50 (MQ=255)
                                             tgGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGGc   <  1:1803935/50‑1 (MQ=255)
                                              ggCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGGCa  >  1:5192156/1‑50 (MQ=255)
                                               |                                                
ACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGGCA  >  NC_000913/3816532‑3816627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: