Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,621,890 T→C S232P (TCA→CCA)  purM → phosphoribosylformylglycinamide cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,621,8900TC100.0% 129.5 / NA 36S232P (TCA→CCA) purMphosphoribosylformylglycinamide cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/16);  total (20/16)

AGCGGTCCACACTCGAACGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCA  >  NC_000913/2621755‑2622020
                                                                                                                                       |                                                                                                                                  
aGCGGTCCACACTCGAACGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGt                                                                                                                                >  2:501265/1‑140 (MQ=255)
   ggTCCACACTCGAACGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCt                                                                                                                             <  2:907741/140‑1 (MQ=255)
      ccACACTCGAACGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGa                                                                                                                          >  2:909787/1‑140 (MQ=255)
      ccACACTCGAACGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGa                                                                                                                          >  2:891922/1‑140 (MQ=255)
               aaCGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGttgattga                                                                                                                 >  2:219025/1‑140 (MQ=255)
               aaCGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGttgattga                                                                                                                 >  1:566263/1‑140 (MQ=255)
               aaCGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTAATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGttgattga                                                                                                                 >  1:772722/1‑140 (MQ=255)
                aCGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGaa                                                                                                                <  1:909787/140‑1 (MQ=255)
                      aTTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTc                                                                                                          >  2:45469/1‑140 (MQ=255)
                      aTTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTc                                                                                                          >  2:895451/1‑140 (MQ=255)
                      aTTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTc                                                                                                          >  2:393061/1‑140 (MQ=255)
                       ttCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCg                                                                                                         <  1:219025/140‑1 (MQ=255)
                         cGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGAt                                                                                                       <  1:891922/140‑1 (MQ=255)
                                                           tgATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCt                                                                     >  2:406275/1‑140 (MQ=255)
                                                            gATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCtt                                                                    <  1:406275/140‑1 (MQ=255)
                                                                  cAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGa                                                              >  1:128/1‑140 (MQ=255)
                                                                   aaaCCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGaa                                                             >  1:742672/1‑140 (MQ=255)
                                                                                 gATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGt                                               >  2:390587/1‑140 (MQ=255)
                                                                                     gTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTg                                           <  2:742672/140‑1 (MQ=255)
                                                                                       aaGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGcc                                         <  2:128/140‑1 (MQ=255)
                                                                                       aaGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGcc                                         <  1:895451/140‑1 (MQ=255)
                                                                                       aaGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGcc                                         <  1:45469/140‑1 (MQ=255)
                                                                                       aaGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGcc                                         <  1:393061/140‑1 (MQ=255)
                                                                                        aGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCa                                        >  2:139992/1‑140 (MQ=255)
                                                                                                    atcatcTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACt                               >  2:185907/1‑137 (MQ=255)
                                                                                                    atcatcTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACt                               <  1:185907/137‑1 (MQ=255)
                                                                                                       atcTGCTGGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCa                         <  1:139992/140‑1 (MQ=255)
                                                                                                         ctgctgGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACt                               <  2:183713/132‑1 (MQ=255)
                                                                                                         ctgctgGCACCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACt                               >  1:183713/1‑132 (MQ=255)
                                                                                                                 aCCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATtgatg               >  2:211470/1‑140 (MQ=255)
                                                                                                                 aCCGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATtgatg               >  2:559029/1‑140 (MQ=255)
                                                                                                                  ccGACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATtgatga              >  1:866452/1‑140 (MQ=255)
                                                                                                                    gACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAAt            >  1:913950/1‑140 (MQ=255)
                                                                                                                    gACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAAt            <  2:566263/140‑1 (MQ=255)
                                                                                                                    gACCCGCATTTACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAAt            <  1:390587/140‑1 (MQ=255)
                                                                                                                              tACGTGAAGCCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCa  <  1:211470/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
AGCGGTCCACACTCGAACGGCTATTCGCTGGTGCGCAAAATTCTTGAAGTCAGCGGTTGTGATCCGCAAACCACCGAACTTGATGGTAAGCCATTAGCCGATCATCTGCTGGCACCGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCA  >  NC_000913/2621755‑2622020

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