Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,283,900 A→G S61P (TCC→CCC)  purU ← formyltetrahydrofolate deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,283,9000AG100.0% 85.4 / NA 25S61P (TCC→CCC) purUformyltetrahydrofolate deformylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/13);  total (12/13)

TTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATATTGGTAATACGTGCGATCAGA  >  NZ_CP009273/1283764‑1284027
                                                                                                                                        |                                                                                                                               
ttGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGAt                                                                                                                               <  1:304415/139‑1 (MQ=255)
     tttCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCAtt                                                                                                                           >  1:350211/1‑138 (MQ=255)
            aacaaATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATc                                                                                                                    <  1:120012/138‑1 (MQ=255)
                   cGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAg                                                                                                            >  1:553145/1‑139 (MQ=255)
                         ggCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGt                                                                                                      <  2:280478/139‑1 (MQ=255)
                                    tcttTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATaaaa                                                                                            >  1:526193/1‑138 (MQ=255)
                                          gTGACCAGAATCTCTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGAAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                     <  2:553145/139‑1 (MQ=255)
                                          gTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                     <  2:260671/139‑1 (MQ=255)
                                             aCCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                     >  1:389235/1‑136 (MQ=255)
                                             aCCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                     <  2:389235/136‑1 (MQ=255)
                                               cAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                     >  1:243541/1‑134 (MQ=255)
                                               cAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                     <  2:243541/134‑1 (MQ=255)
                                               cAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTa                                                                                >  2:404783/1‑139 (MQ=255)
                                                       cTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                     >  1:547955/1‑126 (MQ=255)
                                                       cTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                     <  2:547955/126‑1 (MQ=255)
                                                            cGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATc                                                                                                                    <  1:488686/90‑1 (MQ=255)
                                                            cGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATc                                                                                                                    >  2:488686/1‑90 (MQ=255)
                                                            cGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATcaaca                                                                    >  1:394418/1‑138 (MQ=255)
                                                                            gATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTAc                                                   <  2:350211/139‑1 (MQ=255)
                                                                                                             cgcTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTa                                                    <  1:535112/105‑1 (MQ=255)
                                                                                                             cgcTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTa                                                    >  2:535112/1‑105 (MQ=255)
                                                                                                             cgcTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATAtt                     <  1:536706/136‑1 (MQ=255)
                                                                                                             cgcTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATAtt                     >  2:536706/1‑136 (MQ=255)
                                                                                                                            ccAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATATTGGTAATACGTGCGAt      >  2:304503/1‑136 (MQ=255)
                                                                                                                             cagcagGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATATTGGTAATACGTGCGATCAGa  <  1:561704/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                               
TTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATATTGGTAATACGTGCGATCAGA  >  NZ_CP009273/1283764‑1284027

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