Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,628,1990CT20.0% 85.5 / 11.3 41W447* (TGG→TAG) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base C (19/13);  new base T (0/8);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 3.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAACG  >  NC_000913/3628057‑3628342
                                                                                                                                              |                                                                                                                                               
cAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACa                                                                                                                                           <  7:239363/149‑1 (MQ=255)
  gTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACata                                                                                                                                         >  7:146789/1‑149 (MQ=255)
  gTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACata                                                                                                                                         <  1:312520/149‑1 (MQ=255)
      cTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATcc                                                                                                                                      <  1:374838/148‑1 (MQ=255)
          ggTGTACCGCTGGCAAGCACTTTTCCGTCGTGCATCATTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTTGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGc                                                                                                                                 <  5:337507/149‑1 (MQ=255)
              tACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCAcc                                                                                                                             >  8:184698/1‑149 (MQ=255)
                 cGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGa                                                                                                                          >  8:339288/1‑149 (MQ=255)
                        aaGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGacac                                                                                                                   <  7:184698/149‑1 (MQ=255)
                         aGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGacac                                                                                                                   >  8:291889/1‑148 (MQ=255)
                          gCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGacac                                                                                                                   >  8:350574/1‑147 (MQ=255)
                            aCTTTTCCGGCGTGCAACAGTGAGATGCGGTCGCAACGTTCCGCTTCGTACATAAAGTCGGAGTAGATTAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGAGCATCAACTGCAAGAACATATCCCTCGCCACCGGATCGACACCAg                                                                                                                >  4:45824/1‑148 (MQ=255)
                              ttttCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGGGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGCTCGACCATCAACTGCCAGAACATATCCCTCGCCCCCGGCTCGACACAAGact                                                                                                             >  3:233103/1‑147 (MQ=255)
                                   cGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGc                                                                                                        >  3:115264/1‑149 (MQ=255)
                                           aTCAGTGAGATGCGGTCGCACCGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCCTCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAg                                                                                                >  1:75218/1‑149 (MQ=255)
                                           aTCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAg                                                                                                >  8:286620/1‑149 (MQ=255)
                                             cAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATc                                                                                                                        <  8:354514/123‑1 (MQ=255)
                                             cAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATc                                                                                                                        >  7:354514/1‑123 (MQ=255)
                                              aGTGAGATGCGGGCGCAACGTTCCGCTTCGGTCATAAAGGGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATTTTTTTTTTTATGATATGGCTAGCAACGTGATCTACACCAGAAGTAGGCTCATCGAGGAt                                                                                             <  6:71004/149‑1 (MQ=255)
                                                  agaTGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAAc                                                                                         <  7:339288/149‑1 (MQ=255)
                                                   gaTGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGAt                                                                                                                         >  5:277648/1‑116 (MQ=255)
                                                   gaTGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGAt                                                                                                                         <  6:277648/116‑1 (MQ=255)
                                                     tGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATTGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACAtc                                                                                      <  3:192346/149‑1 (MQ=255)
                                                        ggTCGCAACGTTCCGCTTCGTTCATACAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATTTTTTTTTTTATGATATGGCTATCAACGTGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTcc                                                                                   <  2:254961/149‑1 (MQ=255)
                                                        ggTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTcc                                                                                   <  7:286620/149‑1 (MQ=255)
                                                        ggTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACAATTTTTTTTTTTATGATATGGCTAGCAACGTGATCTACACCAGAAGTAGGCTCATCGAGGATTAACATCTcc                                                                                   <  2:96888/149‑1 (MQ=255)
                                                                     cGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATTTTTTTTTTTATGATATGGCTAGCAACGTGATCTACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATc                                                                       <  4:22409/148‑1 (MQ=255)
                                                                        ttCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATcaccg                                                                   >  3:12099/1‑149 (MQ=255)
                                                                          cggtCATAAAGTGGGGGGAGATGACGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCTTTTTTTGTTTTATGATATGGCTAGCAACGTGATCTACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATcaccgc                                                                  <  2:75218/145‑1 (MQ=255)
                                                                           gTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGACTcaccgcca                                                                >  2:71804/1‑149 (MQ=255)
                                                                                    gtgggtggAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACatat                                                                                                                                        <  4:335535/68‑1 (MQ=255)
                                                                                    gtgggtggAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACatat                                                                                                                                        >  3:335535/1‑68 (MQ=255)
                                                                                                gAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGCAGTAGGCGCATCGAGGATTAACATATCCGGGCGATGAATcaccgccacc                                                              >  6:381008/1‑130 (MQ=255)
                                                                                                gAAGATAGGCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATcaccgccacc                                                              <  5:381008/130‑1 (MQ=255)
                                                                                                 aaGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATTTTTTTTTTTATGATATGGCTAGCAACGTGATCTACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGc                                           <  8:146789/148‑1 (MQ=255)
                                                                                                     tAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATg                                      <  7:350574/149‑1 (MQ=255)
                                                                                                       gTCACTTTGTCCTGGCGCGAGAGATCGACCATTTATTTTTATATCATATGGCTAGCAACGTGATCTACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGc                                     <  4:115264/148‑1 (MQ=255)
                                                                                                        tCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGcc                                    >  2:113145/1‑148 (MQ=255)
                                                                                                        tCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGcc                                    >  6:65749/1‑148 (MQ=255)
                                                                                                          aCTTTGTCCTGGCGCGAGAGATCGACCATCATTTTTTTTAAGATATGGCTAGCAACGTGATCTACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGa                                  <  4:12099/148‑1 (MQ=255)
                                                                                                                  cTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATAGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAAt                          >  6:299676/1‑148 (MQ=255)
                                                                                                                                         aCTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAAcg  >  3:186291/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAACG  >  NC_000913/3628057‑3628342

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