Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,137,813 A→G 15.4% intergenic (‑164/‑49) PP_4522 ← / → PP_4523 LysR family transcriptional regulator/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,137,8130AG15.4% 58.6 / 3.7 26intergenic (‑164/‑49)PP_4522/PP_4523LysR family transcriptional regulator/agmatinase
Reads supporting (aligned to +/- strand):  ref base A (9/13);  new base G (2/2);  total (11/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 = 9.86e-01

TGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGACTCTCCACCAACCACTGGGCGGCAACGAAATGCCGCGTTTCGGCGGCATCGCCACCATGCTCCGCCTCC  >  NC_002947/5137675‑5137952
                                                                                                                                          |                                                                                                                                           
tgCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCgcggc                                                                                                                                  >  1:631995/1‑150 (MQ=255)
   gcgcAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCggcggc                                                                                                                               >  1:291693/1‑150 (MQ=255)
               tGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCAGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaac                                                                                                                   >  1:528867/1‑150 (MQ=255)
                   ctACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaac                                                                                                               <  1:478230/150‑1 (MQ=255)
                       ccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAGTGGCGCGCGGCGGCACAGGTTCGaacaaa                                                                                                                <  2:254955/145‑1 (MQ=255)
                       ccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAGTGGCGCGCGGCGGCACAGGTTCGaacaaa                                                                                                                >  1:254955/1‑145 (MQ=255)
                            tttGGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACAcaacaacc                                                                                                      <  1:256021/150‑3 (MQ=255)
                                      gaaTTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACtga                                                                                            >  2:151516/2‑150 (MQ=255)
                                          tCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAATAATCGACCGACTGATGAg                                                                                        <  2:334695/150‑1 (MQ=255)
                                          tCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAg                                                                                        <  2:291693/150‑1 (MQ=255)
                                           ccTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCATGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAAAACAACAATCGACCGACTGATGAgg                                                                                       <  2:339447/150‑1 (MQ=255)
                                                 cTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCAcc                                                                                 <  1:380679/150‑1 (MQ=255)
                                                      tATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGa                                                                            <  1:151516/150‑1 (MQ=255)
                                                        tttGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACa                                                                          >  1:115957/1‑150 (MQ=255)
                                                              ggTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGActc                                                                    <  2:631995/150‑1 (MQ=255)
                                                               gTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGActct                                                                   >  1:549499/1‑150 (MQ=255)
                                                                 aGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGACTCTcc                                                                 <  2:281900/150‑1 (MQ=255)
                                                                               ggCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAGTGGCGCGCTGCGGCACAGGTTCGAACAAACACAACAATcga                                                                                                    >  1:92321/1‑101 (MQ=255)
                                                                               ggCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAGTGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcga                                                                                                    <  2:92321/101‑1 (MQ=255)
                                                                                            gCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGACTCTCCACCAACCACTGGGCGGCAACGAAATGc                                      >  1:185497/1‑150 (MQ=255)
                                                                                                      ccGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGACTCTCCACCAACCACTGGGCGGCAACGAAATGCCGCGTTTcgg                            <  2:185497/150‑1 (MQ=255)
                                                                                                          aaTGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACAcaacaa                                                                                                        <  2:306970/70‑1 (MQ=255)
                                                                                                          aaTGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACAcaacaa                                                                                                        >  1:306970/1‑70 (MQ=255)
                                                                                                                                ccACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGACTCTCCACCAACCACTGGGCGGCAACGAAATGCCGCGTTTCGGCGGCATCGCCACCATGCTCCGCCTcc  <  2:115957/150‑1 (MQ=255)
                                                                                                                                    acGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGACTCTCCACCAACCACTGGGCGGCAACGa                                           >  1:711159/1‑105 (MQ=255)
                                                                                                                                    acGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGACTCTCCACCAACCACTGGGCGGCAACGa                                           <  2:711159/105‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
TGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGACTGATGAGGCCCACCGTGGACAAGACTCTCCACCAACCACTGGGCGGCAACGAAATGCCGCGTTTCGGCGGCATCGCCACCATGCTCCGCCTCC  >  NC_002947/5137675‑5137952

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