Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,505,5000CT66.7% 7.0 / 12.2 18L72L (CTG→TTG) nanY2,7‑anhydro‑N‑acetylneuraminate hydratase
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base T (7/5);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCCAAATTATCTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTTTGTGAAAAACCAATTGCATTAAGTTATGAAGATTGTGTGGATATGGTCAAAGCGTGTAAAGAAGCTGGTGTGACCTTTATGGCCGGGCATATTATGAATTTTTTCAATGGGGTTCAATATGCACGGAAGTTAATTAAAGAAGGTGTTATCGGCGAAATATTATCATGTCATACTA  >  NC_000913/4505282‑4505725
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
gCTTTATGAGTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCCAATTTATCTGCATAAAGAACCATTAATTAAAGCAGCAAAg                                                                                                                                                                                                   <  1:99759/251‑1 (MQ=255)
                             gttggTTACTTTGGAGCTGAATTAGCTCGTTTTATGAAGATGCATGATAATGCAAAAAATACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGGCAAGCTTGGATTCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCCAAATTATTTGCATAAAGAACCATTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTTTGTGaaaaan                                                                                                                                                                       <  1:58649/250‑2 (MQ=255)
                                          gcgcTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCCAAATTATTTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGGTGTTTGTTAAAAACCAATTGCATTAAGt                                                                                                                                                         <  1:270036/251‑1 (MQ=255)
                                               gAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGCTCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCCAAATTATCTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTTTGTGAAAAACCAATTGCATTAAGTTATg                                                                                                                                                     >  1:102288/1‑250 (MQ=255)
                                                          ttttATGAATATGCATGATAATGCAANAATTACNNGTNTATACGATCCTGAAAATGGNGANNNTNNNGNACGTGNATNNCAGTGTATCNATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCNTAGCCCCCCNNAATTATCTGCATAAAGAACCANTAANNANNGNAGCnaan                                                                                                                                                                                                   >  2:73187/1‑192 (MQ=255)
                                                                                 cAAAAATTACATGTGTATACGATACTGAAAATAGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCCAAATTATTTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGc                                                                                                                                                                                         >  1:116581/1‑180 (MQ=255)
                                                                                                        ccTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCCAAATTATCTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTTTGTGAAAAACCAATTGCATTAAGTTATGAAGATTGTGTGGATATGGTCAAAGCGTGTAAAGAAGCTGGTGTGACCTTTATGGCCg                                                                                            <  1:73187/250‑1 (MQ=255)
                                                                                                                             gCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTAGATTGCGTGATCGTAGCCACCCCCAATTATTTGCATAAAGAACCAGTAATTAAagccgaaaag                                                                                                                                                                                                   >  1:288175/1‑119 (MQ=255)
                                                                                                                                        ncAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATAGTAGCCACCCCAAATTATTTGc                                                                                                                                                                                                                                >  1:37480/2‑86 (MQ=255)
                                                                                                                                         cAGTGTATCAATATGTCAATCTTGGATTCTTTAGTCTCAAGTAAATTAGTCGATTTCGTGATCGTAGCCACCCCAAATTATCTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTGTGTGAAAAACCAATTGCATTAAGTTATg                                                                                                                                                     <  1:109024/160‑1 (MQ=255)
                                                                                                                                                                     cTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCCAAATTATTTGCGTAAAGAACCAGTAATTAAAGCAGAAAAGACTAAGACGAATGTTTTTTGTGCAAAAACAATTGCATTAAGTTATGAAGATTGTGTGGATATGGTCAAAGCGTGTAAAGAAGCTgg                                                                                                             >  1:170260/1‑172 (MQ=255)
                                                                                                                                                                                     aTTAGTCGATTGCGTGATCGTAGCCTCCCCAAATTATTGGCATAAAGAACCTGTAATTAAAGCAGCAAAGAATTAGAAGTATGGTTGTTGTGAAAGAGCAATTGCATTAATTTGTGATGATTGGGTGGATATGGGGAAATCTTGTGAAGAGTCTTCTGTGACCTTTATGTTCGGGCATATTATGAATTTTTTCAATGGGGTTCAatat                                                         <  1:165618/208‑1 (MQ=255)
                                                                                                                                                                                                gCGTGATCGTAGCCACCCCAAATTATCTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTTTGTGAAAAACCAATTGCATTAAGTTATGAAGATTGTGTGGATATGGTCAAAGCGTGTAAAGAAGCTGGTGTGACCTTGATGGCCGGGCAGATTATGAATTTTTTCAATGGGTTTCAATAGGCACGGAAGTTAATTAAAGAAGGTGTTATCggn                        <  1:36630/230‑2 (MQ=255)
                                                                                                                                                                                                  gttcgcGTAGCCACACCAATTTATTTGCATAAAGCACCAGTCATTAAAGCAGCCAAGAATCAGAAGAATGTTTTTTGTGGAAAAACAATGGCATTAAGTTATGAGGATTGTGTGGATATGGTCAAAGCGTGTACAGATGCTGGTGTGACATTTATGGCCGGGCATATTATGAATTTTTTCAATGGGGATCAATATGCACGGACGTTCATTAACGAAGGTGTTATAGGCGCAATATTATCATGTCATACTa  >  1:164903/6‑250 (MQ=255)
                                                                                                                                                                                                   tGATCGTAGCCACCCCAAATTATTTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTTTGTGAAAAACCAATTGCATTAAGTTATGAAGATTGTGTGGATATGGTCAAAGCGTGTAAAGAAGCTGGTGTGACCTTTATGGCCGGGCATATTATGAATTTTTTCAATGGGGTTCAATATGCACGGAAGTTAATTAAAGAAGGTGTTATCGGCg                       <  1:13297/228‑1 (MQ=255)
                                                                                                                                                                                                   tGATCGTAGCCACCCCAAATTATTTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTTTGTGAAAAACCAATNNCATTAANnt                                                                                                                                                        >  2:13297/1‑99 (MQ=255)
                                                                                                                                                                                                              aCCCCAAAGTATTTGCATAAAGAGCCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTATTGTTAAAAACCAATTGCAGTAAGTTAGGAAGATTGTGTGGATAGGATCAAAGCGTGTAAAGGAGCTTGTGTGACCTTGAGGGCGGGGCATATTAGGAGTTTTTTCAATGGGGTTCAATATGCACGGAAGTTAATTAAAGAAGGTGtt                              <  1:315879/210‑1 (MQ=255)
                                                                                                                                                                                                               ncccAAATTATTTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTTTGTGAAAAACCAATTGCATTAAGTTATGACGATTGTGTGGATATGGTCAAAGCGTGTAAAGAAGCTGGTGTGACCTTTATGGCCGGGCATATTATGAATTTTTTCAATGGGGTTCa                                                             >  1:62444/2‑178 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
GCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCCAAATTATCTGCATAAAGAACCAGTAATTAAAGCAGCAAAGAATAAGAAGCATGTTTTTTGTGAAAAACCAATTGCATTAAGTTATGAAGATTGTGTGGATATGGTCAAAGCGTGTAAAGAAGCTGGTGTGACCTTTATGGCCGGGCATATTATGAATTTTTTCAATGGGGTTCAATATGCACGGAAGTTAATTAAAGAAGGTGTTATCGGCGAAATATTATCATGTCATACTA  >  NC_000913/4505282‑4505725

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