Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,732,7010TC53.3% ‑2.4 / 34.4 30A1205A (GCT→GCC) lhrATP‑dependent helicase Lhr
Reads supporting (aligned to +/- strand):  ref base T (9/5);  new base C (9/7);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCTCAACATCGCCAGGCGTTTATGAGGGCGTAGATGGCGTGATGCGGGTGATCGAACAGCTTGCCGGAGTCGGTTTACCCGCCTCACTCTGGGAAAGCCAGATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTT  >  NC_000913/1732463‑1732945
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
cctaagaaccgccgggagtttaggtgggCGTAGGTGGGGTGAGGCGGGTGGGCGACCAGCTTGCCGGGGTCGTTTTCCCCGCCTCACTCTGGGAAAGCCAGATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTc                                                                                                                                                                                                                                          <  2:87044/226‑1 (MQ=255)
                             gTAGATGGCGTGATGCGGGCGATCGAACAGCTTGCCGGAGTCGGTTTACCCGCCTCACTCTGGGAAAGCCAGATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCg                                                                                                                                                                                                             >  2:132330/1‑251 (MQ=255)
                              tAGATGGCGTGATGCGGGTGATCGAACAGCTTGCCGGAGTCGGTTTACCCGCCTCACTCTGGGAAAGCCAGATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGt                                                                                                                                                                                                            >  2:314226/1‑251 (MQ=255)
                                                              ggatgaggcgggttACCCGACTCACTAAGGGACATCCAGATCCTGCATACCCGCGTACGCGTCTATTTGTCAGACATGCTATATGAATTACTGGCACACGGGGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGCAGACGACGGCCCGGGGGCAATGCCTCTCCAGGAATCAGCCGCAGACACGTTCACGCCGGCCGAAGCGGATCAGGCGAATCGATCGGCGCCGCCACAAGAGAGCGTCGCTggtgt                                                                                                                                                                            >  2:117920/13‑247 (MQ=255)
                                                                gggaggctgTTTCCCCGCCTCACTCTGGGAAAGCCAGATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTgg                                                                                                                                                                          <  1:166289/243‑1 (MQ=255)
                                                                     ttggtttccccgccgccatctGGGAAAGCCAGATCCTCCCTGCCCGCTTACGCTACTATTCGTCAGAAGTGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGTTGGCACTGCAACTACAGAAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGa                                                                                                                                                                      <  1:172844/232‑1 (MQ=255)
                                                                        gTTTACCCGCCTCACTCTGGGAAAGCCAGATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACggg                                                                                                                                                                  >  1:230940/1‑250 (MQ=255)
                                                                            ccccGCCTCACTCTGGGAAAGCCAGATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACggagga                                                                                                                                                               <  1:314226/249‑1 (MQ=255)
                                                                             cccGCCTCCCTCTGGGAAAGCCAGATCTTGCCTGCCCGCGTCCGCGACTCTTAGTCAGAAATGCTCGATGAATTACTTGCACCCGGTGCGGTCATCTGGCCGGGGAAACAAAAACTGGGTGAAGATGACGGCGTGGTGACACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAACGGATCAGGCACATCATTCGGCGCTGCAACATGCGATAGTCGCTGTTCTTGCtaacaggaaagc                                                                                                                                                             >  2:305219/1‑240 (MQ=255)
                                                                                                  cAGATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATc                                                                                                                                          >  2:89613/1‑249 (MQ=255)
                                                                                                   agagcctgcctaCCCGCGCACGCGACAATTCGTCAGAAATGATCGATGAATTACTGGCAAACGGTGAGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTACAGAATAGTTCACTCCCGACGAAGCGGATCAGGCGAATCGTTCGGCGCGGCACCAAGCGATAGACGCTGTTCGGGCTGACGGAGGGGCCTGGGTTGccacaccaatcagc                                                                                                                                       >  2:194315/5‑238 (MQ=255)
                                                                                                    gATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGcc                                                                                                                                      <  1:132330/251‑1 (MQ=255)
                                                                                                        ctgcctgcCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGACGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACg                                                                                                                                  >  1:237472/1‑251 (MQ=255)
                                                                                                            ctgccCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGAt                                                                                                                               >  1:77755/1‑250 (MQ=255)
                                                                                                              gccCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGTGCAAAAAAAACTGGGGGAAGATGACGGCCTGGGGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAAGCGTTCGGCGCTGCAACAAGCGATAGTCGCGGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAGCAAATCAGCCAACGGATa                                                                                                                              <  2:206291/249‑1 (MQ=255)
                                                                                                                                   tCAGAATTGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGt                                                                                                       >  1:85431/1‑251 (MQ=255)
                                                                                                                                                      tACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAACAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATACTGCAGAATCGTTCACTCCCGCCGAAGCGGCTCAGGCGAATCGTTCGGCGCTGCAACACGCGATAGTCGCTGTTCTGGCTGACGGAGGATCCTTGTTTGCACAACCAATCAGCCAACGGATACGCGACAAAATCGGCGCATCGGTAGATCTCTCGGCCCtcca                                                                                     >  2:92739/1‑247 (MQ=255)
                                                                                                                                                                          aTCTGGTAGGGTCAAAAAAAACTGGTTGAAGATGACGGCCTGGTGGTACTGCATCGACATGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTgg                                                                 <  2:77755/250‑1 (MQ=255)
                                                                                                                                                                               gTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAAAAAATCAGCCAACGGATACGCGCCAAAATCGTCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTggg                                                           >  2:61722/1‑250 (MQ=255)
                                                                                                                                                                                        aaaaaaaCTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAAGCGtca                                                  <  1:89613/251‑1 (MQ=255)
                                                                                                                                                                                          aaaaaCTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGtcatc                                                >  1:121878/1‑251 (MQ=255)
                                                                                                                                                                                                tGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGc                                          >  1:155892/1‑251 (MQ=255)
                                                                                                                                                                                                          tGACGGCCTGGTGGCACTGCATCTACAGGAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTggg                                >  2:32967/1‑251 (MQ=255)
                                                                                                                                                                                                            aCGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGc                               >  2:131642/1‑250 (MQ=255)
                                                                                                                                                                                                                    ggggCACGGCTTCTACAGGAATATGCTGCGGAATCGTCCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTAcg                      <  2:91059/249‑1 (MQ=255)
                                                                                                                                                                                                                              aTCTACAGGACTATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCCGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGCGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCAcc            >  2:189816/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggAATATGCCGCAGAATCGTCCCCTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCgcagca     <  2:6938/251‑1 (MQ=255)
                                                                                                                                                                                                                                        atatGCCGCAGAATCGTTCCCTCCCGACGAAGCGGATCAGGCGAAACGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCTTTACGCGCCCTCACCCGCAGCAGtt  <  2:230940/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tatGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTa                                                                           >  1:182116/1‑177 (MQ=255)
                                                                                                                                                                                                                                         tatGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAAGCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTa                                                                           <  2:182116/177‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
CCTCAACATCGCCAGGCGTTTATGAGGGCGTAGATGGCGTGATGCGGGTGATCGAACAGCTTGCCGGAGTCGGTTTACCCGCCTCACTCTGGGAAAGCCAGATCCTGCCTGCCCGCGTACGCGACTATTCGTCAGAAATGCTCGATGAATTACTGGCAACCGGTGCGGTTATCTGGTCGGGGCAAAAAAAACTGGGTGAAGATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTT  >  NC_000913/1732463‑1732945

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