Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,232,4680AG54.5% ‑4.2 / 20.8 22intergenic (+14/+32)umuC/dsbBDNA polymerase V catalytic protein/protein thiol:quinone oxidoreductase DsbB
Reads supporting (aligned to +/- strand):  ref base A (5/5);  new base G (12/0);  total (17/5)
Fisher's exact test for biased strand distribution p-value = 9.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGC  >  NC_000913/1232405‑1232529
                                                               |                                                             
tttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGaa                                                              >  1:1936104/1‑65 (MQ=255)
 ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGgaa                                                             >  1:1713666/1‑65 (MQ=255)
    accacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGaaaaaa                                                          <  1:1441746/65‑1 (MQ=255)
     ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGaaaaaaa                                                         >  1:1215622/1‑64 (MQ=255)
          ttATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTc                                                    <  1:546499/65‑1 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAGAAAGCccgcccgc                                                >  1:2542295/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAGCGCTCCCGc                                                >  1:2315426/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAGCGCGCCCGc                                                >  1:981221/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAGCGCCCCCGc                                                >  1:2507495/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAGCACTccccc                                                >  1:1570140/1‑63 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCcccccc                                                >  1:646889/1‑63 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCTccccc                                                >  1:2087063/1‑63 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCTCCCGc                                                >  1:1525481/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCTCCCGc                                                >  1:1774880/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCTCCCGc                                                >  1:1099643/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCGccccc                                                >  1:2409381/1‑63 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCCACCGc                                                >  1:2769184/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGc                                                >  1:2326394/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGAGCCCCCGc                                                >  1:2336292/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAACGCcccccc                                                >  1:683129/1‑63 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAACGCTCCCGc                                                >  1:1663057/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAACGCTCCCGc                                                >  1:152958/1‑65 (MQ=255)
              aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAACGCTCCCGc                                                >  1:1287289/1‑65 (MQ=255)
                   gcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGTTGCCGCTGGaa                                           >  1:1220114/1‑64 (MQ=255)
                        gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCCCCCGCAGGAGGGCCg                                      >  1:1192516/1‑65 (MQ=255)
                            tGATTTACTGCGGGTCAAACAAATATAGCGGCAGAAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                  <  1:647497/65‑1 (MQ=255)
                             gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATgg                                 <  1:2524130/65‑1 (MQ=255)
                               tttACTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCTCCCGCAGGAGCGCCGAATGGac                               >  1:1472723/1‑64 (MQ=255)
                               tttACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                               <  1:1452448/65‑1 (MQ=255)
                                  aCTGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCTCCCGCAGGAGCGCCGAATGGATTAg                            >  1:579331/1‑65 (MQ=255)
                                    tGAGGGTCAAATAAATATAGCGGCAGGGAAAAAACGCCCCCGCACGAGCGCCCAATGGATTAGCg                          >  1:1703879/1‑65 (MQ=255)
                                         gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCCAGAGCGCCGAAAGGCTTAACGAccca                     >  1:94762/1‑63 (MQ=255)
                                         gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGa                     <  1:551963/65‑1 (MQ=255)
                                         gTCAAATAAATATAGCGGCAGGAAAAAAGCGCCCCCGCAGGAGCGCCGAATGGATTAGCGACCGa                     >  1:1228069/1‑65 (MQ=255)
                                              ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGa                <  1:1232401/65‑1 (MQ=255)
                                              ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGa                <  1:96787/65‑1 (MQ=255)
                                              ataaataTAGCGGCAGGAAAAAAACGCTCCCGCACGAGCGCCGACAGGATTAGCGACCGAACAGa                >  1:2073364/1‑65 (MQ=255)
                                              ataaataGAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGa                <  1:328349/65‑1 (MQ=255)
                                               taaataTAGCGGCAGGGAAAAAGCGCTCCCGCCGGAGCGCCGAATGGATTAGCGACCCAACacaa               >  1:480979/1‑62 (MQ=255)
                                               taaataTAGCGGCAGGGAAAAAACGCTGCCGCCGGAGCGCCGAATGGATTAGCGACCGAACAGAt               >  1:424191/1‑65 (MQ=255)
                                               taaataTAGCGGCAGGAAAAAAGCGCTCCCGCCGGGGCGCCGAATGGATTAACGAACGAACAGAt               >  1:1919928/1‑65 (MQ=255)
                                                            aGGAAAATAGCGCTCCCGCAGGAGCGGCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc  <  1:1331220/65‑1 (MQ=255)
                                                               |                                                             
TTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGC  >  NC_000913/1232405‑1232529

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