Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,174,765	A→C	10.9%	intergenic (‑71/+48)	glmU ← / ← atpC	bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/ATP synthase subunit epsilon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,174,765	0	A	C	10.9%	130.0 / 9.2	55	intergenic (‑71/+48)	glmU/atpC	bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/ATP synthase subunit epsilon
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base C (3/3);  total (28/27)
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.52e-01

GTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGATT  >  NC_002947/6174620‑6174914                                                                                                                                                  |                                                                                                                                                      gTGCAACACTTTGGGCATCGCCGAGCGCATGGGGGTGGCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGATACCAATCCTGGGCGTTCAGCGTACGCCGCCGGGACAAAACGGATACAGTAAAAGGGTAGCCGAGGCt                                                                                                                                                   >  2:87779/1‑150 (MQ=255)       cacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCACGTCAGACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCtt                                                                                                                                             <  1:87779/150‑1 (MQ=255)        acTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCttt                                                                                                                                            <  2:332327/150‑1 (MQ=255)              ggCAGCGCCGAGCGCATGCGGGTGTCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCa                                                                                                                                      <  1:1126681/150‑1 (MQ=255)               gCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCaa                                                                                                                                     >  1:600443/1‑150 (MQ=255)                   cgcCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACGAATCCTGGGCGGTCAGTGAAAGAACGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGc                                                                                                                                 >  1:185504/1‑150 (MQ=255)                           gcATGCGGGTGCCTTGGCCGGCGGGGAGAATAACGATAGCGAGTGACAATACTTGAGACCAATCCTGGGACGTCAGTTAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAg                                                                                                                         <  2:338805/150‑1 (MQ=255)                               gCGGGTGCCTTGACCGACGGCGAGAATAACGATATAGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCTGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTAt                                                                                                                     <  2:185504/150‑1 (MQ=255)                                   gTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACTCCTTGACGCAATCGCAACGGAAGTTATCAGc                                                                                                                 >  2:373262/1‑150 (MQ=255)                                       cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCggg                                                                                                             >  2:948062/1‑150 (MQ=255)                                         tGGCCGGGGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGc                                                                                                           <  2:789011/150‑1 (MQ=255)                                                  gcgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGAGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACtt                                                                                                  >  2:445467/1‑150 (MQ=255)                                                       aaTAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCtt                                                                                             >  1:565531/1‑150 (MQ=255)                                                         tAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTg                                                                                           >  1:1103967/1‑150 (MQ=255)                                                          aaCGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCtttt                                                                                                                                           >  1:945863/1‑101 (MQ=255)                                                          aaCGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCtttt                                                                                                                                           <  2:945863/101‑1 (MQ=255)                                                          aaCGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgc                                                                                          <  1:553102/150‑1 (MQ=255)                                                                 tCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGcc                                                                                                                <  2:67298/121‑1 (MQ=255)                                                                 tCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGcc                                                                                                                >  1:67298/1‑121 (MQ=255)                                                                  cGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGc                                                                                  <  2:405305/150‑1 (MQ=255)                                                                   gAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCttt                                                                                                                                            <  2:347714/91‑1 (MQ=255)                                                                   gAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCttt                                                                                                                                            >  1:347716/1‑91 (MQ=255)                                                                   gAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCttt                                                                                                                                            >  1:347714/1‑91 (MQ=255)                                                                   gAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCttt                                                                                                                                            <  2:347716/91‑1 (MQ=255)                                                                     gTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGa                                                                                                                                                       >  1:890506/1‑78 (MQ=255)                                                                     gTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGa                                                                                                                                                       <  2:890506/78‑1 (MQ=255)                                                                       gACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCttt                                                                                                                                            <  1:1398128/87‑1 (MQ=255)                                                                       gACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCttt                                                                                                                                            >  2:1398128/1‑87 (MQ=255)                                                                          aTTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCtt                                                                                                                                             >  2:906717/1‑83 (MQ=255)                                                                          aTTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCtt                                                                                                                                             <  1:906717/83‑1 (MQ=255)                                                                          aTTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGt                                                                          >  1:1131275/1‑150 (MQ=255)                                                                                 gCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGc                                                                   >  2:694989/1‑150 (MQ=255)                                                                                   tACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCATCTCCCATTTGCTGGACGGTACGCAGCTg                                                                 <  1:1108087/150‑1 (MQ=255)                                                                                    aCCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCtt                                                                                                                                             <  1:4827/73‑1 (MQ=255)                                                                                    aCCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCtt                                                                                                                                             >  2:4827/1‑73 (MQ=255)                                                                                       aaTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCt                                                             <  1:730418/150‑1 (MQ=255)                                                                                          ccTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAGAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCg                                                          <  1:373262/150‑1 (MQ=255)                                                                                             gggCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGcc                                                       >  1:1015836/1‑150 (MQ=255)                                                                                                  gTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGc                                                  <  1:445467/150‑1 (MQ=255)                                                                                                     aGTGACCGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAg                                               >  2:1274537/1‑150 (MQ=255)                                                                                                       tGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGAc                                             <  2:600443/150‑1 (MQ=255)                                                                                                       tGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGAc                                             >  2:1076408/1‑150 (MQ=255)                                                                                                            gACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCNAGGCTACCCNTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcg                                        >  1:629208/1‑150 (MQ=255)                                                                                                            gACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcg                                        <  2:1103967/150‑1 (MQ=255)                                                                                                                            agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGt                        >  1:1163639/1‑150 (MQ=255)                                                                                                                              aaaaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCg                      >  1:503961/1‑150 (MQ=255)                                                                                                                                aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgc                                         >  2:340115/1‑129 (MQ=255)                                                                                                                                aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgc                                         <  1:340115/129‑1 (MQ=255)                                                                                                                                aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGaa                    >  2:121862/1‑150 (MQ=255)                                                                                                                                     aaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgc                                         >  2:280233/1‑124 (MQ=255)                                                                                                                                     aaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgc                                         <  1:280233/124‑1 (MQ=255)                                                                                                                                      aaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcgc              <  1:694989/150‑1 (MQ=255)                                                                                                                                          ggTAGCCTCGGCTGCCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCa                                                <  1:817772/112‑1 (MQ=255)                                                                                                                                          ggTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCa                                                >  2:1082126/1‑112 (MQ=255)                                                                                                                                          ggTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCa                                                >  2:817772/1‑112 (MQ=255)                                                                                                                                          ggTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCa                                                <  1:1082126/112‑1 (MQ=255)                                                                                                                                                 gAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGAt   <  2:676582/150‑1 (MQ=255)                                                                                                                                                  aGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGAtt  <  2:1131275/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGATT  >  NC_002947/6174620‑6174914

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