Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,634,054	T→C	G86G (GGT→GGC)	ytjC →	putative phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,634,054	0	T	C	97.5%	144.3 / ‑5.0	40	G86G (GGT→GGC)	ytjC	putative phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (17/22);  minor base G (0/1);  total (17/23)
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

CCCGTGCCAAAGAGCTTGGCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGG  >  NC_000913/4633918‑4634190                                                                                                                                         |                                                                                                                                         cccGTGCCAAAGAGCTTGGCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGcgtg                                                                                                                                       <  2:96123/140‑1 (MQ=255) cccGTGCCAAAGAGCTTGGCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGcgtg                                                                                                                                       <  1:500350/140‑1 (MQ=255)   cGTGCCAAAGAGCTTGGCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGGCTGCGTGAATTAAACATGGGGGTGCt                                                                                                                                     <  2:125292/140‑1 (MQ=255)    gTGCCAAAGAGCTTGGCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTg                                                                                                                                    >  1:973686/1‑140 (MQ=255)       ccAAAGAGCTTGGCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGGCATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGaa                                                                                                                                 <  2:1049774/140‑1 (MQ=255)             agCTTGGCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGa                                                                                                                           <  1:365597/140‑1 (MQ=255)               cTTGGCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACa                                                                                                                         >  1:1186717/1‑140 (MQ=255)                  ggCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACatat                                                                                                                      >  1:98026/1‑140 (MQ=255)                      ttACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGAt                                                                                                                  <  1:212392/140‑1 (MQ=255)                         cGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATtct                                                                                                               <  2:98026/140‑1 (MQ=255)                              atcatcAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGAc                                                                                                          <  2:281238/140‑1 (MQ=255)                              atcatcAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGAc                                                                                                          <  1:339294/140‑1 (MQ=255)                              atcatcAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGAc                                                                                                          <  2:1072071/140‑1 (MQ=255)                                   caGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCgaag                                                                                                     >  1:1057120/1‑140 (MQ=255)                                         cGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGaagaagaag                                                                                               >  1:531647/1‑140 (MQ=255)                                                      aCCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCg                                                                                  >  1:138119/1‑140 (MQ=255)                                                      aCCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCg                                                                                  <  2:1186717/140‑1 (MQ=255)                                                           gCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGc                                                                             >  1:191344/1‑140 (MQ=255)                                                             gTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTg                                                                           >  1:879841/1‑140 (MQ=255)                                                             gTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTg                                                                           <  2:651294/140‑1 (MQ=255)                                                                  gcggAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCaa                                                                      <  2:973686/140‑1 (MQ=255)                                                                       aaTCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCa                                                                 <  2:191344/140‑1 (MQ=255)                                                                               cccAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGAc                                                         <  2:531647/140‑1 (MQ=255)                                                                                 cAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACgg                                                       <  2:58441/140‑1 (MQ=255)                                                                                         cGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTc                                               >  1:526583/1‑140 (MQ=255)                                                                                              gtgACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGc                                                     >  1:482709/1‑129 (MQ=255)                                                                                              gtgACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGc                                                     <  2:482709/129‑1 (MQ=255)                                                                                              gtgACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGaa                                          >  1:895180/1‑140 (MQ=255)                                                                                                            gATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCa                            >  2:826911/1‑140 (MQ=255)                                                                                                            gATTCTCGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAAAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCa                            <  2:1057120/140‑1 (MQ=255)                                                                                                                 tcGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTc                                  >  1:555597/1‑129 (MQ=255)                                                                                                                 tcGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTc                                  <  2:555597/129‑1 (MQ=255)                                                                                                                  cGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTc                                  >  1:533474/1‑128 (MQ=255)                                                                                                                  cGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTc                                  <  2:533474/128‑1 (MQ=255)                                                                                                                  cGCCTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCt                      >  1:806466/1‑140 (MQ=255)                                                                                                                    ccTGCGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCa                    <  2:895180/140‑1 (MQ=255)                                                                                                                        cGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGc                             <  2:160637/127‑1 (MQ=255)                                                                                                                        cGTGAATTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGc                             >  1:160637/1‑127 (MQ=255)                                                                                                                             aTTAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATcgcg           >  1:469807/1‑140 (MQ=255)                                                                                                                              ttAAACATGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGt          <  2:840119/140‑1 (MQ=255)                                                                                                                                    aTGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGc    <  2:451635/140‑1 (MQ=255)                                                                                                                                    aTGGGCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGc    <  1:826911/140‑1 (MQ=255)                                                                                                                                      gggCGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCgg  <  2:806466/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         CCCGTGCCAAAGAGCTTGGCATTACGCATATCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGG  >  NC_000913/4633918‑4634190

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