Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,649,397 G→C 8.0% intergenic (+874/+817) 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,3970GC 8.0% 156.8 / 4.7 50intergenic (+874/+817)PP_3214/PP_3216hypothetical protein/cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base C (2/2);  total (22/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 = 4.10e-01

GACCACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTTGCT  >  NC_002947/3649250‑3649535
                                                                                                                                                    |                                                                                                                                          
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                                                                                                                               <  8:1372/149‑1 (MQ=255)
            cTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa                                                                                                                               <  2:147809/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                                                                                                                         <  1:316161/91‑1 (MQ=255)
                                                                            aGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGAc                                                                                                                         >  2:316161/1‑91 (MQ=255)
                                                                                      gcgcCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCggg                                                     <  3:80499/149‑1 (MQ=255)
                                                                                           ccGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTAcc                                                                                    >  7:202439/1‑113 (MQ=255)
                                                                                           ccGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTAcc                                                                                    <  8:202439/113‑1 (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)
                                                                                                                                                  gggCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTtgct  <  7:171019/141‑1 (MQ=255)
                                                                                                                                                  gggCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTtgct  >  8:171019/1‑141 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GACCACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTTGCT  >  NC_002947/3649250‑3649535

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: ‑