Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,118,719 G→A E69K (GAA→AAA)  zapB → cell division factor ZapB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,118,7190GA95.2% 51.3 / NA 21E69K (GAA→AAA) zapBcell division factor ZapB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/11);  minor base C (1/0);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCC  >  NC_000913/4118473‑4118958
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
cGAGGATTTTTATCCATCAACACCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAATAACAGCAGAACGGCTGGCAGAAACg                                                                                                                                                                                                                                             <  1:140392/251‑1 (MQ=255)
                                     gagGTATGACAATGTCATTAGAAGTGTTTGAAAACCTGGCAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGCAATCGCAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAgg                                                                                                                                                                                                             >  1:430854/1‑246 (MQ=255)
                                                                        cTGGANNCANNNGTNNAGCAGGCGATTGATACCATCACTNTGTTGCAGATGGAAATCGAAGAGCTGAAANANANAAACNNCTCNCTGTCGCAGGNAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAANCNTCTGAAAGAACAGCAGAACGGCTGGCAGAAANGTCTGCAGGccc                                                                                                                                                                                                                                  >  2:364854/1‑190 (MQ=255)
                                                                        cTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAATAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGccc                                                                                                                                                                                                                                  <  1:364854/190‑1 (MQ=255)
                                                                                aaaaGTACAGCAGGCGATTGATACCCTCACTCTGTTGCAGATGGAAATCGCAGCACTGGAAGACAAAAACAACTCACTGTCGCAGGAAGTTCAACATGAACAGCATCAGCGCGCAGAGCTGGAGCGTGAGAACAAACATCTGAAAGAACAGCAGAACGACTGGCAGAAACGTCTGCa                                                                                                                                                                                                                                       >  1:449570/1‑177 (MQ=255)
                                                                                            ngCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGAt                                                                                                                                                                                                       >  1:316733/2‑197 (MQ=255)
                                                                                                                     taGATGGAAATCAAAGAGCGGAAATAAAAAAACAACTCACGGTGGCAGGAAGTTCAAAACGCCCCGCATCAGCGCGAAGATCTTGAGCGTTTGAAAAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGc                                                                                                                                                                                                                                    <  1:359573/142‑1 (MQ=255)
                                                                                                                           gAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGGCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATc                                                                                                                  <  1:206501/251‑1 (MQ=255)
                                                                                                                           gAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTAGCAGGAAGTTCTAAGTGCCAAGAATCAGCGGGAAGAGCTCGAGCGTGAGCCCAACCATCTGATAGCACAGCAGAACGGCTGGCAGAAACGTCTGCCGGACctgctg                                                                                                                                                                                                                             >  1:329079/1‑144 (MQ=255)
                                                                                                                                       nnnnnnnnnnaannNNNCTCNCTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGAt                                                                                                                                                                                                       <  2:316733/144‑1 (MQ=255)
                                                                                                                                         gaaagaaaAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGCACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGTAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTa                                                                                                                                                                          >  1:499074/1‑181 (MQ=255)
                                                                                                                                                          cACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGNNAGAAACNNc                                                                                                                                                                                                                                           >  2:39093/1‑99 (MQ=255)
                                                                                                                                                          cACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACtg                                                                                                                                                                       <  1:39093/167‑1 (MQ=255)
                                                                                                                                                                   aGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCGGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCGTGACTGGGTAGGCAAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCGGAAAGAATAATCCCggg                                                                          <  1:145812/251‑2 (MQ=255)
                                                                                                                                                                                          cAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACTTCTGCAGGCCCTGCTGGGTTGCATGGAAGATGTGTTATTACGGTTCATCTTCCGCGTAATTTCTTTACTGGGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTttatta                                                                                                                    <  1:451325/186‑1 (MQ=255)
                                                                                                                                                                                           aGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGCCCTGCTgg                                                                                                                                                                                                                            <  1:377997/81‑1 (MQ=255)
                                                                                                                                                                                            gcgcGAAGGGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGGAGAAACTTCTTCAGGCCCTGCTGGGTCGTATTGAATAGGGCTGATTACTCTTCATCTTCCGCGTAATTTCTTGACGGTGTAGGCCAACAGGTAACGCa                                                                                                                                                  <  1:452791/154‑1 (MQ=255)
                                                                                                                                                                                                   gagCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAgg                                                                                                                                                                 <  1:444107/132‑1 (MQ=255)
                                                                                                                                                                                                            ggTGAGAACAACCATCTGAAAGAACAGCAGAACGGGTGGCAGCAACGTCTGCAGGCCCTGCTGGGGCGCATGGAAGAGGGCTGATTACGCTTCATCTTCAGCGTAATTTCTGTACTGTGTAGGCCTCCAGGTCACGCAGAGACAAGGCCCCTTGCGGTGCCTGTCTTATCATTCAATATCCAGCGGATCTTCTGAACGGCTGATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGCAGACGGTGaca                                  >  1:88681/2‑250 (MQ=255)
                                                                                                                                                                                                                                           aCGGCTGGCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGcc  >  1:4818/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 ggCAGAAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTc                                                                    >  1:213220/1‑179 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
CGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCC  >  NC_000913/4118473‑4118958

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