Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,307,019	A→G	7.5%	V78A (GTG→GCG)	recC ←	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,307,019	0	A	G	7.5%	148.5 / 2.3	53	V78A (GTG→GCG)	recC	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base G (2/2);  total (25/28)
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

GCTCGATGAGCGCTGGCAACAGGCGCATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCAGAGGCGCCAGCGGGTAGCGACGCATC  >  NC_002947/5306885‑5307165                                                                                                                                       |                                                                                                                                                    gcTCGATGAGCGCTGGCAACAGGCGCATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAgcg                                                                                                                                       <  2:330194/149‑1 (MQ=255)   tCGATGAGCGCTGGCAACAGGCGCATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAgcgg                                                                                                                                      >  5:154529/1‑148 (MQ=255)     gATGAGCGCTGGCAACAGGCGCATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAgcgg                                                                                                                                      <  6:154529/146‑1 (MQ=255)         aGCGCTGGCAACAGGCGCATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGc                                                                                                                                <  2:44843/148‑1 (MQ=255)                     aGGCGCATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTg                                                                                                                   <  3:304831/149‑1 (MQ=255)                      ggCGCATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTg‑g                                                                                                                 <  1:177588/149‑1 (MQ=255)                        cgcATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTc                                                                                                               <  4:271084/149‑1 (MQ=255)                         gcATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCa                                                                                                              >  7:303880/1‑149 (MQ=255)                           aTCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATc                                                                                                            >  3:148929/1‑149 (MQ=255)                           aTCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATc                                                                                                            >  2:85600/1‑149 (MQ=255)                                gcgcCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAGGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATc                                                                                                            >  2:62513/1‑144 (MQ=255)                                gcgcCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAACGCTACCGGGCAGTTGCACAACGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATc                                                                                                            <  1:62513/144‑1 (MQ=255)                                   ccAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCg                                                                                                    <  1:279226/149‑1 (MQ=255)                                       gTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATc                                                                                                <  6:81648/149‑1 (MQ=255)                                         cAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCtt                                                                                                                    <  8:193628/128‑1 (MQ=255)                                         cAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCtt                                                                                                                    >  7:193628/1‑128 (MQ=255)                                           ggggTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCa                                                                                            <  1:85600/149‑1 (MQ=255)                                            gggTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAg                                                                                           >  4:74081/1‑149 (MQ=255)                                              gTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGc                                                                                          >  6:174770/1‑148 (MQ=255)                                                 ccTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAgg                                                                                      >  1:172262/1‑149 (MQ=255)                                                           cAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGGTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGa                                                                            <  3:208955/149‑1 (MQ=255)                                                              agaCACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCTGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGAt                                                                                                 <  4:301890/125‑1 (MQ=255)                                                              agaCACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGAt                                                                                                 >  3:301890/1‑125 (MQ=255)                                                                 cacTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCAtt                                                                      <  5:174770/149‑1 (MQ=255)                                                                  acTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCGCATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGc                                                                                                     >  4:1631/1‑117 (MQ=255)                                                                  acTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCGCATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGc                                                                                                     <  3:1631/117‑1 (MQ=255)                                                                    ttCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGcac                                                                                                                                                    <  8:314665/69‑1 (MQ=255)                                                                        gggATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAAt                                                               <  7:135442/149‑1 (MQ=255)                                                                        gggATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAAt                                                               >  5:136402/1‑149 (MQ=255)                                                                          gATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGc                                                             >  3:277269/1‑149 (MQ=255)                                                                                gCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGc                                                       >  5:270899/1‑149 (MQ=255)                                                                                      gCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCa                                                 >  8:186913/1‑149 (MQ=255)                                                                                        cAAGCACGCTGCGGTAGAGCTGCCACATAATGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCAc                                                >  1:30686/1‑148 (MQ=255)                                                                                         aaGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCa                                              <  4:148929/149‑1 (MQ=255)                                                                                           gCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCAtt                                                         >  8:200902/1‑136 (MQ=255)                                                                                           gCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCAtt                                                         <  7:200902/136‑1 (MQ=255)                                                                                            cACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAggg                                           <  7:186913/149‑1 (MQ=255)                                                                                                  gCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAAttt                                     >  1:334797/1‑149 (MQ=255)                                                                                                    ggTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTc                                    <  8:35217/148‑1 (MQ=255)                                                                                                    ggTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTc                                    >  7:35217/1‑148 (MQ=255)                                                                                                      tAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGtt                                 <  6:136402/149‑1 (MQ=255)                                                                                                        gagCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGttt                                <  4:277269/148‑1 (MQ=255)                                                                                                              ccACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCagag                         >  5:71062/1‑149 (MQ=255)                                                                                                               cacaTAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTGAGTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCagag                         >  1:274332/1‑149 (MQ=255)                                                                                                                acaTAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAAt                                                               <  4:102609/109‑1 (MQ=255)                                                                                                                acaTAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAAt                                                               >  3:102609/1‑109 (MQ=255)                                                                                                                           ccGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCAGAGGCGCCAGCgg               <  2:274332/146‑1 (MQ=255)                                                                                                                           ccGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCAGAGGCGCCAGCGGGt             <  6:270899/148‑1 (MQ=255)                                                                                                                           ccGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCAGAGGCGCCAGCGGGTa            <  2:30686/149‑1 (MQ=255)                                                                                                                            cGGGCAGTTGCGCATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑Gt                                                                                                                >  1:284881/1‑48 (MQ=255)                                                                                                                            cGGGCAGTTGCGCATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑Gt                                                                                                                <  2:284881/48‑1 (MQ=255)                                                                                                                            cGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGc                                                                          >  2:32176/1‑86 (MQ=255)                                                                                                                             gggCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTc                              <  7:258645/129‑1 (MQ=255)                                                                                                                             gggCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTc                              >  8:258645/1‑129 (MQ=255)                                                                                                                                      cacaTCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCAGAGGCGCCAGCGGGTAGCGACGCATc  >  6:299150/1‑148 (MQ=255)                                                                                                                                      cacaTCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCAGAGGCGCCAGCGGGTAGCGACGCATc  <  2:172262/148‑1 (MQ=255)                                                                                                                                       acaTCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCAGAGGCGCCAGCGGGTAGCGACGc     >  3:258219/1‑144 (MQ=255)                                                                                                                                       |                                                                                                                                                    GCTCGATGAGCGCTGGCAACAGGCGCATCAGGCGCCAGGTCAGGGGTGCCTTGTCGAGCAGAGACACTTCGGGGATTTCGCTACGGCCAAGCACGCTGCGGTAGAGCTGCCACATAAAGCTACCGGGCAGTTGCACATCGATTGCAGCGGCGATGCCGCAGCCGCCTTG‑GTCATCTTCCAGCGGATCTTCAGCCAGGGCCAGCTTGAGCCATTGGGCAATGCCATTGCTTTGCACCAGGGCAATTTCGTTTTCCAGAGGCGCCAGCGGGTAGCGACGCATC  >  NC_002947/5306885‑5307165

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑