Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730587,9570GT74.2% 38.5 / 22.6 31A621S (GCG→TCG) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base T (9/14);  total (11/20)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAAACTCTAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGAACACGTTGCAG  >  CP000730/587821‑588092
                                                                                                                                        |                                                                                                                                       
aGCAAACTCTAAATTAGATGAACATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGAAAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGAc                                                                                                                                     <  2:332656/141‑1 (MQ=255)
aGCAAACTCTAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGAc                                                                                                                                     <  2:332651/141‑1 (MQ=255)
  cAAAATTTAAATTAGAAGAAAAAGGTCCTTTCATTGAAGAATAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAg                                                                                                                                   <  1:287014/141‑1 (MQ=255)
                gATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTc                                                                                                                     <  1:222602/141‑1 (MQ=255)
                                tCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACt                                                                                                     <  1:351985/141‑1 (MQ=255)
                                tCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTAAAAAAAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACt                                                                                                     <  2:118441/141‑1 (MQ=255)
                                            aaGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGAGACACCATCCATACCATACTTAAAAAACCATCACAAAAACCGTGCAt                                                                                         >  1:429478/1‑141 (MQ=255)
                                                   aTGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATggg                                                                                  <  2:394861/141‑1 (MQ=255)
                                                           tCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACa                                                                          <  2:401094/141‑1 (MQ=255)
                                                              gTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGc                                                                       >  1:417429/1‑141 (MQ=255)
                                                                  tAACAATACAGTTATGGCGAAAAAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACg                                                                   <  2:210489/141‑1 (MQ=255)
                                                                     cAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                                                <  2:50800/141‑1 (MQ=255)
                                                                          cAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATGCCATTCGTAGAAAATGAGGACTCAAACCGGGCATTGATGTGTGCGAACATGCAACGTCAAGCAg                                                           <  2:37785/141‑1 (MQ=255)
                                                                          cAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAg                                                           <  2:312979/141‑1 (MQ=255)
                                                                                tGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCtt                                                     <  2:6248/141‑1 (MQ=255)
                                                                                 ggCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCttt                                                    >  1:64498/1‑141 (MQ=255)
                                                                                      aaGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTtgatga                                               >  2:240283/1‑141 (MQ=255)
                                                                                             aaTGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGa                                        <  2:256117/141‑1 (MQ=255)
                                                                                                 gATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCa                                    >  2:269565/1‑141 (MQ=255)
                                                                                                  aTTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCAc                                   >  2:239162/1‑141 (MQ=255)
                                                                                                  aTTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCAc                                   <  1:171549/141‑1 (MQ=255)
                                                                                                          gATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATttgtt                           <  1:36813/141‑1 (MQ=255)
                                                                                                                tCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACa                     >  1:120252/1‑141 (MQ=255)
                                                                                                                       aGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATgg              >  1:263250/1‑141 (MQ=255)
                                                                                                                        gCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGa             <  2:179591/141‑1 (MQ=255)
                                                                                                                         cAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGaa            >  2:202577/1‑141 (MQ=255)
                                                                                                                         cAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGaa            >  2:372514/1‑141 (MQ=255)
                                                                                                                         cAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTTATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGCAGCACCATTTGTTGGTACAGGTATGGaa            <  2:429478/141‑1 (MQ=255)
                                                                                                                           aGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGAaca          <  1:76467/141‑1 (MQ=255)
                                                                                                                               gttTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGAACACGtt      <  1:132193/141‑1 (MQ=255)
                                                                                                                                 ttCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGAACACGTTgc    >  2:377692/1‑141 (MQ=255)
                                                                                                                                   cagcagCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGAACACGTTgcag  <  1:184579/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AGCAAACTCTAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGAACACGTTGCAG  >  CP000730/587821‑588092

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