Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,433,665 C→T intergenic (‑105/‑25) smf ← / → def protein Smf/peptide deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,433,6650CT95.5% 64.8 / NA 22intergenic (‑105/‑25)smf/defprotein Smf/peptide deformylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/14);  minor base G (0/1);  total (7/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTT  >  NC_000913/3433421‑3433869
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                            
ccaaaTGTAAGCCCTTTTTGCTGAAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCTTATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCAGCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCATAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGa                                                                                                                                                                                                        <  2:168041/248‑1 (MQ=255)
                        aaTACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCAGTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTTTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCa                                                                                                                                                                                <  1:257839/251‑1 (MQ=255)
                                  caccAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAg                                                                                                                                                                      <  2:204371/248‑1 (MQ=255)
                                         atGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCAGCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTAc                                                                                                                                                                <  2:84044/250‑1 (MQ=255)
                                            cGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCACTCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGACTAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACatat                                                                                                                                                            >  1:97477/1‑251 (MQ=255)
                                                tGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGCTTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCg                                                                                                                                                        >  2:188000/1‑251 (MQ=255)
                                                   tttGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGa                                                                                                                                                      <  1:91513/250‑1 (MQ=255)
                                                                    gATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCCGAGGGGGATTTGTCTAGTATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGt                                                                                                                                    >  1:267838/1‑251 (MQ=255)
                                                                          gACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTa                                                                                                                               >  2:92918/1‑250 (MQ=255)
                                                                             cATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAAcc                                                                                                                           <  1:92918/251‑1 (MQ=255)
                                                                                  catcGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTa                                                                                                                       <  1:63115/250‑1 (MQ=255)
                                                                                  catcGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATCTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTag                                                                                                                      <  1:234918/251‑1 (MQ=255)
                                                                                        cGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCGTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTCGCTAAACCGGTAGaagaag                                                                                                                <  1:269447/251‑1 (MQ=255)
                                                                                              aGCTGCTGATACTAATTAAACGCAGCAAAATGTCTGTAGCGACTATCCTTATCTCCCTGTCAGAAGCAGCCTTAGCAATCTTTGCGAGTGGTCAGTGATGCGGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGGTTTGCAAGTGTTACTTATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATg                                                                                                          <  1:175142/251‑1 (MQ=255)
                                                                                                        aCTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTc                                                                                                 <  1:32379/250‑1 (MQ=255)
                                                                                                                                         cctCCTTATCACCCTGCCGTAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGGCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATTTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGTTAGAAGAAGTGAATGCTGAAATTCAGCGTATCGTCGATGATATTTTCGAGACGATGTa                                                               <  2:113248/249‑1 (MQ=255)
                                                                                                                                          aTCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTGGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAAGTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCTGTAGAAGAAGTGAATGGAGAAATTCAGCTTATCGTCGATGATATGTTCGAGACGATGTAc                                                              <  2:97477/251‑1 (MQ=255)
                                                                                                                                                                         cAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAAccc                               >  1:2123/1‑251 (MQ=255)
                                                                                                                                                                               ttGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGtt                          <  1:133343/250‑1 (MQ=255)
                                                                                                                                                                                          cAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCa               >  1:87158/1‑250 (MQ=255)
                                                                                                                                                                                                      cAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGCAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGtt  >  2:226497/1‑251 (MQ=255)
                                                                                                                                                                                                          cAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCg                                                                                                                          <  1:109306/127‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                            
CGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTT  >  NC_000913/3433421‑3433869

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