Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,439,602 (T)7→6 intergenic (‑93/+14) yhdN ← / ← rplQ conserved protein YhdN/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5960T.84.2% 20.2 / ‑1.8 19intergenic (‑87/+20)yhdN/rplQconserved protein YhdN/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base . (1/15);  total (3/16)
Fisher's exact test for biased strand distribution p-value = 5.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

CACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAAACCC‑‑GCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCG  >  NC_000913/3439377‑3439851
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                          
cACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAAACCC‑‑GCCGGGGCGGG‑‑‑TTTTACGTTGCTTCAGATTACTCTGCagcagc                                                                                                                                                                                                                                         <  1:1447659/251‑1 (MQ=255)
       ttCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAACCC‑‑GCCGGGGCGGGTTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                                                                                                                                                   <  1:552858/298‑1 (MQ=255)
                            atcatcCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGTATAAAA‑‑‑‑‑C‑‑GCGGGGCGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTGTCGCCTGCa                                                                                                                                                                         >  1:1447665/1‑278 (MQ=255)
                                                                gTTATCAAAACTCACCTTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAAACCACATAAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAACCC‑‑GCCGGGGCGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAACCAGCTCg                                                                                                                                                                                                        <  1:552866/224‑1 (MQ=255)
                                                                                             ttCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAAACCC‑‑GCCGGGCGGG‑‑‑‑TTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGCGCCTGCACGGAAGCCACACTTCAg                                                                                                                                                       <  1:1447670/235‑1 (MQ=255)
                                                                                               cTGCTATATGGCGCTCTGCCCACTGGTCAAGTAAACCACATAAAATACTCCAGGGGATGAAACAAAAAAGAGTACAGCTATTGTAGATAAAGT‑‑GGGGATACTACGGGTATAAAAAACCC‑‑GCCGGGGCGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAACGATCAAACCAGCTCGATGTAAAGCCATCGGCGCGTTGTCGCCTGCACGGAAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGCCCAGTTCGTTAAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCaa              <  1:552871/378‑1 (MQ=255)
                                                                                                  cTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAAACCC‑‑GCCGGGGCGG‑‑‑‑TTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAAcg                                                                                                                 <  1:1447671/266‑1 (MQ=255)
                                                                                                         gcgcTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAACCC‑‑GCCGGGGCGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAACACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAgac                                             <  1:552872/328‑1 (MQ=255)
                                                                                                              ctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAACCC‑‑GCCGGGGCGGG‑‑TTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAg                                                                <  1:552873/305‑1 (MQ=255)
                                                                                                                         tCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAACCC‑‑GCCGGGGCGGG‑‑TTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTGCCACGATCTCGTTATCACGATACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCtt                   <  1:552874/336‑1 (MQ=255)
                                                                                                                                                      gATGAACAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAAACCC‑‑GCCGGGGCGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTGTCGCCTGCACGGAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGAACGAtttagca                                <  1:552875/296‑6 (MQ=255)
                                                                                                                                                                            cTATTGTAGATAAGTGAGGGGATACGTACGGGTAATCAAAAAACCCCGGCCGGGGCGGG‑‑‑TTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGc                                  <  1:1447675/280‑1 (MQ=255)
                                                                                                                                                                                 gTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAAACCC‑‑GCCGGGGCGGGtttttta                                                                                                                                                                                                                                                                      >  1:1447678/1‑54 (MQ=255)
                                                                                                                                                                                        aaGT‑‑GGGGATACTACGGGTAATAAAAAACCC‑‑GGCCGGGCGGG‑‑‑TTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                                                                                                                                                       <  1:1447677/152‑1 (MQ=255)
                                                                                                                                                                                                          ggtaatcaaaaaacCCG‑‑GCCGGGGCGGG‑‑‑TTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGAtct                                                                             <  1:1447679/202‑1 (MQ=255)
                                                                                                                                                                                                                 taaaaaaCCC‑‑GCCGGGGCGGG‑‑‑TTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAAACGCGGGCCCAGTTCGTTaa                                                                                               <  1:552878/183‑1 (MQ=255)
                                                                                                                                                                                                                   aaaaaCCC‑‑GCCGGGGCGGGTTTTTTTATCGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGAcgacga                                       >  1:1447682/1‑239 (MQ=255)
                                                                                                                                                                                                                      aaCCC‑‑GCCGGGGCGGG‑‑TTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAAcg                                                                                                                 <  1:552879/163‑1 (MQ=255)
                                                                                                                                                                                                                        ccc‑‑GCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGTCTCGTTATCACGATGTACGGGCGAATGCCGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCg  >  1:1447684/1‑269 (MQ=255)
                                                                                                                                                                                                                         cccggccgGGGCGGG‑‑‑TTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCg                                                                           <  1:1447681/194‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                          
CACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGT‑‑GGGGATACTACGGGTATAAAAAAACCC‑‑GCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCG  >  NC_000913/3439377‑3439851

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