Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,741,285 2 bp→CT 100% intergenic (+476/+275) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,741,2850AC87.5% 38.9 / ‑0.8 16intergenic (+476/+276)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (7/7);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0029474,741,2860CT87.5% 37.3 / ‑1.2 16intergenic (+477/+275)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (7/7);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGAACCTCCGCCAACTGCATTCGTGAACCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAC  >  NC_002947/4741172‑4741404
                                                                                                                   ||                                                                                                                        
ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGTTATAGCTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCTGGacac                                                                                      <  2:252091/150‑1 (MQ=255)
                        caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGTTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                              <  2:250507/150‑1 (MQ=255)
                                     ggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                 <  1:612227/150‑1 (MQ=255)
                                       cgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                               <  2:633886/150‑1 (MQ=255)
                                                     gCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                 >  2:712435/1‑150 (MQ=255)
                                                        tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTAAATGCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCATCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                              >  2:154515/1‑150 (MQ=255)
                                                             aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                  >  1:749210/1‑125 (MQ=255)
                                                             aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                  <  2:749210/125‑1 (MQ=255)
                                                                        ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt              <  1:629360/150‑1 (MQ=255)
                                                                         gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc                                                                 <  1:277588/98‑1 (MQ=255)
                                                                         gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc                                                                 >  2:277588/1‑98 (MQ=255)
                                                                         gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg             >  1:668860/1‑150 (MQ=255)
                                                                             aGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTAACACCACGATAGCCCTTGTACGGGGAcc         >  2:418852/1‑150 (MQ=255)
                                                                               ctctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGCCTTGCATTCGTGAATCCAGGCACAGGCCGCGGACACGTCCTCCGGCAGCGCCTTGACCTTGGCCAGCGGCGCGTTGATCACGGTCAAACCGCGATAGGCCTTGTGCTTGCACCCg       >  2:331201/1‑150 (MQ=255)
                                                                                tctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCgg      >  1:259321/1‑150 (MQ=255)
                                                                                    gTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAc  >  2:135733/1‑150 (MQ=255)
                                                                                                                   ||                                                                                                                        
CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGAACCTCCGCCAACTGCATTCGTGAACCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAC  >  NC_002947/4741172‑4741404

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