Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,528,4810CA20.0% 83.1 / 14.7 35Y671* (TAC→TAA) ‡yrfFputative RcsCDB‑response attenuator, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (14/14);  new base A (3/4);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGAATAACACGTTTGCG  >  NC_000913/3528340‑3528614
                                                                                                                                             |                                                                                                                                     
cGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCg                                                                                                                                >  6:72550/1‑149 (MQ=255)
 gAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGCTACTCAGCATATTGTCCTTCATCCGATCGCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGCTGCTGGGTAGCGCCATTGTAAATGGCgc                                                                                                                               <  7:31039/149‑2 (MQ=255)
  aGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTa                                                                                                                              <  3:92929/149‑1 (MQ=255)
               aaaaaTCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGcgtc                                                                                                                 >  6:108481/1‑149 (MQ=255)
               aaaaaTCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGcgtc                                                                                                                 <  5:108481/149‑1 (MQ=255)
               aaaaaTCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGc                                                                                                                    <  6:126606/146‑1 (MQ=255)
               aaaaaTCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGc                                                                                                                    >  5:126606/1‑146 (MQ=255)
                        cACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCa                                                                                                          <  5:72550/147‑1 (MQ=255)
                                                       cTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAAt                                                                                                                                    >  1:146095/1‑90 (MQ=255)
                                                       cTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAAt                                                                                                                                    <  2:146095/90‑1 (MQ=255)
                                                       cTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTa                                                                                                             >  6:41302/1‑113 (MQ=255)
                                                       cTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTa                                                                                                             <  5:41302/113‑1 (MQ=255)
                                                            tCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCtat                                                                     >  1:39989/1‑148 (MQ=255)
                                                                 tCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGa                                                                >  8:14134/1‑148 (MQ=255)
                                                                           ttCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGTGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGc                                                         <  3:32281/145‑1 (MQ=255)
                                                                           ttCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGAGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGc                                                         >  4:32281/1‑145 (MQ=255)
                                                                                       gcgcTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAAc                                            <  8:120974/146‑1 (MQ=255)
                                                                                       gcgcTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAAc                                            >  7:120974/1‑146 (MQ=255)
                                                                                                 aGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTc                               <  2:39989/149‑1 (MQ=255)
                                                                                                      cacaTTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGcc                                                        <  2:4576/119‑1 (MQ=255)
                                                                                                      cacaTTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGcc                                                        >  1:4576/1‑119 (MQ=255)
                                                                                                                  aCTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGaaa                          >  8:20670/1‑137 (MQ=255)
                                                                                                                  aCTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGaaa                          <  7:20670/137‑1 (MQ=255)
                                                                                                                   cTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGAAt              <  7:14134/148‑1 (MQ=255)
                                                                                                                             tGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAAc                                            >  8:50437/1‑108 (MQ=255)
                                                                                                                             tGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAAc                                            <  7:50437/108‑1 (MQ=255)
                                                                                                                              gggTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGAATAACACGTTTgcg  >  6:72225/1‑149 (MQ=255)
                                                                                                                               ctgtacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa               >  4:62122/6‑135 (MQ=255)
                                                                                                                               ctgtacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa               <  3:60031/130‑1 (MQ=255)
                                                                                                                               ctgtacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa               <  3:62122/130‑1 (MQ=255)
                                                                                                                               ctgtacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa               >  4:60031/6‑135 (MQ=255)
                                                                                                                                  tacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa               <  1:51461/130‑1 (MQ=255)
                                                                                                                                  tacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa               >  2:51461/3‑132 (MQ=255)
                                                                                                                                    cgcCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAgg                                                                          <  7:118801/71‑1 (MQ=255)
                                                                                                                                    cgcCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAgg                                                                          >  8:118801/1‑71 (MQ=255)
                                                                                                                                             |                                                                                                                                     
CGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGAATAACACGTTTGCG  >  NC_000913/3528340‑3528614

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