Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,307,460	G→A	9.3%	intergenic (‑161/+25)	recC ← / ← PP_5SG	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit/5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,307,460	0	G	A	9.3%	109.8 / 3.0	43	intergenic (‑161/+25)	recC/PP_5SG	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (15/24);  new base A (2/2);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

GAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGG  >  NC_002947/5307313‑5307569                                                                                                                                                    |                                                                                                              gAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgcgtgcg                                                                                                             >  1:541463/1‑150 (MQ=255)        tGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCGAAAGCAGGGATTCGCGCAGGTGGAGGTTTTTGCCCACAAAGACAAAACTCCTACCTGCGTGCGCAGATAg                                                                                                      >  1:443567/1‑150 (MQ=38)        ctaTGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAg                                                                                                      <  1:219705/148‑1 (MQ=255)           tGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAAGGACATCAGGCTGCGGGAGCGGCCGTGCGCGCTCCAGCAGCGATTCGCGCAGGTCGAGTTTTTTGCGCACAAAGACACAACCCCTATCCGCGCACGCAGGTAggag                                                                                                   >  2:219705/1‑148 (MQ=255)              cTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAgggg                                                                                                   >  1:482065/1‑147 (MQ=255)              cTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAgggg                                                                                                   <  2:482065/147‑1 (MQ=255)              atGAGAGGGGACCGATATCACGCGTACATGAAACGGCCGAATGCCAGCAGGCTGCGGGAGCGGCCGTGCCCGCTACCGCAGCGATTCGAGCAGGTCGAGTTTTTTACCCACAGAAACAAAACCCCTACCTGCTTGCGCAGATAGGGGttt                                                                                                <  2:1252094/149‑1 (MQ=255)                     gggACCGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaaa                                                                                         <  1:405064/150‑1 (MQ=255)                          cGATCTCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACACAGACCAACCCCATACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAAt                                                                                    <  1:35630/150‑1 (MQ=37)                          cGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGAGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAAt                                                                                    <  1:1296133/150‑1 (MQ=255)                                cACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGAc                                                                              <  2:541463/150‑1 (MQ=255)                                 aCGCGTACATGAAACGGCCGAATGCCATAAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCCCAGGTCGAGTTTTTTGCCCACAAAGACACAACCCCTCCCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATccgatag                                                                             >  2:35630/1‑144 (MQ=37)                                         aTGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTa                                                                     >  1:169166/1‑150 (MQ=255)                                         aTGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTa                                                                     >  1:702744/1‑150 (MQ=255)                                             aaCGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAggg                                                                                                    <  1:302981/115‑1 (MQ=38)                                             aaCGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAggg                                                                                                    >  2:302981/1‑115 (MQ=38)                                                   cGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCt                                                                                  <  2:1078286/127‑1 (MQ=255)                                                   cGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCt                                                                                  >  1:1078286/1‑127 (MQ=255)                                                    gAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                                                          <  2:160689/150‑1 (MQ=255)                                                     aaTGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATggg                                                         <  2:417555/150‑1 (MQ=255)                                                      aTGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGTGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgggg                                                        <  1:917768/150‑1 (MQ=255)                                                          cATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAAt                                                                                    <  1:331809/118‑1 (MQ=255)                                                          cATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAAt                                                                                    >  2:331809/1‑118 (MQ=255)                                                          cATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCTCATGGGGaaag                                                    >  2:324720/1‑149 (MQ=255)                                                          cATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                                                          <  1:95081/144‑1 (MQ=255)                                                          cATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                                                          >  2:95081/1‑144 (MQ=255)                                                           aTCAGGCTGCGGTAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGa                                                                               <  1:440416/122‑1 (MQ=255)                                                           aTCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGa                                                                               >  2:440416/1‑122 (MQ=255)                                                                  tGCGGGAGCGGGCACGGCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACt                                            >  1:354674/1‑150 (MQ=17)                                                                            gCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAt                                                                                        >  2:896691/1‑96 (MQ=255)                                                                            gCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAt                                                                                        <  1:896691/96‑1 (MQ=255)                                                                             ccGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGa                                 <  1:612801/150‑1 (MQ=255)                                                                              cGTGCCCGCCCCCGCAGCGATTCGCGCAGGTCCAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTCGCGAAATGAATCTTGACGACGACCTACTCTCACATGCGGAAACCCCACACTACCATCGGCGAt                                <  1:324720/150‑1 (MQ=14)                                                                               gTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATg                               <  2:986471/150‑1 (MQ=255)                                                                               gTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATg                               <  2:718643/150‑1 (MQ=255)                                                                               gTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATg                               <  2:169166/150‑1 (MQ=255)                                                                                     gCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGCCGATggtagtg                         <  1:713767/150‑7 (MQ=255)                                                                                         ccGCAGCGATTCGCGCAGGTCGAGTACTGTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCa                     <  1:1236657/150‑1 (MQ=255)                                                                                           gCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGAc                                                                              <  1:1034635/91‑1 (MQ=255)                                                                                           gCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGAc                                                                              >  2:1034635/1‑91 (MQ=255)                                                                                           gCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCac                    >  1:1144605/1‑149 (MQ=255)                                                                                           gCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCac                    <  2:1144605/149‑1 (MQ=255)                                                                                                gATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGa              <  2:702744/150‑1 (MQ=40)                                                                                                            tCGACTTTTTCGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGCCCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCCTCGTTTTACTCCTGCGTTCgggatggg  >  2:812395/1‑150 (MQ=11)                                                                                                                                                    |                                                                                                              GAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGG  >  NC_002947/5307313‑5307569

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