Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,260,9040AG28.0% 32.7 / 14.1 25intergenic (‑101/+24)dauA/prsC4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (10/8);  new base G (4/3);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAA  >  NC_000913/1260760‑1261052
                                                                                                                                                |                                                                                                                                                    
agagCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCgg                                                                                                                                                  >  1:137959/1‑149 (MQ=255)
             gCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCgg                                                                                                                                     <  2:137959/149‑1 (MQ=255)
                                                                             aaCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAt                                                                                                                              <  4:41201/92‑1 (MQ=255)
                                                                             aaCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAt                                                                                                                              >  3:41201/1‑92 (MQ=255)
                                                                             aaCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacga                                                                                 <  8:20387/137‑1 (MQ=255)
                                                                             aaCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacga                                                                                 >  7:20387/1‑137 (MQ=255)
                                                                                                   gCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCa                                                                                                            <  2:78634/88‑1 (MQ=255)
                                                                                                   gCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCa                                                                                                            >  1:78634/1‑88 (MQ=255)
                                                                                                     aTAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAg                                                                <  2:17271/130‑1 (MQ=255)
                                                                                                     aTAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAg                                                                >  1:17271/1‑130 (MQ=255)
                                                                                                     aTAGGTCATACAAATGGATATTACAGACAAAAAAACACGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATg                                                                                                               >  5:96287/1‑83 (MQ=255)
                                                                                                     aTAGGTAATACAAATGGATATTACAGACAAAAAATCCCGCGGGAGCGGGTCTTTGAGGCGGGTTCGATTAGTGTTCGTACATg                                                                                                               <  6:96287/83‑1 (MQ=255)
                                                                                                              aCAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGt                                                 >  3:87053/1‑136 (MQ=255)
                                                                                                              aCAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGt                                                 <  4:87053/136‑1 (MQ=255)
                                                                                                                      atatTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAg                            >  3:135337/1‑149 (MQ=255)
                                                                                                                       tatTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGt                           >  3:45924/1‑149 (MQ=255)
                                                                                                                                 tcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTc                                                <  3:103111/116‑1 (MQ=255)
                                                                                                                                 tcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTc                                                >  4:103111/3‑118 (MQ=255)
                                                                                                                                 aaaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt                 <  4:135337/149‑1 (MQ=255)
                                                                                                                                  caaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttga                    <  2:117707/144‑1 (MQ=255)
                                                                                                                                  caaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttga                    >  1:117707/2‑145 (MQ=255)
                                                                                                                                   aaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCa                                               >  8:55380/1‑117 (MQ=255)
                                                                                                                                   aaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCa                                               <  7:55380/117‑1 (MQ=255)
                                                                                                                                   aaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCa               >  3:96826/1‑149 (MQ=255)
                                                                                                                                         ccgccgCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCt          >  7:43564/1‑148 (MQ=255)
                                                                                                                                                aGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGaa  <  8:43564/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAA  >  NC_000913/1260760‑1261052

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: