Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,433,900 C→T L71L (CTA→TTA)  def → peptide deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,433,9000CT84.9% 65.5 / 4.3 34L71L (CTA→TTA) defpeptide deformylase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (15/13);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

CCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGCTGAAAC  >  NC_000913/3433667‑3434143
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
ccTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATGCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTAGCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCGGGCGGCAACCCAGGTGGATATCCATCAACGTATCATTGTTATGGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCa                                                                                                                                                                                                                                     <  1:21306/251‑1 (MQ=255)
             ctGGAGATTTATGTCAGTTTTGCAAGTGTTCCATATTCCGGACGATCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGAGGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAAGCAATCCAGAGCTTTTAg                                                                                                                                                                                                                           <  1:2080/248‑1 (MQ=255)
               ggAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGaa                                                                                                                                                                                                                         <  1:179681/248‑1 (MQ=255)
                gagaTTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAAAGGTTAGTGTTAATCAATCCAGAGCTTTTAGaaaaa                                                                                                                                                                                                                      >  2:85868/1‑250 (MQ=255)
                                                  ccGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGAAGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGCAGAAGGTTGCCTGTc                                                                                                                                                                                   >  1:96972/1‑251 (MQ=255)
                                                          gCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGATGCCTGTCGATCCCTg                                                                                                                                                                           >  1:153804/1‑251 (MQ=255)
                                                                 cGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGGTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTAGTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAAGCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGaacaac                                                                                                                                                                     <  2:159939/250‑1 (MQ=255)
                                                                  gCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCTGATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGTCGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGt                                                                                                                                                                   <  2:82076/251‑1 (MQ=255)
                                                                         tGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTa                                                                                                                                                             <  1:270387/250‑1 (MQ=255)
                                                                                  ggTAGAAGAAGTGAAGGCAGAAATTCAGCTTTTCGTCGATGATATGTACGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAATCCGTGACGAACTGTTAGTGTTAATCAAACCAGAGCGTTTAGAAAAAAGCGGCGAATCTGGCATTGAAGAAGGTTGCCTGGCGATCGCTTAACAACGTGCTTTAGTGCcgcgcg                                                                                                                                                   <  2:199005/251‑1 (MQ=255)
                                                                                    tAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCgcgcgc                                                                                                                                                  >  1:54517/1‑250 (MQ=255)
                                                                                               aaTGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAGACGGCGAAACAGGAATTGAAGAAGGTTGCCTGTCGATCCCTGAACACCGGGCTTTAGTGCCGCACGCAGAGAAAGTTa                                                                                                                                      >  2:104165/1‑250 (MQ=255)
                                                                                               aaTGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTa                                                                                                                                      >  1:354889/1‑250 (MQ=255)
                                                                                                     gaaaagcagctgatcttcgATGATATGGTCGAGACTATGTCAGCAGAAGAAGGTATTGTCCTGGCGGCACCCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGTAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTATTATAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTAGCCTGTACATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGGTAAAATTcg                                                                                                                               <  2:117785/235‑1 (MQ=255)
                                                                                                              cGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGa                                                                                                                      <  1:380996/251‑1 (MQ=255)
                                                                                                               gTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTaa                                                                                                                                                                                                                                             <  1:228398/132‑1 (MQ=255)
                                                                                                               gTATCGTCGATGATATGTTAGAGACGATGTACGCAGAAGAAGGTATTAGCCTGGCGGCAACCCAGGATGTTATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTaa                                                                                                                                                                                                                                             >  2:228398/1‑132 (MQ=255)
                                                                                                                  ccgtcgATGATAGGGTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAATCCAGGTTGATATCCATAAACGTATCATTGTTATTGAGGTTTCGGAAAACCGTGACGAACGGGTAGGGTTAATCAATCCAGAGCTTTTAGAAGAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGGCGATCCCTGAACAACGTGCGTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACcgc                                                                                                                  <  2:171485/250‑1 (MQ=255)
                                                                                                                    gtcgATGATATGTTCGAGACGATGTACGCAGAAGAAGGAATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGAGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGAGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCACGCAGAGAAAGTTAAAATTCGCGCACTTGAccccgc                                                                                                                >  2:26978/1‑247 (MQ=255)
                                                                                                                        atgatATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTATCGGAAAACCGTGACTAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGGTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCAGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGTGACGGt                                                                                                            <  2:50203/251‑1 (MQ=255)
                                                                                                                                   gagaCGATGTACGCAGAAGACGGTCTTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAACGCGGCGAAACAGGCATTGTAGACGGTTGCCTGTCGATCCCTGCACAACGTGCTTTAGTGCCGCGCGCAGGACAAGTTAAAATTCGGGCCATTGACCGCGACGGTAAACCATTTg                                                                                                  >  2:131802/1‑250 (MQ=255)
                                                                                                                                       cGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTa                                                                                                                                                                                                                            <  1:161588/125‑1 (MQ=255)
                                                                                                                                       cGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTa                                                                                                                                                                                                                            >  2:161588/1‑125 (MQ=255)
                                                                                                                                               gCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGACGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAAAAACGTGCTTTAGGGCGGGGCGCAGGGAAGGTTAAAATTCGCGCCTTTGACCGCGACGGTAAACAATTTGAACTGGAAGc                                                                                        >  1:3456/1‑248 (MQ=255)
                                                                                                                                                 agaataAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACg                                                                                   <  2:100574/251‑1 (MQ=255)
                                                                                                                                                             tGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGATAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGcc                                                                        >  1:209295/1‑250 (MQ=255)
                                                                                                                                                                           cccAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCAt                                                         >  1:282100/1‑251 (MQ=255)
                                                                                                                                                                                  tGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATgg                                                  >  1:14335/1‑251 (MQ=255)
                                                                                                                                                                                    atatCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGAt                                                >  1:221596/1‑251 (MQ=255)
                                                                                                                                                                                                         ttATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAg                                                                                      >  1:85772/1‑192 (MQ=255)
                                                                                                                                                                                                         ttATTGAGGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCGGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAg                                                                                      <  2:85772/192‑1 (MQ=255)
                                                                                                                                                                                                                  tttCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTAt                  >  1:46854/1‑251 (MQ=255)
                                                                                                                                                                                                                              gTGACGAACGGTTAGTTTTAATCATGCCAGAGCTTTTAGAAAAAAGCTGCGAAACAGGCATGGAAGAAGGTTGCCTGTCGATCCCTTAACAACGTGCTTTAGTGTCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGTAATCAGACGGTCTGTTATCCATCTGTATTCAGCATGAGAGGGATCACCTGGTCTGCAAACTTTTTATGGATTATCTGTCACCGCTg      <  2:39863/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cGAACGGCGAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGTTGAaac  <  2:96972/251‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
CCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGCTGAAAC  >  NC_000913/3433667‑3434143

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