Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,649,401 C→G 7.7% intergenic (+878/+813) 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,4010CG 7.7% 166.3 / 4.9 52intergenic (+878/+813)PP_3214/PP_3216hypothetical protein/cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (21/27);  new base G (2/2);  total (23/29)
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.97e-01

CACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTTGCT  >  NC_002947/3649253‑3649535
                                                                                                                                                     |                                                                                                                                      
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                                                                                                                               <  8:213566/119‑1 (MQ=255)
                                       gcTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCa                                                                                                                               >  7:213566/1‑119 (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)
                                                                                                                                               gggCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTtgct  <  7:171019/141‑1 (MQ=255)
                                                                                                                                               gggCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTtgct  >  8:171019/1‑141 (MQ=255)
                                                                                                                                                     |                                                                                                                                      
CACAGCGACCTGCCAGCCCATCTGTGCCAGCCACGAGGCGCTCATGTTGGCGCGCACGCCATCGTCGTCCACCAGCGCGATGCGCGCCCCGCGC‑ACACTGGCGACGTGGTCGGTCTCTTGAACCAGTTGCCCGCCGGGCACCGGGCGGCTGCCCGGCAGATGACCGGCGCTGTATTCCTCGGTCGTACGCGGTCTGCTACCGCCCGGGCCTGGCCGGCGGCGTTGGCGCGGGTAACCGCGGTGACCTCGGCAAAGCGGCGGTCCTGACCATGAGCCAGTTGCT  >  NC_002947/3649253‑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: ‑