Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,4950AT23.1% 46.8 / 11.4 26L15Q (CTG→CAG) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base A (10/10);  new base T (3/3);  total (13/13)
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.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTT  >  NC_000913/3629355‑3629621
                                                                                                                                             |                                                                                                                              
gggCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGTTCCCGTCCCGCCTGTcgcg                                                                                                                          >  3:37683/1‑148 (MQ=255)
      cGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGTTCCCGTCCCGCCTGTCGCGCAACTgg                                                                                                                   <  4:37683/149‑1 (MQ=255)
        gAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCTGTCGCGCAACTGGcgg                                                                                                                >  3:10016/1‑149 (MQ=255)
          aaTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCTGTCGCGCAACTGGcgggc                                                                                                              <  4:10016/149‑2 (MQ=255)
                      caaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                  >  6:5120/1‑149 (MQ=255)
                          cTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCTGTCGCGCAACTGGcgggcgt                                                                                                            >  6:27698/1‑132 (MQ=255)
                          cTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCTGTCGCGCAACTGGcgggcgt                                                                                                            <  5:27698/135‑4 (MQ=255)
                            cGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAt                                                                                            <  3:70975/149‑1 (MQ=255)
                             gACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGc                                                                                                                  <  4:30417/126‑1 (MQ=255)
                             gACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGc                                                                                                                  >  3:30417/1‑126 (MQ=255)
                             gACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGa                                                                                             <  5:5120/147‑1 (MQ=255)
                             gACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATg                                                                                           >  8:1224/1‑149 (MQ=255)
                                      aCGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGc                                                                                  >  1:6523/1‑149 (MQ=255)
                                                 gCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcac                                                                       <  7:1224/149‑1 (MQ=255)
                                                                 aTACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                                                >  5:108818/1‑140 (MQ=255)
                                                                 aTACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                                                <  6:108818/140‑1 (MQ=255)
                                                                                         gTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAcc                                                                           >  2:59721/1‑105 (MQ=255)
                                                                                         gTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAcc                                                                           <  1:59721/105‑1 (MQ=255)
                                                                                             tatTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAccc                           <  2:6523/149‑1 (MQ=255)
                                                                                                        gcAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                <  5:21464/117‑1 (MQ=255)
                                                                                                        gcAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                >  6:21464/1‑117 (MQ=255)
                                                                                                          aaCGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacggga                                                                                                         >  7:77594/1‑58 (MQ=255)
                                                                                                          aaCGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacggga                                                                                                         <  8:77594/58‑1 (MQ=255)
                                                                                                          aaCGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCgg              >  1:114701/1‑149 (MQ=255)
                                                                                                                                       cccgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt  >  7:59543/1‑133 (MQ=255)
                                                                                                                                       cccgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt  <  8:59543/133‑1 (MQ=255)
                                                                                                                                             |                                                                                                                              
GGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTT  >  NC_000913/3629355‑3629621

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