Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,861,949 T→C 8.0% I16V (ATC→GTC)  purM ← phosphoribosylformylglycinamidine cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,861,9490TC 8.0% 140.7 / 3.2 50I16V (ATC→GTC) purMphosphoribosylformylglycinamidine cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

TCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTGGTCTGTGAAG  >  NC_002947/1861801‑1862094
                                                                                                                                                      |                                                                                                                                                  
tCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGAt                                                                                                                                                    <  3:176412/149‑1 (MQ=255)
  ggTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGt                                                                                                                                                  >  4:74248/1‑149 (MQ=255)
      ccGGAGACCAACACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTa                                                                                                                                               <  5:50636/148‑1 (MQ=255)
         gagaCCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGAt                                                                                                                                                    >  8:128644/1‑140 (MQ=255)
         gagaCCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGAt                                                                                                                                                    <  7:128644/140‑1 (MQ=255)
              cAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTacaagatag                                                                                                                                       >  3:182441/1‑142 (MQ=255)
              cAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTac                                                                                                                                              <  4:182441/141‑1 (MQ=255)
              cAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGt                                                                                                                                      >  4:99051/1‑149 (MQ=255)
              cAGCACCGGCTGCTTGTAGCCGGCAGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGAGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCAGGCGTCGATCTCTACACCGGCga                                                                                                                                      >  7:255499/1‑148 (MQ=255)
                gCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTcc                                                                                                                                    <  6:269441/149‑1 (MQ=255)
                gCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATCACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTCCGCCGGCGTCGATGTCTACACCGGCGTcc                                                                                                                                    <  1:24326/149‑1 (MQ=255)
                   ccGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATTACTTCAGGGCGTGCGGTGCGTT‑TCAGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCtt                                                                                                                                  <  1:96389/149‑1 (MQ=255)
                   ccGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTg                                                                                                                                 <  3:299635/149‑1 (MQ=255)
                    cGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGt                                                                                                                                <  6:191315/149‑1 (MQ=255)
                     ggCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTa                                                                                                                               <  2:264732/149‑1 (MQ=255)
                          ctTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCATT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCa                                                                                                                          <  2:60201/149‑1 (MQ=255)
                                      cggGATCTCGCAGAGGGCGCCGAATCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgct                                                                                                              >  7:47069/1‑149 (MQ=255)
                                         gATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgcttgc                                                                                                           <  3:74248/149‑1 (MQ=255)
                                         gATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGACGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgcttgc                                                                                                           >  1:75510/1‑149 (MQ=255)
                                            ctcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCg                                                                                                                                       >  4:240474/1‑118 (MQ=255)
                                            ctcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCg                                                                                                                                       <  3:240474/118‑1 (MQ=255)
                                               gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCAtt                                                                                                      <  4:183466/148‑1 (MQ=255)
                                               gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCAt                                                                                                       <  8:47069/147‑1 (MQ=255)
                                               gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTg                                                                                                     <  8:173320/149‑1 (MQ=255)
                                                      gcgcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGagctaca                                                                                              <  3:99051/149‑7 (MQ=255)
                                                        gcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGTGTAGACATCGACGCCgg                                                                                                                                         >  5:150227/1‑104 (MQ=255)
                                                        gcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGTGTAGACATCGACGCCgg                                                                                                                                         <  6:150227/104‑1 (MQ=255)
                                                              gccgccCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTag                                                                                       >  1:30421/1‑148 (MQ=255)
                                                                 gccCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGA‑CAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAggaggg                                                                                  >  4:157206/1‑149 (MQ=255)
                                                                            ccATGACTTCAGGGCGTGCGGTGCGTTGT‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgc                                                                                                               >  1:177875/1‑111 (MQ=255)
                                                                            ccATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGCCGAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCt                                                                         >  3:247337/1‑148 (MQ=255)
                                                                            ccATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgc                                                                                                               <  2:177875/110‑1 (MQ=255)
                                                                            ccATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCt                                                                         >  8:66294/1‑148 (MQ=255)
                                                                                   ttCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAAc                                                                  >  3:91978/1‑148 (MQ=255)
                                                                                      aGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGATTCGCCGGCGTCGATGTCTACACCGGCGTCATTGTAGCTCAGGGAGGGTTGCTTGCTAATTGATCCAGGCCTGTAGGAGGGAGGGATTCTTAAAAACGAc                                                               >  5:175876/1‑148 (MQ=255)
                                                                                         gCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCAt                                                            >  3:309876/1‑148 (MQ=255)
                                                                                                              aCGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGa                                                                <  6:146766/125‑1 (MQ=255)
                                                                                                              aCGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGa                                                                >  5:146766/1‑125 (MQ=255)
                                                                                                                gCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATaa                                      <  7:66294/149‑1 (MQ=255)
                                                                                                                           gTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGt                           >  6:153309/1‑149 (MQ=255)
                                                                                                                               gACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTg                       >  8:186852/1‑148 (MQ=255)
                                                                                                                                aCCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGaa                     <  4:309876/149‑1 (MQ=255)
                                                                                                                                 ccAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGaaa                    <  2:75510/149‑1 (MQ=255)
                                                                                                                                 ccAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGaaa                    <  5:153309/149‑1 (MQ=255)
                                                                                                                                  cAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAAc                   >  8:178033/1‑149 (MQ=255)
                                                                                                                                    gTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGt                 >  3:195564/1‑149 (MQ=255)
                                                                                                                                     aaggacGCCGGTGTAGACATCGACGCCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCt                                           >  7:296156/6‑123 (MQ=38)
                                                                                                                                     aaggacGCCGGTGTAGACATCGACGCCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCt                                           <  8:296156/118‑1 (MQ=38)
                                                                                                                                       cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGGTTGGTCATTGATCCAGGGCGTTAGGAGGGAGGGATTCTTAAAATCGTCCATGACGGCCAAGGCCGGCTGATAAGCGGCGGGGGTGCTGaagggt                 >  1:272816/1‑143 (MQ=255)
                                                                                                                                       cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAg              <  7:186852/149‑1 (MQ=255)
                                                                                                                                       cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAg              <  4:91978/149‑1 (MQ=255)
                                                                                                                                          cGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTgg           >  7:247000/1‑149 (MQ=255)
                                                                                                                                                   cGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTGGTCTGTGAAg  >  2:186192/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTGGTCTGTGAAG  >  NC_002947/1861801‑1862094

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