Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4960TA27.0% 72.9 / 13.9 37intergenic (+301/‑155)dapB/carA4‑hydroxy‑tetrahydrodipicolinate reductase/carbamoyl phosphate synthetase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (23/4);  new base A (0/10);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 2.87e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGA  >  NC_000913/29367‑29631
                                                                                                                                 |                                                                                                                                       
ttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                               >  1:103336/1‑140 (MQ=255)
                   tCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTa                                                                                                            >  1:1382505/1‑140 (MQ=255)
                   tCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTa                                                                                                            >  1:1297845/1‑140 (MQ=255)
                        ttttttTCATGCATTCTTCTTGCTGCGCAAGCGTTTCCAAAAACGGGTTAGATGTTTTTTTTGTCGTTTAACGCCTGAAAAACTTGCATGACCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                       <  2:509824/140‑1 (MQ=255)
                        ttttttTCATGAATTTTTCTTCCTGCGAAAGCTTTCTCCAAAACGGTTTAGAAAATTTTTTTGTCGTGAAGTGCGTTAAAAACATGCATGACCCAAAAAATAATAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                       <  2:369481/140‑1 (MQ=255)
                        ttttttTCATGAAATTTTCTTGCTGCGCAAGCTTTTCCCAAAACAGGTTAGATGATTTTTTTGTCGCTTAATCCTTAAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                       <  2:653307/140‑1 (MQ=255)
                           ttttCATCCATTCTTCTTGCGGCGAAAGCTTTTACAAAAACAGGTTAGATGATTTTTTTGTCGCTAAATCCTTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                    <  2:109608/139‑1 (MQ=255)
                                 tGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                              >  1:53679/1‑140 (MQ=255)
                                         tcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                      >  2:958733/1‑140 (MQ=255)
                                         tcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                      >  1:592864/1‑140 (MQ=255)
                                           ttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                    >  1:270067/1‑140 (MQ=255)
                                               tgcGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                <  1:982165/140‑1 (MQ=255)
                                                gcgcAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                               >  2:1388504/1‑140 (MQ=255)
                                                   cAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                            >  2:700286/1‑140 (MQ=255)
                                                                    gggTTAGATGTTTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                           <  2:1085794/139‑1 (MQ=255)
                                                                      gTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                         >  1:1273560/1‑140 (MQ=255)
                                                                      gTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                         >  1:1085794/1‑140 (MQ=255)
                                                                             tatttttttGTCGTTTAATGCCTGAAAAACATGAATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGc                                                  <  1:880672/136‑1 (MQ=255)
                                                                                cTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                               >  1:1407008/1‑140 (MQ=255)
                                                                                  ttttGTCGCTTAATGCCTGAAAAACATGCATGAGCAAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa                                             <  1:1048610/140‑1 (MQ=255)
                                                                                   tttGTCGTTTATTCCCTAAAAAACATGCATGCCCCAAAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAt                                            <  2:1100446/140‑1 (MQ=255)
                                                                                     tGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                          >  1:1438514/1‑140 (MQ=255)
                                                                                             aaTGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                  <  2:121796/140‑1 (MQ=255)
                                                                                                ccctaaaaaaCATGCAGGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                               <  2:270067/134‑1 (MQ=255)
                                                                                                         aCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTaa                      >  2:638251/1‑140 (MQ=255)
                                                                                                             gcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                  >  2:960875/1‑140 (MQ=255)
                                                                                                             gcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                  >  2:919751/1‑140 (MQ=255)
                                                                                                             gcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                  >  2:272019/1‑140 (MQ=255)
                                                                                                             gcaggAGCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                  <  2:592864/140‑1 (MQ=255)
                                                                                                               atgAGCCACAAAGTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                >  2:654183/1‑140 (MQ=255)
                                                                                                                 gacccAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa              <  2:652768/137‑1 (MQ=255)
                                                                                                                  aGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt             >  1:652768/1‑140 (MQ=255)
                                                                                                                   gCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa            >  1:756342/1‑140 (MQ=255)
                                                                                                                    ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa           >  1:1102175/1‑140 (MQ=255)
                                                                                                                    ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa           >  1:1098042/1‑140 (MQ=255)
                                                                                                                    caaaaaaaTAAAAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa           <  2:54597/136‑1 (MQ=255)
                                                                                                                        aaaaaaaaaaaaaaaaCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg       <  2:1102175/140‑1 (MQ=255)
                                                                                                                        aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg       <  1:272019/140‑1 (MQ=255)
                                                                                                                          aataatAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgag     <  1:1103319/140‑1 (MQ=255)
                                                                                                                            taataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtt   >  1:111419/1‑139 (MQ=255)
                                                                                                                             aataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga  <  1:700286/140‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                       
TTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGA  >  NC_000913/29367‑29631

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑