Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,892,0290GC42.3% ‑5.2 / 12.5 27intergenic (+43/‑477)PP_RS08755/napAtRNA‑Asp/Periplasmic nitrate reductase
Reads supporting (aligned to +/- strand):  ref base G (7/8);  new base C (1/10);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 8.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATA‑TTCCACTTCCAGGGCCCACTGAACACCCGGAAGT‑AAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATC  >  KT2440/1891904‑1892160
                                                                                                                               |                                                                                                                                   
tataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATA‑TTCCACTTCCAGGGCCCACTGCACACCCGGAAGT‑AAAAAAGAAAGCGAcc                                                                                                                            >  7:20122/1‑135 (MQ=17)
 ataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGGCCGCGGCGGCCAA‑TTACACATGCAGGGCCCACAGCACACCCCGAAGT‑AAAAAATAAAaccaccg                                                                                                                           >  2:313244/1‑134 (MQ=11)
    gcgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGGCCGCTGCGCCCTA‑ATACACATGCAGGGCCCACAGCACCCCCGGAAGT‑AAAAAAGAAAGCGACCATAg                                                                                                                        >  6:117094/1‑135 (MQ=11)
           ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATA‑TTCCACTTCCAGGGCCCACTGCACACCCGGAAGT‑AAAAAAGAAAGCGACCTTAGGGTCGcg                                                                                                                 >  3:72647/1‑134 (MQ=21)
           ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCGGGGGCATA‑TTGCACATCCACGGCCCCCGGACCACCCCGAAGT‑AAAAAAGAAATCGACCGTAGGGTCggg                                                                                                                 >  2:85102/1‑133 (MQ=11)
           ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCGGCGGCATA‑TTCCACTTCCAGGGCCCCCAGACCACCCGGAAGT‑AAAAAAGAAAGCGACCGTAGGGTCggg                                                                                                                 >  5:165604/1‑133 (MQ=11)
             tAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGGCCGCTGCGGCCTA‑TTCCACGTCCAGGGCCCACTGAACACCCCGAACT‑AAAAAAGAAAGCGACCTTAGGGgcgcggg                                                                                                               >  2:252247/1‑132 (MQ=11)
             tAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGGTCGCGTCCCGTGCGCTGCGCCATA‑TTACACTTCCAGGGCCCACTGCCCCCCCCGTAGT‑AAAAAAGAAAGCGACCTTAGGGTCggggg                                                                                                               >  7:267860/1‑131 (MQ=11)
                       ggtgataaTCCCGGCCGGCCACGCCGGGGGGCGCGGGGTCGCTCCCCGGCCGCTGCGCAATA‑TTCCACTACCAGGGCCCACTGACCACCCGCACGT‑AAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCt                                                                                                     <  1:294126/129‑1 (MQ=255)
                        gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATA‑TTACACTTCCAGGGCCCCCTGCACACCCGGAAGT‑AAAAAACAAACCGACCTTAGGGGCGCGTGttttttttttt                                                                                                    >  1:224901/1‑132 (MQ=11)
                              aacccGGCCTGCCACGCGGGGGGGCGCGGGGGCGCGGCCCGCCCCCTCCGCCATA‑TTCCACTACCCGCGCCCACTCACCACCCGCACTA‑AAAAAACAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGcatcat                                                                                              <  7:339333/132‑1 (MQ=255)
                                                 ggggGCGCGGGTTCGAGTCCCGTCCGCTGCGCCATA‑TTCCACTTCCAGGGCCCACTGAACACCCGGAAGT‑AAAAAAGAAAGCGACCTTCGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGc                                                                           <  1:116074/135‑1 (MQ=255)
                                                   gggcgcgGGTTCGAGGCCCGTCCGCTGCGCCATA‑TTCCACTTCCAGGGCCCACTTAACACCCGTAAGA‑AAAAAAGAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGcaa                                                                         <  4:155389/132‑1 (MQ=18)
                                                           gTTCGAGTCCCGTCCGCTGCGCCATA‑TTCCACTTCCAGGGCCCACTGAACACCCGGAAGT‑AAAAAACAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGc                                                                 <  4:143188/135‑1 (MQ=255)
                                                               gcggcccGGCCGGTCCGCCAAT‑TTCCACTACCCGGCCCCACTGAACACCCGAAATA‑AAAAAACACCGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAg                                                              <  5:244558/130‑1 (MQ=255)
                                                                 gccccGGCCGCTGCGCAATA‑TCCCACTTCCCGGCCCCACTCAACACCCGCACTA‑AAAAACCAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTg                                                           <  1:192821/133‑1 (MQ=255)
                                                                  gcccgcccgCGGCG‑CATATTTCCACTCCCCGCCCCCACAAAAC‑CCCGTAAATAAAAAAACAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTgc                                                          <  1:322242/134‑1 (MQ=255)
                                                                                         ccctcccAGGCCCCACAGACCCCCCCGAACT‑AAAACACAACGCGACCTTAGGGTCGCGTTTTTTATTTCTGCATCATGCAATCGATTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGt                                    <  1:98467/130‑1 (MQ=255)
                                                                                         cACTTCCAGGGCCCACTCAACACCCGCAAGT‑AAAAAAGACAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGt                                    <  3:125512/135‑1 (MQ=255)
                                                                                         cACTCCCCGGCCCCACTCACCACCCGCAATA‑AAAAAACAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCTCGCAGCGt                                    <  1:313244/135‑1 (MQ=255)
                                                                                          aCTTCCAGGGCCCACTGAACCCCCGGAAGT‑AAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGt                                    <  4:72647/134‑1 (MQ=255)
                                                                                               cAGGGCCCACTGAACACCCGGAATT‑AAAAAAGAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCa                               <  2:224901/134‑1 (MQ=255)
                                                                                                   ccccACTCACCCCCCGCAAGA‑AAAAAACAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCa                          <  2:94776/134‑1 (MQ=255)
                                                                                                    cccACTGACCACCCGCACTT‑AAAAAACACCGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCaa                         <  1:209896/135‑1 (MQ=255)
                                                                                                                     aGT‑AAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCAcc        <  8:20122/135‑1 (MQ=255)
                                                                                                                      gT‑AAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTATACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCg       <  4:87508/135‑1 (MQ=255)
                                                                                                                       t‑AAAACACAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCgg      <  2:136871/135‑1 (MQ=255)
                                                                                                                            acagaacgAGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCGACTCTTCTTCCCGCAGCGTCAGCACCTCAACACCGTCCTCGATCACCGGCAtc  <  6:332488/133‑1 (MQ=255)
                                                                                                                               |                                                                                                                                   
TATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATA‑TTCCACTTCCAGGGCCCACTGAACACCCGGAAGT‑AAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATC  >  KT2440/1891904‑1892160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑