Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,313,871 T→G 22.8% 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,8710TG22.8% 68.3 / 20.6 41intergenic (+33/‑126)hemB/proP5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease
Reads supporting (aligned to +/- strand):  ref base T (12/19);  new base G (4/5);  total (16/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 = 7.93e-01

CGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAA  >  NC_002947/3313724‑3314015
                                                                                                                                                   |                                                                                                                                                
cGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGc                                                                                                                                                >  2:825636/1‑150 (MQ=255)
  gTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAg                                                                                                                                              <  2:487577/150‑1 (MQ=255)
      cATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGc                                                                                                                                          <  2:859366/150‑1 (MQ=255)
       aTCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCGATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGACGCTGCCATTGCAGCGGCt                                                                                                                                         >  1:339813/1‑150 (MQ=255)
        tCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCtt                                                                                                                                        <  2:1042099/150‑1 (MQ=255)
          gacgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATAAAAGCCGCTGCAATGGCAGCGGCtttc                                                                                                                                      <  2:873326/150‑2 (MQ=255)
            cgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCt                                                                                                                                    <  1:1106350/150‑1 (MQ=255)
             gaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttcat                                                                                                                                    <  2:897516/149‑4 (MQ=255)
             gaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttcat                                                                                                                                    >  1:897516/1‑146 (MQ=255)
                       tggtgCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGt                                                                                                                         <  1:282511/150‑1 (MQ=255)
                       tggtgCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGt                                                                                                                         >  2:282511/1‑150 (MQ=255)
                         gtgCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTg                                                                                                                       <  1:334428/150‑1 (MQ=255)
                         gtgCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTg                                                                                                                       <  1:177895/150‑1 (MQ=255)
                         gtgCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTg                                                                                                                       <  1:177857/150‑1 (MQ=255)
                                aaaGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCAGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGgtggtg                                                                                                                >  2:688324/1‑150 (MQ=255)
                                                 aaGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACaaaaa                                                                                               >  2:128546/1‑150 (MQ=255)
                                                     gTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACa                                                                                           <  2:902296/150‑1 (MQ=255)
                                                                  ccTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAat                                                                              >  2:951588/1‑150 (MQ=255)
                                                                    tGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAatta                                                                            <  1:480581/150‑1 (MQ=255)
                                                                        ccTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTa                                                                        <  1:951588/150‑1 (MQ=255)
                                                                        ccTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTa                                                                        <  2:708192/150‑1 (MQ=255)
                                                                         cTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAg                                                                       >  2:1084422/1‑150 (MQ=255)
                                                                         cTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAg                                                                       >  1:26634/1‑150 (MQ=255)
                                                                             ccaaCTACGCAATGGACCTGGCGCGCGAATGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCAGCTTTTCTTTGTGCGCCGTTGGGTGGTGAAAAAAAAAAAAAAAAAAACAGTAAGCGCTGAATTATTTAGTGTc                                                                   <  2:660091/147‑1 (MQ=255)
                                                                                  ttCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt                                                              >  2:818788/1‑150 (MQ=255)
                                                                                           gACCTGGCGCGCGAAGGTATCTGAAGCCAGGCGGAAGGCATGAAAGCCGCTGCCACAGCAGCGAAATTTCTTTATGAACCGTTGGTTTGTAAAAAAAAAAAAAAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGAt                                                     <  2:320629/150‑1 (MQ=255)
                                                                                                ggCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCtt                                                >  1:892076/1‑150 (MQ=255)
                                                                                                 gcgcgcGAAGATACCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTa                                               <  2:923422/150‑1 (MQ=255)
                                                                                                              tCTGAAGCCAGGCCTTGGTCATGAAATCCGCTGCAATGGCAGCGGc                                                                                                                                          >  2:480577/1‑46 (MQ=21)
                                                                                                              tCTGAAGCCAGGCCTTGGTCATGAAATCCGCTGCAATGGCAGCGGc                                                                                                                                          <  1:480577/46‑1 (MQ=21)
                                                                                                                     ccAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGCTGAGTGGTGCCGGCCTCTTCCAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGAt                           <  1:688324/150‑1 (MQ=255)
                                                                                                                      cAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTATTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTa                                               <  2:342679/129‑1 (MQ=255)
                                                                                                                      cAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTATTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTa                                               >  1:342679/1‑129 (MQ=255)
                                                                                                                         tccatGGTCATGAAAGCCTCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCACAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCACCGACAGGGATACgg                       <  2:339813/146‑1 (MQ=255)
                                                                                                                          ccTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACggg                      <  1:1084422/150‑1 (MQ=255)
                                                                                                                                  aTGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATGGACAGGGATACGGGTATGAAAc              >  2:559284/1‑150 (MQ=255)
                                                                                                                                      aaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTg                                                                                                                       >  2:745221/1‑41 (MQ=17)
                                                                                                                                      aaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTg                                                                                                                       <  1:745221/41‑1 (MQ=17)
                                                                                                                                        gCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCa        >  1:395138/1‑150 (MQ=255)
                                                                                                                                         ccGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCaa       <  1:818788/150‑1 (MQ=255)
                                                                                                                                              gCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGaaaa  <  1:128546/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAA  >  NC_002947/3313724‑3314015

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: