Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,639,596 C→T 15.2% intergenic (‑186/‑125) czcR‑III ← / → PP_1439 response regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,639,5960CT15.2% 119.5 / 12.0 55intergenic (‑186/‑125)czcR‑III/PP_1439response regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base T (4/4);  total (29/26)
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.04e-01

TAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGACCAAATACCTGCTG  >  NC_002947/1639452‑1639744
                                                                                                                                                |                                                                                                                                                    
tAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTcgcg                                                                                                                                                  <  4:215897/149‑1 (MQ=37)
tAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                                                   <  1:335722/148‑1 (MQ=255)
tAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                                                   <  1:340243/148‑1 (MQ=255)
 gatctcATCAGGGTCATGCGCGCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTcgcg                                                                                                                                                  <  6:62664/143‑1 (MQ=25)
 aGTGACATCAGGGTCAAGCGCCACAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTcgcg                                                                                                                                                  <  8:125168/148‑1 (MQ=25)
       aTCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                                                   <  8:190796/141‑1 (MQ=255)
       aTCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                                                   >  7:190796/1‑141 (MQ=255)
       aTCAGGGTCAAGCGCCCCAATGTTAAGGCATGAAAAAAGGCCTGCGGAGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCACTTCATACAGGCGCGGGTTTACTCGCGAagatac                                                                                                                                           >  5:62664/1‑146 (MQ=25)
          aGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGcac                                                                                                                                         >  3:299861/1‑146 (MQ=255)
             gTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                                                   >  6:296916/1‑135 (MQ=255)
             gTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                                                   <  5:296916/135‑1 (MQ=255)
                 aGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTCCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                                                 <  7:214883/149‑1 (MQ=255)
                  gcgcCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCGTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                                                 >  3:342518/1‑148 (MQ=255)
                  gcgcCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                                                 >  1:197375/1‑148 (MQ=255)
                     caccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACTGATTTGTAATTCTGCTGTTATCTGGCTGGCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAATAGGCCCGAACAGTCcaca                                                                                                                             <  6:170553/147‑1 (MQ=255)
                        cAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacac                                                                                                                            >  6:129559/1‑147 (MQ=255)
                             ttACGGCTTGAAAAACTGCCTGCGGTGCATATAAGGTTTTTTTATTTATGGTGTTATGTCGCTGTCTGTTGGGGCCTGTTCGCGAGTAAACCCGGGCGTGTAGGTACGGGTTAACCGGCTATGAGGCCCGAACAGTCCACACAAGGcc                                                                                                                      <  6:337894/148‑1 (MQ=255)
                              tACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTGca                                                                                                                               >  3:154307/1‑138 (MQ=255)
                                cGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACAt                                                                                                                  <  5:129559/149‑1 (MQ=255)
                                        aaaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGAc                                                                                                          <  1:90661/149‑1 (MQ=255)
                                              gCCTGAGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACACGCGCCTGTAGGAACGGGTTTACACGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAAt                                                                                                    >  5:281632/1‑149 (MQ=255)
                                                        cAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCAGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGCCTGTAATCCCCGTCTCa                                                                                          >  2:98789/1‑149 (MQ=255)
                                                         aGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCAc                                                                                         >  6:97726/1‑149 (MQ=255)
                                                               cGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCCCCTTGcc                                                                                   >  3:18980/1‑149 (MQ=255)
                                                               cGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACACGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGcc                                                                                   >  7:328351/1‑149 (MQ=255)
                                                                 gATTTCTAATTCTGCTGTTATCTGGCGGTCTGTTCGGGCCGCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCCAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGt                                                                                 <  4:124334/149‑1 (MQ=255)
                                                                      gTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCca                                                                                                                               >  8:350772/1‑98 (MQ=255)
                                                                          ttCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagccagcc                                                                        >  4:325714/1‑149 (MQ=255)
                                                                               ctgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCg                                                                    <  5:97726/148‑1 (MQ=255)
                                                                                tgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTTCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCATCGcg                                                                  >  8:344740/1‑148 (MQ=255)
                                                                                tgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGc                                                                   >  2:50485/1‑148 (MQ=255)
                                                                                             tctgtTCGGGCCTCTTCGCGCGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGCAGAGGCCCGAACAGTGCCCCCACGGCCACATTGCAGGACTGTAATCCCCGCCTGCCCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgtgt                                                     >  7:33534/1‑149 (MQ=255)
                                                                                               tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATCGCCGTGTCt                                                   >  1:316935/1‑149 (MQ=255)
                                                                                                 ttCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGa                                                 >  7:26629/1‑149 (MQ=255)
                                                                                                       ccTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgtgt                                                     >  1:108774/1‑139 (MQ=255)
                                                                                                       ccTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgtgt                                                     <  2:108774/139‑1 (MQ=255)
                                                                                                           ttCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttatta                                       >  5:292710/1‑149 (MQ=25)
                                                                                                             cgcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTAc                                      <  3:325714/148‑1 (MQ=255)
                                                                                                              gcgGGTAAACCAGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGa                                    <  6:292710/149‑1 (MQ=37)
                                                                                                              gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGa                                    <  2:197375/149‑1 (MQ=255)
                                                                                                                gggTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAgg                                  >  4:168682/3‑149 (MQ=37)
                                                                                                                 aGTAAACCCGCGCATGTAGGAACGGGTTTACCCGCGAAGAGGCCCTAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc                                 <  8:328351/149‑1 (MQ=255)
                                                                                                                       cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc                                 >  5:265902/1‑143 (MQ=255)
                                                                                                                       cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc                                 <  6:265902/143‑1 (MQ=255)
                                                                                                                        ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac                           >  3:251956/1‑148 (MQ=255)
                                                                                                                        ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac                           >  4:42680/1‑148 (MQ=255)
                                                                                                                         cgcgCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAGCGGTTATTACGAGGCAAACACGa                         <  6:281632/149‑1 (MQ=255)
                                                                                                                              cTGTAGGAACTGGTTTACCCTCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAAc                    <  4:299861/149‑1 (MQ=255)
                                                                                                                              cTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAAc                    <  8:26629/149‑1 (MQ=255)
                                                                                                                                    gAAGGGGTTTACGCGCGAAGAGGCCCGAACAGTCCACTCAAGGCCACTTTGCTGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGAc               <  7:344740/148‑1 (MQ=255)
                                                                                                                                    gAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGAt                        >  6:223639/1‑139 (MQ=255)
                                                                                                                                    gAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGAt                        <  5:223639/139‑1 (MQ=255)
                                                                                                                                      aCGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGAACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCGCGCTAGGATGGCCGTGTCTGATCGGTTATTACCAGGCAAACAGGATGACCCTGACCa             >  1:92947/1‑148 (MQ=255)
                                                                                                                                       cGGGTTTACCCGCCATTAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGACCaaa           <  4:18980/149‑1 (MQ=255)
                                                                                                                                          gTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGACCAAATAc        >  5:290381/1‑149 (MQ=255)
                                                                                                                                                ccGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGACCAAATACctgctg  >  5:342964/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGACCAAATACCTGCTG  >  NC_002947/1639452‑1639744

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: