Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,703,2520AT27.3% 21.3 / 12.4 22intergenic (+18/+16)add/ydgJadenosine deaminase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (7/9);  new base T (3/3);  total (10/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 = 1.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCC  >  NC_000913/1703105‑1703397
                                                                                                                                                   |                                                                                                                                                 
taCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaaa                                                                                                                                                   <  1:8137/148‑1 (MQ=255)
     ttGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAcc                                                                                                                                             <  6:141505/149‑1 (MQ=255)
             ccAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                                                                                          >  1:131942/1‑144 (MQ=255)
             ccAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                                                                                          <  2:131942/144‑1 (MQ=255)
                     tgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCACCATCCAttt                                                                                                                                     <  1:73086/141‑1 (MQ=255)
                     tgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCACCATCCAttt                                                                                                                                     >  2:73086/1‑141 (MQ=255)
                                                              gTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCACCATCCAttttgtta                                                                                                                                >  1:70709/1‑101 (MQ=255)
                                                              gTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCACCATCCAttttgtta                                                                                                                                <  2:70709/105‑5 (MQ=255)
                                                               tCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGc                                                                                   >  1:10522/1‑149 (MQ=255)
                                                                           tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAcc                                                                                                                                             >  6:123519/1‑79 (MQ=255)
                                                                           tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAcc                                                                                                                                             <  5:123519/79‑1 (MQ=255)
                                                                             tCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCGCAATCAACTCCa                                                                      <  2:110285/148‑1 (MQ=255)
                                                                              ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACaaa                                                                                                                    >  6:11533/1‑101 (MQ=255)
                                                                              ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACaaa                                                                                                                    <  5:11533/101‑1 (MQ=255)
                                                                                tCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTAc                                                                  >  2:22310/1‑149 (MQ=255)
                                                                                         aGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGc                                                         <  2:10522/149‑1 (MQ=255)
                                                                                           gAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCt                                                       <  1:22310/149‑1 (MQ=255)
                                                                                                                                   aaaaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGcc                                                                                       <  2:174516/77‑1 (MQ=255)
                                                                                                                                   aaaaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGcc                                                                                       >  1:174516/1‑77 (MQ=255)
                                                                                                                                    aaaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCa              >  7:128672/1‑149 (MQ=255)
                                                                                                                                              tGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCAc          <  4:3218/143‑1 (MQ=255)
                                                                                                                                              tGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCAc          >  3:3218/1‑143 (MQ=255)
                                                                                                                                                cAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGcc  <  8:128672/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCC  >  NC_000913/1703105‑1703397

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: