Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,633,810	A→G	Y5C (TAC→TGC)	ytjC →	putative phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,633,810	0	A	G	96.4%	90.3 / ‑7.0	28	Y5C (TAC→TGC)	ytjC	putative phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (16/11);  minor base C (0/1);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCTTGGCAT  >  NC_000913/4633690‑4633939                                                                                                                         |                                                                                                                                  tGCAGAATGGCCTGACTTTTAGCGAGATGGGTGGTCGCACAAACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCa                                                                                                                       <  2:167629/133‑1 (MQ=255)  gCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCAc                                                                                                                      >  2:410847/1‑133 (MQ=255)    aGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGG‑‑AAAAAGCATGTTACAAGTAGGCATAGTCCCCCACGGTg                                                                                                                  >  1:207577/1‑133 (MQ=255)                       cGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGcc                                                                                                 <  1:66510/133‑1 (MQ=255)                       cGGGATTGGTAGTCGGACAGACAAATTGATGCATAATCAGCTTTACTCAGAAAATTAACGTTACAGCAGTATACGG‑‑AAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGa                                                                                               <  2:207577/133‑1 (MQ=255)                         ggATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGa                                                                                               <  1:410847/133‑1 (MQ=255)                          gATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAg                                                                                              <  2:491375/133‑1 (MQ=255)                              ggtggtCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGAc                                                                                          >  1:467666/1‑133 (MQ=255)                               gtggtCGCACTGACATCTTGGTGCTTCATCAGTATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACg                                                                                         <  2:233934/133‑1 (MQ=255)                                                 tGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGTGACATAGTCAGGGCCAGTCTGa                                                                       >  1:163222/1‑133 (MQ=255)                                                         aaTCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGc                                                               >  2:46057/1‑133 (MQ=255)                                                           tCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTg                                                             <  1:85642/133‑1 (MQ=255)                                                                   aCTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCTGCTGACCGCCaa                                                     >  2:114816/1‑133 (MQ=255)                                                                    cTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCaaa                                                    >  2:257050/1‑133 (MQ=255)                                                                       aGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGt                                                 >  2:54568/1‑133 (MQ=255)                                                                              tAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAAc                                          >  1:64459/1‑133 (MQ=255)                                                                              tAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAAc                                          >  1:488465/1‑133 (MQ=255)                                                                                   ttACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCg                                     >  2:362563/1‑133 (MQ=255)                                                                                            tataCGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTg                            >  1:66947/1‑133 (MQ=255)                                                                                              taCGGAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGc                          >  1:45414/1‑133 (MQ=255)                                                                                                 ggAAAAAAAGAATGTTACAGGTATCCCTTGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAAc                       <  2:163222/133‑1 (MQ=255)                                                                                                  gAAAAAAAGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAAcc                      >  1:27365/1‑133 (MQ=255)                                                                                                   aaaaaaaGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAAccc                     >  2:277424/1‑133 (MQ=255)                                                                                                      aaaaGCATGTTACAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTg                  >  2:10549/1‑133 (MQ=255)                                                                                                                 aCAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCtt       <  2:64459/133‑1 (MQ=255)                                                                                                                 aCAGGTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGCGCCAAAGAGCtt       <  2:488465/133‑1 (MQ=255)                                                                                                                    ggTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCTGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCTTGGc    <  1:114816/133‑1 (MQ=255)                                                                                                                     gTATGCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCTTGGCa   <  2:27365/133‑1 (MQ=255)                                                                                                                      tatgCCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCTTGGCAt  <  1:54568/133‑1 (MQ=255)                                                                                                                         |                                                                                                                                  TGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCTTGGCAT  >  NC_000913/4633690‑4633939

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