Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,649,393 C→G 9.2% intergenic (+870/+821) PP_3214 → / ← PP_3216 hypothetical protein/cysteine dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,649,3930CG 9.2% 162.4 / 5.4 52intergenic (+870/+821)PP_3214/PP_3216hypothetical protein/cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (22/25);  new base G (2/3);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCGACCACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAG  >  NC_002947/3649245‑3649530
                                                                                                                                                     |                                                                                                                                         
ccGTCGACCACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCAc                                                                                                                                           >  7:1372/1‑149 (MQ=255)
     gACCACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGggc                                                                                                                                      >  1:147809/1‑149 (MQ=255)
        cacaGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGggcggc                                                                                                                                   >  6:36577/1‑149 (MQ=255)
              gACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCg                                                                                                                             <  6:222081/149‑1 (MQ=255)
               aCCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCgg                                                                                                                            <  5:36577/149‑1 (MQ=255)
                 cTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa                                                                                                                          <  2:147809/149‑1 (MQ=255)
                 cTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa                                                                                                                          <  8:1372/149‑1 (MQ=255)
                    ccAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGTCGGGCACCGGGCGGCTGCCCGGCAGAt                                                                                                                       >  3:295573/1‑149 (MQ=255)
                    ccAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCAGCCGCCCGGTGCCCGGcgg                                                                                                                         <  2:123321/147‑3 (MQ=255)
                    ccAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGAt                                                                                                                       >  3:94068/1‑149 (MQ=255)
                        cccATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGAcc                                                                                                                   >  7:164382/1‑149 (MQ=255)
                            tCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGAGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGgcg                                                                                                               >  8:259232/1‑149 (MQ=255)
                                    agccaCTAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTc                                                                                                       <  6:247678/149‑1 (MQ=255)
                                      ccaCGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGTCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCt                                                                                                     <  4:295573/149‑1 (MQ=255)
                                       caCGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTc                                                                                                    <  8:164382/149‑1 (MQ=255)
                                        aCGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGt                                                                                                           >  3:205458/1‑141 (MQ=255)
                                        aCGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGt                                                                                                           <  4:205458/141‑1 (MQ=255)
                                           aGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTc                                                                                                <  7:135427/149‑1 (MQ=255)
                                               gcTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa                                                                                                                          >  7:213566/1‑119 (MQ=255)
                                               gcTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa                                                                                                                          <  8:213566/119‑1 (MQ=255)
                                                   acgatGGCGTGCGCGCCATCGTCGTCCACCAGCGCGATGCGCCCCCCGCGC‑ACACTGGCGACGAGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCCGCGCTGTATTCCTCGGTCCTACGCgg                                                                                        >  1:308120/5‑149 (MQ=255)
                                                        ggCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGc                                                                                   >  5:334435/1‑149 (MQ=255)
                                                         gcgcgcACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCt                                                                                  <  6:334435/149‑1 (MQ=255)
                                                          cgcgcACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTa                                                                                 <  3:74271/149‑1 (MQ=255)
                                                               aCGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGcc                                                                            >  2:64733/1‑149 (MQ=255)
                                                                 gCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCg                                                                          >  5:247546/1‑149 (MQ=255)
                                                                    aTCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGc                                                                       <  6:247546/149‑1 (MQ=255)
                                                                        tcgtcCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTgg                                                                   <  7:259232/149‑1 (MQ=255)
                                                                             caccaGCGCGATGCGCGCCCCGCGCAACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCg                                                                                                   >  1:3337/1‑113 (MQ=255)
                                                                             caccaGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCg                                                                                                   <  2:3337/112‑1 (MQ=255)
                                                                             caccaGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCcggc                                                              <  1:64733/149‑1 (MQ=255)
                                                                                caGCGCGATGCGCGCGCCGCGC‑ACTCTGGGGACGTGGTGGGTCTCTTGAACCAGTTGCCGGCCGGGGAGGGGGCGGCTGCCCGGGAGATGACCGGCGCTGTATTCCTCGGTCGTACGGGGTGTGCTACCGCCCGGGCGTGGCCggcggc                                                           <  2:308120/149‑1 (MQ=255)
                                                                                 aGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGAc                                                                                                                    >  2:316161/1‑91 (MQ=255)
                                                                                 aGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGAc                                                                                                                    <  1:316161/91‑1 (MQ=255)
                                                                                           gcgcCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCggg                                                <  3:80499/149‑1 (MQ=255)
                                                                                                ccGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTAcc                                                                               <  8:202439/113‑1 (MQ=255)
                                                                                                ccGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTAcc                                                                               >  7:202439/1‑113 (MQ=255)
                                                                                                       acacTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCAGCCGCCCGGTGCCCGGcg                                                                                                                          <  4:38634/64‑2 (MQ=21)
                                                                                                       acacTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCAGCCGCCCGGTGCCCGGcg                                                                                                                          >  3:38634/1‑63 (MQ=21)
                                                                                                           tGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCt                                 >  8:43606/1‑149 (MQ=255)
                                                                                                             gCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTc                                <  3:307932/148‑1 (MQ=255)
                                                                                                               gACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTg                                                                    >  6:1270/1‑110 (MQ=255)
                                                                                                               gACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTg                                                                    <  5:1270/110‑1 (MQ=255)
                                                                                                                             ctTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCt               >  8:39791/1‑149 (MQ=255)
                                                                                                                              ttGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGcc                                                                      >  7:146950/1‑93 (MQ=255)
                                                                                                                              ttGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGcc                                                                      <  8:146950/93‑1 (MQ=255)
                                                                                                                              ttGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAgcg                       <  4:206248/140‑1 (MQ=255)
                                                                                                                              ttGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAgcg                       >  3:206248/1‑140 (MQ=255)
                                                                                                                                gAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGa             <  3:88391/148‑1 (MQ=255)
                                                                                                                                     aGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGcccct         >  2:278286/1‑145 (MQ=255)
                                                                                                                                     aGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGa       <  7:39791/149‑1 (MQ=255)
                                                                                                                                          cccGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAg  >  5:154282/1‑149 (MQ=255)
                                                                                                                                             gccgGGCAGCCGCCCGGTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAgcg                       >  4:329142/1‑125 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
CCGTCGACCACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAG  >  NC_002947/3649245‑3649530

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