Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,925,007	C→T	intergenic (‑38/‑101)	yobA ← / → holE	copper‑binding protein/DNA polymerase III subunit theta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,925,007	0	C	T	97.0%	94.7 / ‑4.3	33	intergenic (‑38/‑101)	yobA/holE	copper‑binding protein/DNA polymerase III subunit theta
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (11/21);  total (11/22)
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.54e-01

GGCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGCATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACA  >  NC_000913/1924767‑1925247                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 ggCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTg                                                                                                                                                                                                                                         <  1:444767/250‑1 (MQ=255)          aTTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCt                                                                                                                                                                                                                               <  2:433464/251‑1 (MQ=255)             tttGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGa                                                                                                                                                                                                                             <  2:93832/250‑1 (MQ=255)              ttGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGTAGTGGTCATGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTGTTGTAATAACTTTTTTACAGAGTATATCCTTGTCTAATGTCTGAg                                                                                                                                                                                                                            <  1:421098/250‑1 (MQ=255)                                            ccGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCt                                                                                                                                                                                             <  1:312581/251‑1 (MQ=255)                                                gaAGTTTAAAGTGCTTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGCTGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGCGGATCATCAATTCCGGCTTGCCATCCTGGCt                                                                                                                                                                                         >  2:628836/1‑251 (MQ=255)                                                gaAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGAGCGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGCATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCt                                                                                                                                                                                         <  1:232541/251‑1 (MQ=255)                                                     ttAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACAAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGTGCGTGCATTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTGACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCAct                                                                                                                                                                                     <  1:588140/250‑1 (MQ=255)                                                      tAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCActc                                                                                                                                                                                    <  1:129983/250‑1 (MQ=255)                                                                     gggTGCAGCTGTCACTTGCGCGTTTGCGGCAGGATACGGATGCTTTAAATGCGCATGTGCCCAGACGGAAGGTTTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGccc                                                                                                                                                                    <  2:341738/250‑1 (MQ=255)                                                                                     tGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATgagga                                                                                                                                                      <  1:360407/249‑1 (MQ=255)                                                                                      gcgcGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGAt                                                                                                                                                   <  2:16804/251‑1 (MQ=255)                                                                                               ccGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGc                                                                                                                                          >  2:161555/1‑251 (MQ=255)                                                                                                                  tAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATTAGGAGAGTAAGAATGCTGAATAATCTGGCTAAAc                                                                                                                        <  2:253353/250‑1 (MQ=255)                                                                                                                  tAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACt                                                                                                                       >  1:158850/1‑251 (MQ=255)                                                                                                                   aaaTGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACt                                                                                                                       <  2:166223/250‑1 (MQ=255)                                                                                                                     aTGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTgg                                                                                                                     >  1:19564/1‑250 (MQ=255)                                                                                                                         gcATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCa                                                                                                                 >  1:587856/1‑250 (MQ=255)                                                                                                                          cATGTGCCCAGACGGAAGGTGTCACCAACGCAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCaa                                                                                                                >  2:538163/1‑250 (MQ=255)                                                                                                                             gtgCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAAc                                                                                                              >  2:126193/1‑249 (MQ=255)                                                                                                                                                 aCCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTGTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGATGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTc                                                                                        <  1:62026/251‑1 (MQ=255)                                                                                                                                                          gttgTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTTAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTGTTGCTCGCGAATGATGAGGAGATTAAGAATGCTTAAGAATCTGGCTAAACGGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTggcg                                                                               <  2:593346/248‑1 (MQ=255)                                                                                                                                                          gtgttCAGTATTGCCAGCGCGTAGCGAAGTGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATATCATTGTCTAATGTCTGAGTCGATGATCATCAATTCCGGCTTGCCATCCGGGCTCACTCGTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTggcg                                                                               <  1:402637/247‑1 (MQ=255)                                                                                                                                                            ggTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTggcggc                                                                             >  1:149960/1‑251 (MQ=255)                                                                                                                                                                    ttGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCgggg                                                                     >  2:149760/1‑251 (MQ=255)                                                                                                                                                                       ccAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATGTGGCGGCGGCCGGGGTgg                                                                  <  2:569506/251‑1 (MQ=255)                                                                                                                                                                       ccAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCGCACGCTTAGTAACTTTTGCCCGCGAATGGTGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCTGCCGGGGTgg                                                                  <  1:415153/251‑1 (MQ=255)                                                                                                                                                                                             ggCAGTTGAAGCCATATTATCTAATCCTTTTTTTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAATTTTTGCCCGCGAATGATGAGGAGATTAAGAAGGCTGAAGAATCTGGCTAAACTGGATCGAACAGAAATGGATAAAGTGAATGTCGATTGGGCGGCGGCTGGGTGGGCGTTGAAAGAActcgccaatat                                            <  2:329303/250‑11 (MQ=255)                                                                                                                                                                                                                            tGTAATAACTTTTTTACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGa              <  2:19564/250‑1 (MQ=255)                                                                                                                                                                                                                                      tttttACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTgg                                                                                                                     >  2:294245/1‑137 (MQ=255)                                                                                                                                                                                                                                      tttttACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTgg                                                                                                                     <  1:294245/137‑1 (MQ=255)                                                                                                                                                                                                                                       ttttACAGAGTATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAAc   <  2:158850/250‑1 (MQ=255)                                                                                                                                                                                                                                       ttttACAGAGTATAACATTGTCTAATGTCTGAGTCGAGGATCATAACTTCCGGCTTGCCCTCCTGGCCAAATCTTAGTAAATTTTGCCCGCGAATGCTGAGGAGATTAAGAATGCTGAAGAATACGGCTAAACTGGATCAAACAGAACTGGATAAAGTGAATGTCGATTTGGAGGCGGCCGGGGTGGCGTTTAAAGAAAGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGACCAGCCTGAAcc  >  1:58029/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 GGCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGCATAACCTTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACA  >  NC_000913/1924767‑1925247

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