Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,861,946 C→A 7.5% D17Y (GAC→TAC)  purM ← phosphoribosylformylglycinamidine cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,861,9460CA 7.5% 166.9 / 3.5 53D17Y (GAC→TAC) purMphosphoribosylformylglycinamidine cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base A (2/2);  total (26/27)
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.90e-01

CCGTCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCTGTGAAG  >  NC_002947/1861798‑1862094
                                                                                                                                                      |                                                                                                                                                     
ccGTCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTc                                                                                                                                                       <  4:69856/149‑1 (MQ=255)
  gTCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGa                                                                                                                                                     >  5:269441/1‑149 (MQ=255)
   tCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGAt                                                                                                                                                    <  3:176412/149‑1 (MQ=255)
     ggTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGt                                                                                                                                                  >  4:74248/1‑149 (MQ=255)
         ccGGAGACCAACACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTa                                                                                                                                               <  5:50636/148‑1 (MQ=255)
            gagaCCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGAt                                                                                                                                                    >  8:128644/1‑140 (MQ=255)
            gagaCCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGAt                                                                                                                                                    <  7:128644/140‑1 (MQ=255)
                 cAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTacaagatag                                                                                                                                       >  3:182441/1‑142 (MQ=255)
                 cAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTac                                                                                                                                              <  4:182441/141‑1 (MQ=255)
                 cAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGt                                                                                                                                      >  4:99051/1‑149 (MQ=255)
                 cAGCACCGGCTGCTTGTAGCCGGCAGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGAGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCAGGCGTCGATCTCTACACCGGCga                                                                                                                                      >  7:255499/1‑148 (MQ=255)
                   gCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTcc                                                                                                                                    <  6:269441/149‑1 (MQ=255)
                   gCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATCACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTCCGCCGGCGTCGATGTCTACACCGGCGTcc                                                                                                                                    <  1:24326/149‑1 (MQ=255)
                      ccGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATTACTTCAGGGCGTGCGGTGCGTTTCAGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCtt                                                                                                                                  <  1:96389/149‑1 (MQ=255)
                      ccGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTg                                                                                                                                 <  3:299635/149‑1 (MQ=255)
                       cGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGt                                                                                                                                <  6:191315/149‑1 (MQ=255)
                        ggCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTa                                                                                                                               <  2:264732/149‑1 (MQ=255)
                             ctTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCATTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCa                                                                                                                          <  2:60201/149‑1 (MQ=255)
                                         cggGATCTCGCAGAGGGCGCCGAATCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgct                                                                                                              >  7:47069/1‑149 (MQ=255)
                                            gATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgcttgc                                                                                                           <  3:74248/149‑1 (MQ=255)
                                            gATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGACGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgcttgc                                                                                                           >  1:75510/1‑149 (MQ=255)
                                               ctcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCg                                                                                                                                       <  3:240474/118‑1 (MQ=255)
                                               ctcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCg                                                                                                                                       >  4:240474/1‑118 (MQ=255)
                                                  gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCAtt                                                                                                      <  4:183466/148‑1 (MQ=255)
                                                  gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCAt                                                                                                       <  8:47069/147‑1 (MQ=255)
                                                  gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTg                                                                                                     <  8:173320/149‑1 (MQ=255)
                                                         gcgcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGagctaca                                                                                              <  3:99051/149‑7 (MQ=255)
                                                           gcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGTGTAGACATCGACGCCgg                                                                                                                                         >  5:150227/1‑104 (MQ=255)
                                                           gcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGTGTAGACATCGACGCCgg                                                                                                                                         <  6:150227/104‑1 (MQ=255)
                                                                 gccgccCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTag                                                                                       >  1:30421/1‑148 (MQ=255)
                                                                    gccCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAggaggg                                                                                  >  4:157206/1‑149 (MQ=255)
                                                                               ccATGACTTCAGGGCGTGCGGTGCGTTGTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgc                                                                                                               >  1:177875/1‑111 (MQ=255)
                                                                               ccATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGCCGAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCt                                                                         >  3:247337/1‑148 (MQ=255)
                                                                               ccATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgc                                                                                                               <  2:177875/110‑1 (MQ=255)
                                                                               ccATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCt                                                                         >  8:66294/1‑148 (MQ=255)
                                                                                      ttCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAAc                                                                  >  3:91978/1‑148 (MQ=255)
                                                                                         aGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGATTCGCCGGCGTCGATGTCTACACCGGCGTCATTGTAGCTCAGGGAGGGTTGCTTGCTAATTGATCCAGGCCTGTAGGAGGGAGGGATTCTTAAAAACGAc                                                               >  5:175876/1‑148 (MQ=255)
                                                                                            gCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCAt                                                            >  3:309876/1‑148 (MQ=255)
                                                                                                                 aCGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGa                                                                >  5:146766/1‑125 (MQ=255)
                                                                                                                 aCGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGa                                                                <  6:146766/125‑1 (MQ=255)
                                                                                                                   gCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATaa                                      <  7:66294/149‑1 (MQ=255)
                                                                                                                              gTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGt                           >  6:153309/1‑149 (MQ=255)
                                                                                                                                  gACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTg                       >  8:186852/1‑148 (MQ=255)
                                                                                                                                   aCCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGaa                     <  4:309876/149‑1 (MQ=255)
                                                                                                                                    ccAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGaaa                    <  5:153309/149‑1 (MQ=255)
                                                                                                                                    ccAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGaaa                    <  2:75510/149‑1 (MQ=255)
                                                                                                                                     cAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAAc                   >  8:178033/1‑149 (MQ=255)
                                                                                                                                       gTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGt                 >  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)
                                                                                                                                          cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAg              <  7:186852/149‑1 (MQ=255)
                                                                                                                                          cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAg              <  4:91978/149‑1 (MQ=255)
                                                                                                                                             cGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTgg           >  7:247000/1‑149 (MQ=255)
                                                                                                                                                      cGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCTGTGAAg  >  2:186192/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCGTCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTTTTGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCTGTGAAG  >  NC_002947/1861798‑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: