Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,265,018	C→T	P3S (CCT→TCT)	prfA →	peptide chain release factor RF1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,265,018	0	C	T	96.3%	76.6 / ‑5.4	27	P3S (CCT→TCT)	prfA	peptide chain release factor RF1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/14);  minor base G (0/1);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

AGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTGAGATGGCGCAGG  >  NC_000913/1264771‑1265249                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        agCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCtt                                                                                                                                                                                                                                      >  1:328164/1‑251 (MQ=255)    aaGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGtcttct                                                                                                                                                                                                                                    >  2:94203/1‑250 (MQ=255)     aGTTCGCAATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTAt                                                                                                                                                                                                                                  >  1:167137/1‑251 (MQ=255)               gAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTTGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCaa                                                                                                                                                                                                                         <  1:433882/250‑1 (MQ=255)                            aaGCGTTAGCGGCCCTTGAGCAGGGCGGGGACGGGCAAGCGATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAGCTAAATCGCTTCAACAGGCCTCCCGTGACGGGGATAACGAACGCCTGAATATTCTTCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTGTTTTACAGGGTGCATTTACGCCTATGAAGTCTTTTATCGTTGCCAAACGGGAAGCCCTGc                                                                                                                                                                                                           <  2:257226/251‑1 (MQ=255)                                            tGAGCATGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATAGTCTGCGCGACAGCCTCGGGCTGGAGTATCAGTACATCAGTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATgaag                                                                                                                                                                                            <  1:54464/250‑1 (MQ=255)                                                                                      aTGGAAACTGACTAACCGCGTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCGGAATATTCGGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGGCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTGCAGGCGTTGCTGGGTGACGCTCAAACTATCGTCGACCa                                                                                                                                                 <  2:350008/251‑1 (MQ=255)                                                                                               gACTAACCGCTTGATCCAGGCGCCAACGAAATCACTTCAACAGGCAGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAGCTATCGCCGACCAGGAACGttt                                                                                                                                        <  2:92558/251‑1 (MQ=255)                                                                                                            aTCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCa                                                                                                                           <  1:273642/251‑1 (MQ=255)                                                                                                                       aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAACCTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGCCCAGGAACGTTTTCGCGCATTATCACGCGAATATgc                                                                                                                >  2:35831/1‑251 (MQ=255)                                                                                                                       aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATgc                                                                                                                <  1:97204/251‑1 (MQ=255)                                                                                                                            aaTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGt                                                                                                           <  2:79988/251‑1 (MQ=255)                                                                                                                            aaTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGt                                                                                                           <  1:74197/251‑1 (MQ=255)                                                                                                                                 cTTCAACAGGCCGCCCGGGACGGGGATAACGACAGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGGTGCTGGGAGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTAGCACGCGAATATGCGCAGTTAAgg                                                                                                      >  2:226837/1‑250 (MQ=255)                                                                                                                                      gcAGGCCGAACGTGACGGGGATAACGACCGCCTGGATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCCGTACATCATTTACTTATTTGACAGGATGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAAATTGAAGCCCTGCATGAACGCCCTGAAGGAGTTCAGGCGTTGCTGGGTGACGCGCAAACTCTCGCCGACCAGGAACGTTTTCGCGCTTGATCACGAGAATATGCGACGTTAAGTGgtgt                                                                                                 >  2:306269/2‑251 (MQ=255)                                                                                                                                           ccgccCGAGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGAAGTACATCATTTTCTTTTTTTACCGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGACCGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAACCTATCGCCGACCAGGAACGTTTTCGCGCATTATCATGCGAATATGCGCAGTTATGTGATGtttagc                                                                                            >  2:224879/1‑248 (MQ=255)                                                                                                                                                     cGGGGATAACGAACGCCTGAATATTCTGAGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCCTGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATTCGCAGTTAAGTGATGTTTCGCGATGTTTTAc                                                                                  >  2:381953/1‑251 (MQ=255)                                                                                                                                                                           atTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGaa                                                            <  1:94203/251‑1 (MQ=255)                                                                                                                                                                           atTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGaa                                                            <  2:413311/251‑1 (MQ=255)                                                                                                                                                                                                  ggAGTAGCAGTACATCATTTTCTTTTTTTACAGGCTGCATTTACGCCTATGAAGTCTTCAATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGCAGACGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGGTTTCGCGCATTATCACGCGACTATGCACAGTTACGCGATGTTTCGCCCTGTTTTACCGACTGGCAACCGGTTCAGGAAGATATCGACACCGCCCAgatgac                                     >  2:238159/1‑250 (MQ=255)                                                                                                                                                                                                  ggAGTAGCAGTACATCATTTTCTTTTTTGACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGATACCGCACAgatgat                                     <  2:67058/251‑1 (MQ=255)                                                                                                                                                                                                    aGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGTATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGc                                   <  1:16086/251‑1 (MQ=255)                                                                                                                                                                                                                      tCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATg                  <  1:437687/250‑1 (MQ=255)                                                                                                                                                                                                                       cTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCg                >  1:340872/1‑251 (MQ=255)                                                                                                                                                                                                                           ttttACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTg              >  1:416358/1‑249 (MQ=255)                                                                                                                                                                                                                                  gggTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTGAGATGgcgc     >  2:228813/1‑251 (MQ=255)                                                                                                                                                                                                                                      gCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTGAGATGGCGCAgg  >  1:278325/1‑250 (MQ=255)                                                                                                                                                                                                                                      gCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTGAGATGGCGCAgg  <  1:55670/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        AGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTGAGATGGCGCAGG  >  NC_000913/1264771‑1265249

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