Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,474,996 T→C 8.2% intergenic (‑202/‑104) proA ← / → PP_4812 gamma‑glutamyl phosphate reductase/3‑methyladenine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,474,9960TC 8.2% 137.5 / 3.2 50intergenic (‑202/‑104)proA/PP_4812gamma‑glutamyl phosphate reductase/3‑methyladenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (26/19);  new base C (2/2);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

GCCCAGCGGCGGGTGGCGGATGGTAGAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGACAGCTTTTTCGACCGCGACGCCCAGA  >  NC_002947/5474853‑5475142
                                                                                                                                               |                                                                                                                                                   
gCCCAGCGGCGGGTGGCGGATGGTAGAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTAccc                                                                                                                                                >  4:193968/1‑149 (MQ=255)
  ccAGCGGCGGGTGGCGGATGGTAGAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCcgc                                                                                                                                              >  1:47549/1‑149 (MQ=255)
  ccAGCGGCGGGTGGCGGATGGTAGAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCGTTAGCGGGCGCGCCCGCTCCCACGGGGAGTTGTGTAGGAGCGGGTTTGCCcgc                                                                                                                                              >  4:316054/1‑149 (MQ=255)
      cggcggGTGGCGGATGGTAGAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGa                                                                                                                                            >  2:348406/1‑147 (MQ=255)
         cggGTGGCGGATGGTAGAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCCCTCCCACTGGGATATGTGTTGGAGCGGGTTTACCAGCGAAAggg                                                                                                                                       >  5:41965/1‑149 (MQ=255)
               gCGGATGGTAGAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCGCCCACTGGGAGTTGTGTAGGAGCGGGTTTACCCGCGACAGGGCCTGCa                                                                                                                                 >  4:245302/1‑149 (MQ=255)
                      gTAGAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCaa                                                                                                                          <  7:303642/149‑1 (MQ=255)
                         gAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAt                                                                                                                       >  3:294865/1‑149 (MQ=255)
                          aaaaTGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCACCAt                                                                                                                       >  3:228357/1‑148 (MQ=255)
                                    gggCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAAt                                                                                                           <  2:47549/149‑1 (MQ=255)
                                     ggCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTc                                                                                                               >  1:190868/1‑144 (MQ=255)
                                     ggCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTc                                                                                                               <  2:190868/144‑1 (MQ=255)
                                      gCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTc                                                                                                         <  4:74942/149‑1 (MQ=255)
                                      ccAATGTGTAGCAAAGATGTAACCGGGCGTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTCGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTc                                                                                                         <  4:359601/148‑1 (MQ=255)
                                        aaTGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAAc                                                                                                                         >  8:248217/1‑132 (MQ=255)
                                        aaTGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAAc                                                                                                                         <  7:248217/132‑1 (MQ=255)
                                                   aaGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACAGGGATTTGTGTAGGAGCGCCTATACCCCCGAAAGGGACGGCACTGCCAACAT‑CACCGTACAAATCAGCCTCGATTAAg                                                                                            >  4:285918/1‑149 (MQ=255)
                                                    aGATGTAACCGGGGTTCTCGTGATTTCCATGTCCTGTCGCGGCCCTTTAGCGGGCGCGCCGGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGc                                                                                           <  6:41965/149‑1 (MQ=255)
                                                         tAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTTGTATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCcatt                                                                                       <  4:90413/148‑1 (MQ=255)
                                                             cGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTg                                                                                  <  4:311444/149‑1 (MQ=255)
                                                               gggTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCCCTGGGCTTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCGCATTCAGCCTCGAATAAGCCATTCATTGc                                                                                 >  6:328187/1‑148 (MQ=39)
                                                               gggTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGc                                                                                 >  2:333804/1‑148 (MQ=255)
                                                                  ttATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTatc                                                                             <  3:316054/149‑1 (MQ=255)
                                                                      cGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCCAACAT‑CACCGTCCAATTCATCCTCGATTAAGCCATTCAGTGCTatcatca                                                                         >  7:132917/1‑148 (MQ=255)
                                                                      cGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTa                                                                               >  2:99608/1‑143 (MQ=255)
                                                                      cGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTa                                                                               <  1:99608/143‑1 (MQ=255)
                                                                       gagaTTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATccc                                                                        >  7:280355/1‑149 (MQ=255)
                                                                       gagaTTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGGTTACGCGCGAAAGGGCCGGCACTGGCAACAT‑CTCCGTGCAGTTCAGGCTCGATTAAGCCATTCTTTGCTTTCATccc                                                                        >  7:275246/1‑149 (MQ=255)
                                                                       gagaTTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAACGGGGTTAACCGCGAAACGGCCGGCACGGCCAACATGCACC‑TGCAATTAATCCTAAATTAAGCCAATCAATGGTATAATccc                                                                        >  1:112510/1‑149 (MQ=255)
                                                                        agaTTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCAATTCATTGCTATCATcccc                                                                       >  2:62654/1‑149 (MQ=255)
                                                                          aTTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGc                                                                     <  1:348406/149‑1 (MQ=255)
                                                                          aTTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCACCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGc                                                                     <  5:328187/149‑1 (MQ=255)
                                                                                            ggCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTg                                                   <  1:333804/149‑1 (MQ=255)
                                                                                                cTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGcccc                                               >  6:60122/1‑149 (MQ=255)
                                                                                                           gcgcCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGcc                                    >  3:287510/1‑149 (MQ=255)
                                                                                                            cgcCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGccc                                   <  3:193968/149‑1 (MQ=255)
                                                                                                               ccGCTCCCACTGGGATTTGTGTAGGCGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGc                                <  3:194833/149‑1 (MQ=255)
                                                                                                                     ccACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGAc                           <  3:157955/148‑1 (MQ=255)
                                                                                                                       aCTGGGATTTGTGTAGGAGCGGGTCTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAAc                                                     <  2:247464/120‑1 (MQ=255)
                                                                                                                       aCTGGGATTTGTGTAGGAGCGGGTCTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAAc                                                     <  2:252836/120‑1 (MQ=255)
                                                                                                                       aCTGGGATTTGTGTAGGAGCGGGTCTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAAc                                                     >  1:247464/1‑120 (MQ=255)
                                                                                                                       aCTGGGATTTGTGTAGGAGCGGGTCTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAAc                                                     >  1:252836/1‑120 (MQ=255)
                                                                                                                             aTTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGACAGCTTTTTc                  <  3:285918/149‑1 (MQ=255)
                                                                                                                                 gtgtAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGACAGCTTTTTCGAcc              >  4:173984/1‑149 (MQ=255)
                                                                                                                                  tgtAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGACAGCTTTTTCGACcg             >  2:24561/1‑149 (MQ=255)
                                                                                                                                      ggAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGACAGCTTTTTCGACCGCGAc         >  3:160478/1‑149 (MQ=255)
                                                                                                                                             gTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCCGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGACAGCTTTTTCGACCGCGACGCCCAGa  >  2:85420/1‑149 (MQ=255)
                                                                                                                                             gTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTa                                                                               <  6:177467/72‑1 (MQ=255)
                                                                                                                                             gTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTa                                                                               >  5:177467/1‑72 (MQ=255)
                                                                                                                                             gTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGACAGCTTTTTCGACCGCGACGCCCAGa  >  2:206349/1‑149 (MQ=255)
                                                                                                                                              tttACCCGCGTAACGGCCGGGAATTCCATCAT‑CACCGTCCAATTCAGCCTCAATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGACAGCTTTTTCGACCGCGACGCCCAGa  <  2:112510/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCCCAGCGGCGGGTGGCGGATGGTAGAAAATGGATGGGGCAATGTGTAGGAAAGATGTAACCGGGGTTATCGAGATTTCCATCTCCTGTCGCGGCCCTTTAGCGGGCGCGCCCGCTCCCACTGGGATTTGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAACAT‑CACCGTCCAATTCAGCCTCGATTAAGCCATTCATTGCTATCATCCCCGCCTCTCCAGCCACGAACTGCCCCGCATGCCAGCCCTGCCCGACAGCTTTTTCGACCGCGACGCCCAGA  >  NC_002947/5474853‑5475142

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