Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,313,871 T→G 19.4% intergenic (+33/‑126) hemB → / → proP 5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,313,8710TG19.4% 46.1 / 3.6 24intergenic (+33/‑126)hemB/proP5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease
Reads supporting (aligned to +/- strand):  ref base T (7/12);  new base G (2/3);  total (9/15)
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.00e+00

CATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTC  >  NC_002947/3313730‑3314006
                                                                                                                                             |                                                                                                                                       
cATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGc                                                                                                                                 <  1:520480/150‑1 (MQ=255)
 aTCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCt                                                                                                                                >  2:361296/1‑150 (MQ=255)
                       gcgAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGggtg                                                                                                          <  1:361296/150‑1 (MQ=255)
                                  ggCGCAATCAAGCGTGCTGGTGCAGCCAGGAGCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAAc                                                                                               <  1:309501/150‑1 (MQ=255)
                                   gcgcAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACg                                                                                              >  2:540415/1‑150 (MQ=255)
                                             gCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGAGCCGTTGGGTGGTGCCGCAAACGCACaaaaaaa                                                                                    <  2:587203/150‑1 (MQ=255)
                                                            ccTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAat                                                                     >  2:577049/1‑150 (MQ=255)
                                                                        cTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTcc                                                         <  1:540415/150‑1 (MQ=255)
                                                                               gCCATGGCCCTGGCGCTCTAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGACGCGGCGGTGATCGCGGTGGTGCTTTATGCGCCGTTGGGTGGTGCCGCACACGCTCAATAAAACAGTATGCGCTGAATTATTTAGTGTCCTGTTGTa                                                  <  1:294391/150‑1 (MQ=255)
                                                                                ccATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTAt                                                 >  2:519619/1‑150 (MQ=255)
                                                                                       ccTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATcc                                          <  1:207776/150‑1 (MQ=255)
                                                                                         tGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCt                                        <  1:577049/150‑1 (MQ=255)
                                                                                             gcgcgAAGGTATCTGAAGCCAGGCCTTGGTAATGAGAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTaga                                    <  1:526178/150‑1 (MQ=255)
                                                                                                              gCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGa                   <  1:578015/150‑1 (MQ=255)
                                                                                                              gCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGa                   >  2:578015/1‑150 (MQ=255)
                                                                                                                             tGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTa                                                  >  2:40207/1‑104 (MQ=255)
                                                                                                                             tGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTa                                                  <  1:40207/104‑1 (MQ=255)
                                                                                                                               aaaGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTc  <  2:576645/150‑1 (MQ=255)
                                                                                                                               aaaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTAt          >  2:4745/1‑142 (MQ=255)
                                                                                                                               aaaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTAt          <  1:4745/142‑1 (MQ=255)
                                                                                                                                 aGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAat                                                                     >  2:200965/1‑81 (MQ=255)
                                                                                                                                 aGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAat                                                                     <  1:200965/81‑1 (MQ=255)
                                                                                                                                  gCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGAc                        <  2:21036/125‑1 (MQ=255)
                                                                                                                                  gCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGAc                        >  1:21036/1‑125 (MQ=255)
                                                                                                                                             |                                                                                                                                       
CATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTC  >  NC_002947/3313730‑3314006

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: