Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,683,251 C→G 8.2% intergenic (+60/‑171) recJ → / → PP_1478 ssDNA‑specific exonuclease RecJ/xenobiotic reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,683,2510CG 8.2% 140.2 / 4.7 49intergenic (+60/‑171)recJ/PP_1478ssDNA‑specific exonuclease RecJ/xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base G (2/2);  total (26/23)
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.13e-01

CGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGAT  >  NC_002947/1683115‑1683392
                                                                                                                                        |                                                                                                                                              
cGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                    >  3:79340/1‑149 (MQ=255)
 gTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                   >  1:336817/1‑149 (MQ=255)
   aCAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                        <  8:17033/142‑1 (MQ=255)
   aCAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                        >  7:17033/1‑142 (MQ=255)
        cTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAgg                                                                                                                            >  1:116932/1‑149 (MQ=255)
           gATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGgcgc                                                                                                                         <  7:261518/149‑1 (MQ=255)
           gATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGgcg                                                                                                                          <  2:264063/148‑1 (MQ=255)
                 aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGAtct                                                                                                                                  <  3:117741/134‑1 (MQ=255)
                 aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGAtct                                                                                                                                  >  4:117741/1‑134 (MQ=255)
                  aCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGAtctc                                                                                                                                 <  6:279164/134‑1 (MQ=255)
                  aCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGAtctc                                                                                                                                 >  5:279164/1‑134 (MQ=255)
                   cGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAcc                                                                                                                 >  3:113899/1‑149 (MQ=255)
                    gAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCt                                                                                                                <  4:79340/149‑1 (MQ=255)
                            cgGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCg                                                                                                        >  3:241081/1‑148 (MQ=255)
                             ggTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGaa                                                                                                      >  3:210140/1‑149 (MQ=255)
                                  tGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAAc                                                                                                     <  7:307444/145‑1 (MQ=255)
                                  tGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAAc                                                                                                     >  8:307444/1‑145 (MQ=255)
                                    agatTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTc                                                                                                <  1:45034/148‑1 (MQ=255)
                                       gtgtGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCa                                                                                            >  3:101282/1‑149 (MQ=255)
                                          tgCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATcc                                                                                         <  2:336817/149‑1 (MQ=255)
                                            cAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCCCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCca                                                                                       >  2:268872/1‑149 (MQ=255)
                                                 gatgatTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCaaa                                                                                                                   <  1:350706/117‑1 (MQ=255)
                                                        gcgcATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTccc                                                                            <  2:276022/148‑1 (MQ=255)
                                                          gcATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTccccc                                                                          >  7:188555/1‑148 (MQ=255)
                                                            atatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCccac                                                                       >  1:226313/1‑149 (MQ=255)
                                                                 gAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATc                                                                  >  6:11702/1‑149 (MQ=255)
                                                                      gcgcTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTc                                                             <  8:188555/149‑1 (MQ=255)
                                                                         cTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCCCACAGGCGCTGCAAACCTTTCCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTa                                                           <  1:268872/149‑1 (MQ=255)
                                                                            aGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCGGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTATtc                                                        <  4:241081/148‑1 (MQ=255)
                                                                              tttAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGtc                                                        >  2:255765/1‑146 (MQ=255)
                                                                              tttAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGtc                                                        <  1:255765/146‑1 (MQ=255)
                                                                                   cGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCcca                                                                        >  3:296672/1‑125 (MQ=255)
                                                                                   cGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCcca                                                                        <  4:296672/125‑1 (MQ=255)
                                                                                             ctgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCAcc                                      >  6:319710/1‑149 (MQ=255)
                                                                                              tgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCg                                     >  5:225894/1‑149 (MQ=255)
                                                                                                  ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTa                                                  <  8:200936/132‑1 (MQ=255)
                                                                                                  ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTa                                                  >  7:200936/1‑132 (MQ=255)
                                                                                                     gATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAg                               <  4:113899/148‑1 (MQ=255)
                                                                                                           gcgcTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAt                           <  5:11702/146‑1 (MQ=255)
                                                                                                                  ggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                 <  4:210140/149‑1 (MQ=255)
                                                                                                                   gCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGGCTAGTCTCATTAATGCCCGGCACCGGGCCAGGGATGTGCACGTGAg                >  4:293501/1‑149 (MQ=255)
                                                                                                                    cacaAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt               >  8:292006/1‑149 (MQ=255)
                                                                                                                      caAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                 <  4:305023/145‑1 (MQ=255)
                                                                                                                      caAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                 >  3:305023/1‑145 (MQ=255)
                                                                                                                        aaTCGAGCGCCGCCCGCGCGGCGCTCGATCTAACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTcc           >  1:218128/1‑149 (MQ=255)
                                                                                                                          tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATc                                                                  <  4:254640/92‑1 (MQ=255)
                                                                                                                          tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCAACCCACAGAACCTCCCCCCACCAATc                                                                  >  3:254640/1‑92 (MQ=255)
                                                                                                                          tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCgg         <  3:293501/149‑1 (MQ=255)
                                                                                                                              gcgcCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCg     >  4:101419/1‑149 (MQ=255)
                                                                                                                                 ccgcccgcGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGCACCCCTCGCCATCCCACAGACCCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGAt  >  7:140202/1‑149 (MQ=255)
                                                                                                                                 ccgcccgcGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGc                                  >  4:60356/1‑117 (MQ=255)
                                                                                                                                 ccgcccgcGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGc                                  <  3:60356/117‑1 (MQ=255)
                                                                                                                                       gcgcgGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                 >  7:145973/1‑128 (MQ=255)
                                                                                                                                       gcgcgGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                 <  8:145973/128‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                              
CGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTT‑CGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGAT  >  NC_002947/1683115‑1683392

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