Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	312,189	C→T	G26G (GGC→GGT)	ykgL →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	312,189	0	C	T	97.5%	123.3 / NA	40	G26G (GGC→GGT)	ykgL	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/23);  minor base A (1/0);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACC  >  NC_000913/312053‑312314                                                                                                                                         |                                                                                                                              gTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGa                                                                                                                             <  2:239620/139‑1 (MQ=255)   ttttATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAg                                                                                                                           >  1:661895/1‑139 (MQ=255)           tAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaa                                                                                                                   >  1:668513/1‑139 (MQ=255)            aGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaac                                                                                                                  <  2:705676/139‑1 (MQ=255)                ccAACAGATTTTCCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                              <  2:653896/139‑1 (MQ=255)                ccAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                              <  2:747901/139‑1 (MQ=255)                 cAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                             <  1:647843/139‑1 (MQ=255)                      gATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGTTGAAGCTGGTCAACAAGTCAACAt                                                                                                        <  2:522511/139‑1 (MQ=255)                      gATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACAt                                                                                                        <  2:568019/139‑1 (MQ=255)                        tttACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                      <  2:668513/139‑1 (MQ=255)                        tttACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                      >  2:173965/1‑139 (MQ=255)                        tttACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                      >  2:160900/1‑139 (MQ=255)                                 tttATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATgg                                                                                             >  2:70814/1‑139 (MQ=255)                                          ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                             <  2:414118/114‑1 (MQ=255)                                          ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                             >  2:593750/1‑114 (MQ=255)                                          ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                             >  1:414118/1‑114 (MQ=255)                                          ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                             <  1:593750/114‑1 (MQ=255)                                          ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCgg                                                                                    >  2:12193/1‑139 (MQ=255)                                          ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCgg                                                                                    >  2:496641/1‑139 (MQ=255)                                            aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaa                                                                                                                   >  2:648558/1‑106 (MQ=255)                                            aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaa                                                                                                                   <  1:648558/106‑1 (MQ=255)                                            aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                  >  1:254136/1‑139 (MQ=255)                                                           aaTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAg                                                                   <  2:256012/139‑1 (MQ=255)                                                                     gcgcCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGa                                                                                                                             >  1:190177/1‑71 (MQ=255)                                                                      cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCaaa                                                        >  2:427146/1‑139 (MQ=255)                                                                           cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                   <  1:70814/139‑1 (MQ=255)                                                                           cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                   <  1:423771/139‑1 (MQ=255)                                                                           cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGAAGGCAGGATACAGGGGTGAAGCGGGTAAACAAGTCAACATTAAGATAATGGGATATACGGAGAGAAATGTCAGGCAAATGGCCTCGAATGAAc                                                   <  2:151285/139‑1 (MQ=255)                                                                               gggAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                               <  1:427146/139‑1 (MQ=255)                                                                                          gAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                   >  2:288262/1‑124 (MQ=255)                                                                                          gAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                   <  1:288262/124‑1 (MQ=255)                                                                                           aGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATccc                                   >  2:193606/1‑139 (MQ=255)                                                                                           aGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGAGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATccc                                   >  2:380069/1‑139 (MQ=255)                                                                                                            aaaaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                  <  2:47939/139‑1 (MQ=255)                                                                                                            aaaaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                  <  2:661895/139‑1 (MQ=255)                                                                                                            aaaaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                  <  1:12193/139‑1 (MQ=255)                                                                                                            aaaaGAGGGTAGATGCCAACATAAAGGGGTGAAGCTCGTCAACAAGTCAACATTACGAAAATGTAATATTCGGCGACAAACTTCAGGCAGTTGGCCTCAAATGAACAGGAAAAATATATCCCCAGGAAAATTAATGTTg                  <  1:8930/139‑1 (MQ=255)                                                                                                                        aTGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAaca      >  2:733303/1‑139 (MQ=255)                                                                                                                            cAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACAcc  <  1:741575/139‑1 (MQ=255)                                                                                                                            cAGCATACAGGGGTGAAGATGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACAcc  <  1:733303/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                              GTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACC  >  NC_000913/312053‑312314

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