Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,268,809	C→T	*68* (TAG→TAA)	yjbS ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,268,809	0	C	T	96.7%	90.3 / ‑6.2	30	*68* (TAG→TAA)	yjbS	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/10);  minor base A (0/1);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATAT  >  NC_000913/4268674‑4268939                                                                                                                                        |                                                                                                                                   tAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGttaattaa                                                                                                                                 <  1:916010/139‑1 (MQ=255)                         cAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACttt                                                                                                         >  2:799996/1‑139 (MQ=255)                         cAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACttt                                                                                                         <  1:799996/139‑1 (MQ=255)                            tatCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACt                                                                                                      >  2:835462/1‑139 (MQ=255)                            tatCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACt                                                                                                      <  1:372491/139‑1 (MQ=255)                                   tAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATa                                                                                               <  1:835462/139‑1 (MQ=255)                                       aaGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAAc                                                                                           >  2:464127/1‑139 (MQ=255)                                          gTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTg                                                                                        >  1:746156/1‑139 (MQ=255)                                             cTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATa                                                                                     >  2:503292/1‑139 (MQ=255)                                                     agaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGttaattaa                                                                                                                                 >  2:813188/1‑87 (MQ=255)                                                     agaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGttaattaa                                                                                                                                 <  1:813188/87‑1 (MQ=255)                                                     agaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTtcttct                                                                             >  2:829020/1‑139 (MQ=255)                                                     agaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTtcttct                                                                             >  1:739430/1‑139 (MQ=255)                                                      tagagaccgccTTTAACCTCCTCCGCAAAAAAGGAAAGCGTGGAAATATACATAAAGTCTTTCTTACCCCTTATTTTGTTAAATAAATCAAAAAAAGCGTGAAAAACCTTACTTGCCATAGAACTTGATACTTCTtctc                                                                            <  1:306387/130‑1 (MQ=255)                                                              aCCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAg                                                                         >  1:6034/1‑134 (MQ=255)                                                              aCCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAg                                                                         <  2:6034/134‑1 (MQ=255)                                                                          ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGTAACTTCTTCTCCGGTTCTaaaacgttccctt                                                        >  2:306387/1‑130 (MQ=255)                                                                          ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTg                                                        >  2:411849/1‑139 (MQ=255)                                                                          ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTg                                                        >  1:514293/1‑139 (MQ=255)                                                                                  aaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAAttttt                                                >  2:661011/1‑139 (MQ=255)                                                                                  aaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAAttttt                                                >  2:556297/1‑139 (MQ=255)                                                                                    cAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTaa                                              <  1:829020/139‑1 (MQ=255)                                                                                    cAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTaa                                              >  1:137652/1‑139 (MQ=255)                                                                                       cAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACtg                                           >  2:765187/1‑139 (MQ=255)                                                                                                          tAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTa                                   <  1:352243/128‑1 (MQ=255)                                                                                                          tAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTa                                   >  2:352243/1‑128 (MQ=255)                                                                                                           aaaGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGa                       >  1:824931/1‑139 (MQ=255)                                                                                                             aGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGaaa                     >  1:344915/1‑139 (MQ=255)                                                                                                              gTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAAt                    >  2:17886/1‑139 (MQ=255)                                                                                                               tctttcttACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATc                   <  1:359147/139‑1 (MQ=255)                                                                                                                              aTTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGat    <  2:905073/139‑1 (MQ=255)                                                                                                                              aTTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGat    <  2:514293/139‑1 (MQ=255)                                                                                                                                tttGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGatat  <  2:137652/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   TAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATAT  >  NC_000913/4268674‑4268939

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