Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,639,617 C→T 10.7% intergenic (‑165/‑146) 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,6170CT10.7% 150.1 / 8.4 56intergenic (‑165/‑146)czcR‑III/PP_1439response regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (3/3);  total (29/27)
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.85e-01

TTGAAAGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGAT  >  NC_002947/1639468‑1639764
                                                                                                                                                     |                                                                                                                                                   
ttGAAAGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCcgc                                                                                                                                                     >  1:417564/1‑150 (MQ=255)
       gggACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTa                                                                                                                                              <  1:1141756/150‑1 (MQ=255)
            gATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGCAAc                                                                                                                                         <  2:417564/150‑1 (MQ=255)
                   gACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTa                                                                                                                                  >  2:581223/1‑150 (MQ=255)
                      ccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTAccc                                                                                                                               >  2:807847/1‑150 (MQ=255)
                        tGTAGTGACATCAGGGTCAAGCGACCCAATGTCCGGATCTCAACCTGCGCGTGCCCTATAGATAACGGGTTTATTATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGAAAATCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                             >  1:209817/1‑150 (MQ=255)
                             tGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAaga                                                                                                                        <  1:807847/150‑1 (MQ=255)
                                  tctgggTTAAGCAACCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGccc                                                                                                                   <  2:68743/147‑1 (MQ=18)
                                    aGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGa                                                                                                                 >  2:1047782/1‑150 (MQ=255)
                                     gggTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGaa                                                                                                                <  1:1047782/150‑1 (MQ=255)
                                      ggTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTa                                                                                                                                              <  2:864610/119‑1 (MQ=255)
                                      ggTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTa                                                                                                                                              >  1:864610/1‑119 (MQ=255)
                                        tCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAg                                                                                                             >  2:596473/1‑150 (MQ=255)
                                         cAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTATGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGt                                                                                                            <  1:274143/150‑1 (MQ=255)
                                          aaGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCGGCGAAGAGGCCCTAACAGTc                                                                                                           >  1:205673/1‑150 (MQ=255)
                                                cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacaca                                                                                                     >  1:742053/1‑150 (MQ=255)
                                                cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacaca                                                                                                     >  1:803142/1‑150 (MQ=255)
                                                  cAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAg                                                                                                   <  1:1185167/150‑1 (MQ=255)
                                                          cGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACAtt                                                                                           <  1:581223/150‑1 (MQ=255)
                                                            gCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGc                                                                                         <  2:742053/150‑1 (MQ=255)
                                                              aTGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAg                                                                                       <  2:803142/150‑1 (MQ=255)
                                                                gAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAaga                                                                                                                        <  1:717771/115‑1 (MQ=17)
                                                                gAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAaga                                                                                                                        >  2:717771/1‑115 (MQ=17)
                                                                  aaaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGc                                                                                         <  2:332846/144‑1 (MQ=255)
                                                                  aaaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGc                                                                                         >  1:332846/1‑144 (MQ=255)
                                                                   aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGaa                                                                                                                          >  2:475375/1‑110 (MQ=17)
                                                                   aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGaa                                                                                                                          <  1:475375/110‑1 (MQ=17)
                                                                    aaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTg                                                                                                                                                >  1:581208/1‑87 (MQ=255)
                                                                    aaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTg                                                                                                                                                <  2:581208/87‑1 (MQ=255)
                                                                     aaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTa                                                                                >  2:728931/1‑150 (MQ=255)
                                                                                     aTAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCtt                                                                >  1:764384/1‑150 (MQ=255)
                                                                                          ggATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTg                                                                                  >  2:757838/1‑127 (MQ=255)
                                                                                          ggATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTg                                                                                  <  1:757838/127‑1 (MQ=255)
                                                                                                gTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagcca                                                     <  1:596473/150‑1 (MQ=255)
                                                                                                   aTTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCTCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagccagcc                                                  >  2:52402/1‑150 (MQ=255)
                                                                                                       tgctgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCg                                              >  1:395142/1‑150 (MQ=255)
                                                                                                        gctgTTATCTGGCTGGCTGTTCTGGCCTCTTCGCGAGAGCACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCCGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGc                                             <  2:205673/150‑1 (MQ=255)
                                                                                                          tgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAACCCTCGCTa                                           >  1:76154/1‑150 (MQ=255)
                                                                                                           gTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACaa                                                                                                    <  2:650717/92‑1 (MQ=255)
                                                                                                           gTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACaa                                                                                                    >  1:650717/1‑92 (MQ=255)
                                                                                                                cTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATgg                                     <  2:395142/150‑1 (MQ=255)
                                                                                                                        ctgtTCTGGCCGCTTCACGAGTACACCCGCGCCTGGAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCt                             <  1:911981/150‑1 (MQ=255)
                                                                                                                         tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTa                                           >  2:1036111/1‑135 (MQ=255)
                                                                                                                         tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTa                                           <  1:1036111/135‑1 (MQ=255)
                                                                                                                           ttCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAt                          >  1:98586/1‑150 (MQ=21)
                                                                                                                               gggctctTCGGGAGTAAACCCGCGCCGGTAGGAAAGGTTTTATCCGCGAAGTGGCCCGAACAGTCCACACAAGGCCACATTGCGGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGt                      <  2:209817/147‑1 (MQ=255)
                                                                                                                               ggCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCTCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGt                      <  1:52402/150‑1 (MQ=255)
                                                                                                                                        gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGt                                                           <  2:716631/104‑1 (MQ=255)
                                                                                                                                        gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGt                                                           >  1:716631/1‑104 (MQ=255)
                                                                                                                                        gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAg             >  1:150322/1‑150 (MQ=255)
                                                                                                                                             tAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCaaa        >  2:692625/1‑150 (MQ=255)
                                                                                                                                                aCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAt                          >  2:1011151/1‑129 (MQ=255)
                                                                                                                                                aCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAt                          <  1:1011151/129‑1 (MQ=255)
                                                                                                                                                 cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACg    >  2:267226/1‑150 (MQ=255)
                                                                                                                                                 cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACg    >  1:793234/1‑150 (MQ=255)
                                                                                                                                                  ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAg                                          >  2:1191813/1‑111 (MQ=255)
                                                                                                                                                  ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAg                                          <  1:1191813/111‑1 (MQ=255)
                                                                                                                                                   cgcgCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGAt  <  2:764384/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTGAAAGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGAT  >  NC_002947/1639468‑1639764

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