Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,861,950 G→A 8.2% D15D (GAC→GAT)  purM ← phosphoribosylformylglycinamidine cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,861,9500GA 8.2% 140.8 / 2.8 50D15D (GAC→GAT) purMphosphoribosylformylglycinamidine cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base A (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 = 9.13e-01

GGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTGGTCTGTGAAG  >  NC_002947/1861803‑1862094
                                                                                                                                                     |                                                                                                                                                 
ggTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGt                                                                                                                                                  >  4:74248/1‑149 (MQ=255)
    ccGGAGACCAACACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTa                                                                                                                                               <  5:50636/148‑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                                                                                                                                       <  3:240474/118‑1 (MQ=255)
                                          ctcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCg                                                                                                                                       >  4:240474/1‑118 (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              <  4:91978/149‑1 (MQ=255)
                                                                                                                                     cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAg              <  7:186852/149‑1 (MQ=255)
                                                                                                                                        cGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTgg           >  7:247000/1‑149 (MQ=255)
                                                                                                                                                 cGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTGGTCTGTGAAg  >  2:186192/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTGGTCTGTGAAG  >  NC_002947/1861803‑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: ‑