Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,628,2130CA20.4% 92.4 / 13.1 44A442A (GCG→GCT) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/10);  new base A (1/8);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.67e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAACGTCGTTGAGCTTAAAACG  >  NC_000913/3628067‑3628359
                                                                                                                                                  |                                                                                                                                                  
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)
                                                                                      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                                                     >  6:65749/1‑148 (MQ=255)
                                                                                              tCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGcc                                                     >  2:113145/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)
                                                                                                                                      gAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCg                                    <  7:197254/125‑1 (MQ=255)
                                                                                                                                      gAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCg                                    >  8:197254/1‑125 (MQ=255)
                                                                                                                                         cATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATAACTTCAACGTCGTTGAGc          >  2:231909/1‑148 (MQ=255)
                                                                                                                                          atatCCCTAGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCAGGCGAATGCCGAGCGGCAATGCCTCCGGCAGAATATCTTCAACGTCGTTGAGca         >  1:113560/1‑147 (MQ=255)
                                                                                                                                            atCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAACGTCGTTGAGCTTa       >  1:129235/1‑148 (MQ=255)
                                                                                                                                             tCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAACGTCGTTGAGCTTaa      >  1:50564/1‑148 (MQ=255)
                                                                                                                                              cccTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAACGTCGTTGAGCTTaaaa    >  7:146300/1‑149 (MQ=255)
                                                                                                                                                cTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAACGTCGTTGAGCTTAAAAcg  >  7:232294/1‑149 (MQ=255)
                                                                                                                                                 tCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAACGTCGTTGAGCTTAAAAcg  >  6:258428/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGGAGATGAAGATAGTCACTTTGTCCTGGCGCGAGAGATCGACCATCAACTGCCAGAACATATCCCTCGCCACCGGATCGACACCAGAAGTAGGCTCATCGAGGATTAACATCTCCGGGCGATGAATCACCGCCACCGCCAGCGAAAGCCGCTGGCGAATGCCGAGCGGCAATGACTCCGGCAGAATATCTTCAACGTCGTTGAGCTTAAAACG  >  NC_000913/3628067‑3628359

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