Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,131,519 T→C 7.7% intergenic (‑116/‑128) recQ ← / → PP_4517 ATP‑dependent DNA 3'‑5' helicase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,131,5190TC 7.7% 143.1 / 3.0 52intergenic (‑116/‑128)recQ/PP_4517ATP‑dependent DNA 3'‑5' helicase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base C (2/2);  total (27/25)
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.23e-01

CTTGAGAACGCGCTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATGTCCTTCGCCGAG  >  NC_002947/5131374‑5131661
                                                                                                                                                 |                                                                                                                                              
cTTGAGAACGCGCTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGa                                                                                                                                              >  4:106594/1‑148 (MQ=255)
cTTGAGAACGCGCTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGa                                                                                                                                              >  6:220077/1‑148 (MQ=255)
          cgcTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGa                                                                                                                                   <  1:66296/149‑1 (MQ=255)
          cgcTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGa                                                                                                                                   <  2:2332/149‑1 (MQ=255)
              gAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGTACGAAAAACCCGACCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGta                                                                                                                                >  1:232367/1‑148 (MQ=255)
                   tGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGa                                                                                                                                   >  3:142360/1‑140 (MQ=255)
                   tGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGa                                                                                                                                   <  4:142360/140‑1 (MQ=255)
                        gAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCAc                                                                                                                     <  3:106594/149‑1 (MQ=255)
                          gCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCg                                                                                                                   >  2:85329/1‑149 (MQ=255)
                               tgCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGtgctg                                                                                                              <  3:86334/149‑1 (MQ=255)
                                    cTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACAt                                                                                                         <  1:85329/149‑1 (MQ=255)
                                        aaaTCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACTCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCAc                                                                                                     <  5:220077/149‑1 (MQ=255)
                                                      tCCTGGTTAAAGGCTGGACGAAAACCCCGCCCCGCACACGCCAAATGCTTAATCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATTAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACg                                                                                       <  7:119425/149‑1 (MQ=255)
                                                          ggATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACGGCGGCGATAGACGTTcc                                                                                   >  2:270512/1‑149 (MQ=255)
                                                               aaGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCttt                                                                               <  7:283705/148‑1 (MQ=255)
                                                                   cTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGtgctg                                                                                                              >  7:285930/1‑113 (MQ=255)
                                                                   cTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGtgctg                                                                                                              <  8:285930/113‑1 (MQ=255)
                                                                     ggACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCgcgg                                                                                                >  1:350411/1‑125 (MQ=255)
                                                                     ggACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCgcgg                                                                                                <  2:350411/125‑1 (MQ=255)
                                                                       aCGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCccgctg                                                                      >  8:59556/1‑149 (MQ=255)
                                                                       aCGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCccgctg                                                                      >  4:297923/1‑149 (MQ=255)
                                                                            aaaCCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCccgctgccgct                                                                 >  8:82716/1‑149 (MQ=255)
                                                                                   ccccgcACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGt                                                          <  4:215045/149‑1 (MQ=255)
                                                                                           acGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGAc                                                                                        <  8:171635/111‑1 (MQ=255)
                                                                                           acGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGAc                                                                                        >  7:171635/1‑111 (MQ=255)
                                                                                             gCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTc                                                >  6:122571/1‑149 (MQ=255)
                                                                                                    gCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCg                                         >  2:35796/1‑149 (MQ=255)
                                                                                                           cGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTACACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTc                                  >  8:59362/1‑149 (MQ=255)
                                                                                                                   gCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAgg                           >  3:176215/1‑148 (MQ=255)
                                                                                                                       ggCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGcc                      >  3:38779/1‑149 (MQ=255)
                                                                                                                       ggCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGCTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGcc                      <  2:232367/149‑1 (MQ=255)
                                                                                                                            ccATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACTTGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTa                                                      <  7:68078/112‑1 (MQ=255)
                                                                                                                            ccATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTa                                                      >  8:68078/1‑112 (MQ=255)
                                                                                                                            ccATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGcc                      >  7:195123/1‑144 (MQ=255)
                                                                                                                            ccATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGcc                      <  8:195123/144‑1 (MQ=255)
                                                                                                                            ccATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGt                  >  1:343042/1‑148 (MQ=255)
                                                                                                                             cATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTAACAAAGGTAGCCCTCTaa                >  2:45152/1‑149 (MQ=255)
                                                                                                                               ttGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCg                                               >  4:296870/1‑116 (MQ=255)
                                                                                                                               ttGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCg                                               <  3:296870/116‑1 (MQ=255)
                                                                                                                                  aaaCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATGTcc           <  2:343042/149‑1 (MQ=255)
                                                                                                                                  aaaCCTGCCCCTGGACGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTc                                  >  7:76670/1‑126 (MQ=255)
                                                                                                                                  aaaCCTGCCCCTGGACGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTc                                  >  7:81499/1‑126 (MQ=255)
                                                                                                                                  aaaCCTGCCCCTGGACGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTc                                  <  8:76670/126‑1 (MQ=255)
                                                                                                                                  aaaCCTGCCCCTGGACGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTc                                  <  8:81499/126‑1 (MQ=255)
                                                                                                                                     ccTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATGTCCTTc        <  7:59556/149‑1 (MQ=255)
                                                                                                                                     ccTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATGTCCTTc        <  3:297923/149‑1 (MQ=255)
                                                                                                                                      cTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATGTCCTTCg       >  5:104925/1‑149 (MQ=255)
                                                                                                                                         ccccTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATGTCCTTCGcc     <  4:176215/148‑1 (MQ=255)
                                                                                                                                         ccccTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATGTCCTTCGc      <  4:38779/147‑1 (MQ=255)
                                                                                                                                           ccTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCACAAGGTAGCCCTGTAATGTCCTTCGCCGAg  >  8:293722/1‑149 (MQ=255)
                                                                                                                                            cTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATg              >  4:141925/1‑136 (MQ=255)
                                                                                                                                            cTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATg              <  3:141925/136‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CTTGAGAACGCGCTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGGTAGCCCTGTAATGTCCTTCGCCGAG  >  NC_002947/5131374‑5131661

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