Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,715,3490CA66.7% 33.8 / 40.4 39intergenic (‑29/‑74)yfiE/eamBputative DNA‑binding transcriptional regulator/cysteine and O‑acetylserine exporter
Reads supporting (aligned to +/- strand):  ref base C (5/8);  new base A (12/14);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-02
Rejected as consensus: Frequency below/above cutoff threshold.

AATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGT  >  NC_000913/2715213‑2715443
                                                                                                                                        |                                                                                              
aaTATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaa                                                                                     >  3:123622/1‑148 (MQ=255)
         aGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAATACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAAcattcatt                                                                           <  6:93625/149‑1 (MQ=255)
                       gtgtACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTc                                                              <  4:40207/148‑1 (MQ=255)
                        tgtACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTcc                                                             <  4:10869/148‑1 (MQ=255)
                              gCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg                                                      <  4:123622/149‑1 (MQ=255)
                                      tttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGt                                               <  5:66646/148‑1 (MQ=255)
                                          cgagcttCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTc                                              <  4:50601/139‑1 (MQ=255)
                                          cgagcttCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTc                                              >  3:50601/7‑145 (MQ=255)
                                                   cAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACagag                                 >  8:136011/1‑149 (MQ=255)
                                                           aaCCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg                                                      >  8:59039/1‑120 (MQ=255)
                                                           aaCCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg                                                      <  7:59039/120‑1 (MQ=255)
                                                                tctTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGt                     >  2:23382/1‑148 (MQ=255)
                                                                      accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaaa                                                                   >  1:110441/1‑96 (MQ=255)
                                                                      accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaaa                                                                   <  2:110441/96‑1 (MQ=255)
                                                                      accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCg              >  7:29176/1‑149 (MQ=255)
                                                                      accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACagag                                 >  2:40081/1‑130 (MQ=255)
                                                                      accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACagag                                 <  1:40081/130‑1 (MQ=255)
                                                                        cacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTc                                                        <  8:55321/105‑1 (MQ=255)
                                                                        cacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTc                                                        >  7:55321/1‑105 (MQ=255)
                                                                         acGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGAcc           >  1:147385/1‑149 (MQ=255)
                                                                          cGGTTTTAAGCGTAATAAAACGGAGCAGATCCATAACCCCAAACCTAACGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGAccc          <  6:104685/149‑1 (MQ=255)
                                                                             ttttAATCGTAATAAAACGGCGCAGATCCATAAACCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGgaga                             <  2:10034/127‑1 (MQ=255)
                                                                             ttttAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATAGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGgaga                             >  1:10034/1‑127 (MQ=255)
                                                                               ttAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCtttt      <  2:147385/148‑1 (MQ=255)
                                                                               ttAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGaa                                                                                           >  4:52480/1‑63 (MQ=255)
                                                                               ttAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGaa                                                                                           <  3:52480/63‑1 (MQ=255)
                                                                                 aaGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTaa    <  8:29176/148‑1 (MQ=255)
                                                                                 aaGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg                                                      >  5:110817/1‑98 (MQ=255)
                                                                                 aaGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg                                                      <  6:110817/98‑1 (MQ=255)
                                                                                 aaGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTaa    <  7:136011/148‑1 (MQ=255)
                                                                                  aGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGt  <  1:23382/149‑1 (MQ=255)
                                                                                      taataaAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGt                     >  2:125007/1‑126 (MQ=255)
                                                                                      taataaAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGt                     <  1:125007/126‑1 (MQ=255)
                                                                                        ataaAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCgtg                                                          <  5:80053/87‑1 (MQ=255)
                                                                                        ataaAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCgtg                                                          >  6:80053/1‑87 (MQ=255)
                                                                                               ggCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg                                                      >  4:103781/1‑84 (MQ=255)
                                                                                               ggCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg                                                      <  3:103781/84‑1 (MQ=255)
                                                                                                                         tCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGt  >  7:58022/1‑110 (MQ=255)
                                                                                                                         tCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGt  <  8:58022/110‑1 (MQ=255)
                                                                                                                                        |                                                                                              
AATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGT  >  NC_000913/2715213‑2715443

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: