Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,130,198	A→T	10.5%	intergenic (‑247/+57)	gcvH‑I ← / ← PP_0990	glycine cleavage system protein H/DoxX family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,130,198	0	A	T	10.5%	161.9 / 11.1	57	intergenic (‑247/+57)	gcvH‑I/PP_0990	glycine cleavage system protein H/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base T (2/4);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAG  >  NC_002947/1130053‑1130338                                                                                                                                                   |                                                                                                                                             ttatCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATAACTTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTa                                                                                                                                          <  1:656176/150‑1 (MQ=255)         aGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGAAAAc                                                                                                                                  <  2:706333/150‑1 (MQ=255)           aTTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGAAAACGa                                                                                                                                <  1:573920/150‑1 (MQ=255)           aTTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGa                                                                                                                                >  2:136998/1‑150 (MQ=255)                gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATa                                                                                                                              >  2:412959/1‑147 (MQ=255)                gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATa                                                                                                                              <  1:412959/147‑1 (MQ=255)                          ggAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTg                                                                                                                 >  2:384120/1‑150 (MQ=255)                               aGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTg                                                                                                            <  1:384120/150‑1 (MQ=255)                                   gTATTTATTGTGTGGAATGATATCGTTACAGGGGCCTCGGGGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTGTTCGCTACGATGTAGCGAAAACGATACATGTGGTGCCGCGACTgggta                                                                                                        <  2:572870/150‑4 (MQ=255)                                        tattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGc                                                                                                   >  1:559708/1‑150 (MQ=255)                                             tgtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATg                                                                                              >  2:103754/1‑150 (MQ=255)                                                ggAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCg                                                                                                         <  1:102871/136‑1 (MQ=255)                                                ggAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCg                                                                                                         >  2:102871/1‑136 (MQ=255)                                                ggAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGcc                                                                                           >  2:945399/1‑150 (MQ=255)                                                  aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGccag                                                                                         >  1:1238326/1‑148 (MQ=255)                                                       tatCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTc                                                                                    <  1:285084/150‑1 (MQ=255)                                                        atCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc                                                                                                          <  2:852615/127‑1 (MQ=255)                                                        atCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc                                                                                                          >  1:852615/1‑127 (MQ=255)                                                        atCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCa                                                                                   >  1:375507/1‑150 (MQ=255)                                                          cGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc                                                                                 <  1:136998/150‑1 (MQ=255)                                                            ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtctc                                                                               <  1:103754/150‑1 (MQ=255)                                                               cAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAACACTGTAGGCGGGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACt                                                                            <  1:732400/150‑1 (MQ=255)                                                                  gTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc                                                                                 >  1:388503/1‑142 (MQ=255)                                                                  gTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc                                                                                 <  2:388503/142‑1 (MQ=255)                                                                     ccTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagc                                                                      <  2:559708/150‑1 (MQ=255)                                                                      ctctGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCATGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagca                                                                     >  2:310345/1‑150 (MQ=255)                                                                         tGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAgat                                                                  <  2:1092587/150‑1 (MQ=255)                                                                           agagGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAgatga                                                                >  1:1011023/1‑150 (MQ=255)                                                                               gAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCa                                                            >  2:344754/1‑150 (MQ=255)                                                                               cacaaGGGTTTGATGC‑TCAAACCCTTTTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCa                                                            >  1:108123/4‑150 (MQ=255)                                                                                 aaaGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGc                                                          <  1:344754/150‑1 (MQ=255)                                                                                              cATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCg                                              >  2:287395/1‑150 (MQ=255)                                                                                                    aCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagc                                                                      <  2:107141/120‑1 (MQ=255)                                                                                                    aCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagc                                                                      >  1:107141/1‑120 (MQ=255)                                                                                                      ccTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAg                                                                                                       >  2:1064942/1‑85 (MQ=255)                                                                                                      ccTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAg                                                                                                       <  1:1064942/85‑1 (MQ=255)                                                                                                       cTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCgg                                     <  1:620366/150‑1 (MQ=255)                                                                                                        ttGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCggg                                    >  1:702478/1‑150 (MQ=255)                                                                                                          gtgCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcgc                                                  >  1:233324/1‑134 (MQ=255)                                                                                                          gtgCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcgc                                                  <  2:233324/134‑1 (MQ=255)                                                                                                          gtgCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATACAGCGAAAACGCGCCGGCGCCTCGGGCa                                  >  1:133656/1‑150 (MQ=255)                                                                                                               ccGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAg                             <  2:1011023/150‑1 (MQ=255)                                                                                                                  aCGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGac                          >  2:158124/1‑149 (MQ=255)                                                                                                                          taggCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAaga                  >  2:879742/2‑150 (MQ=255)                                                                                                                           aggCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAagag                 <  2:375507/150‑1 (MQ=255)                                                                                                                                gtgtATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTCCGGTCATCTCACTTTCAGCAGATGTTGTAGCGAAAACGCGCCTGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTg            <  1:158124/150‑1 (MQ=255)                                                                                                                                       gTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTg            <  2:232639/143‑1 (MQ=255)                                                                                                                                       gTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTg            >  1:232639/1‑143 (MQ=255)                                                                                                                                        ttttCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGagca                        >  1:811635/1‑130 (MQ=255)                                                                                                                                        ttttCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGagca                        <  2:811635/130‑1 (MQ=255)                                                                                                                                        ttttCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGag                          >  1:580636/1‑128 (MQ=255)                                                                                                                                        ttttCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGag                          <  2:580636/128‑1 (MQ=255)                                                                                                                                        ttttCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtc    >  2:495295/1‑150 (MQ=255)                                                                                                                                        ttttCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtc    >  2:543637/1‑150 (MQ=255)                                                                                                                                        ttttCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtc    >  1:629717/1‑150 (MQ=255)                                                                                                                                        ttttCGCTACGACGTAGCGAAAACGATACATGTGGTACCTTGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTGATCGGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCCGGAGCAGCAAGAGCCCTGCCcaggtc    >  1:103391/1‑150 (MQ=255)                                                                                                                                          ttCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtcag  <  1:879742/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             TTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGA‑GCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAG  >  NC_002947/1130053‑1130338

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: ‑