Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,360,6450TA26.1% 31.6 / 12.6 23intergenic (+29/‑531)gltD/gltFglutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (8/9);  new base A (3/3);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTT  >  NC_000913/3360499‑3360791
                                                                                                                                                  |                                                                                                                                                  
cAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTgtgt                                                                                                                                                  >  2:47076/1‑149 (MQ=255)
    aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTgg                                                                                                                                              <  8:68647/149‑1 (MQ=255)
    aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTAAGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTgg                                                                                                                                              >  2:77933/1‑149 (MQ=255)
            aaTCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCa                                                                                                                                       >  7:34743/1‑148 (MQ=255)
                  tGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTACGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCGAAttcatt                                                                                                                                    <  1:77933/145‑6 (MQ=255)
                      ggcggcGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGAc                                                                                                                             >  6:115983/1‑148 (MQ=255)
                                                    gtggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCgg                                                                                               <  8:34743/148‑1 (MQ=255)
                                                          aCCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAg                                                                                         <  1:77010/148‑1 (MQ=255)
                                                            cGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTa                                                                                                                                  >  4:106849/1‑105 (MQ=255)
                                                                    ccGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCGGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGGTTTCTGGCAAATGCGGCATGAGCGTTTGACCtt                                                                              >  6:31524/1‑149 (MQ=255)
                                                                             gTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataa                                                                      <  5:115983/148‑1 (MQ=255)
                                                                                aGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAgg                                                                  <  1:47076/149‑1 (MQ=255)
                                                                                     gcAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATaa                                                             >  2:134660/1‑149 (MQ=255)
                                                                                                         gctggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGCCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCttt                                         >  5:81769/1‑149 (MQ=255)
                                                                                                         gctggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCttt                                         <  2:105633/149‑1 (MQ=255)
                                                                                                         gctggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCttt                                         >  1:105633/1‑149 (MQ=255)
                                                                                                                           gTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtgatga                       <  1:134660/149‑1 (MQ=255)
                                                                                                                              aCAATGAATTCGCCAGTCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGt                     <  6:81769/148‑1 (MQ=255)
                                                                                                                               cAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACTAAGTTCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtgat                         <  7:134961/143‑1 (MQ=255)
                                                                                                                               cAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtgat                         >  8:134961/1‑143 (MQ=255)
                                                                                                                                      ttCGCCAGACTCAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGc                                                                                              <  3:115331/67‑1 (MQ=255)
                                                                                                                                         gCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATtgtg                                             >  2:19779/1‑113 (MQ=255)
                                                                                                                                         gCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATtgtg                                             <  1:19779/113‑1 (MQ=255)
                                                                                                                                         gCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATc           >  2:156436/1‑147 (MQ=255)
                                                                                                                                                tgtgTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAAttt  >  8:34281/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTT  >  NC_000913/3360499‑3360791

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: