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.9%	98.6 / NA	32	P3S (CCT→TCT)	prfA	peptide chain release factor RF1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/13);  minor base A (1/0);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTGAGATGGCGCAGGATGAACTG  >  NC_000913/1264781‑1265257                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                cgATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGc                                                                                                                                                                                                                                    >  1:280449/1‑251 (MQ=255)           aCCGCCAAATCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACGTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGaa                                                                                                                                                                                                                          <  1:368317/251‑1 (MQ=255)                    gCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATCATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCAt                                                                                                                                                                                                                 <  1:141406/251‑1 (MQ=255)                                                           aTTATGAAGGAACTGGCAAGGAAACTGTCTAACCGCTTGATCAATGCGCCAACGAAATCACTTAAACAGGCCGCACGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGATGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGACCTACATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGGGAcgcg                                                                                                                                                                          >  2:327240/1‑251 (MQ=255)                                                                            aTGGAAACTGACTAAACGCTTGATCCACGCGCCAACGAAATCCCTTAAACAGGCCGCCCGTGACGGGGATAACGAACACCTAAATCTTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGCACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTAACGCCGAcc                                                                                                                                                          >  2:229728/1‑250 (MQ=255)                                                                                     gactatacGCTTGATCCATGCGCCAACGAAATCACGTCAACAGGCCGCCCGTTTCGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCGGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTGGCCAAACTGGACGCCCTTCAGGAACGCCATGAAGAAGTTCAGGCGTTGCTGGTTGACGCGCAAACTATCACTGCCCAGGAACGt                                                                                                                                                  <  2:138791/242‑1 (MQ=255)                                                                                         agccGCGTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCGTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTcgc                                                                                                                                            <  1:191680/249‑1 (MQ=255)                                                                                            cGCTTGATCCATGCGCCAACGAAATCACTTCAACCGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCgcgc                                                                                                                                          >  1:63537/1‑250 (MQ=255)                                                                                              cTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCAtt                                                                                                                                       >  1:524846/1‑251 (MQ=255)                                                                                                 gATCCAGGCGCCAACGAAATCACTGCACCAGGCCGCACGCGACGGGGATAACGAACGCCCGAATATTCAGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTGCTTTTTTAACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTTGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATc                                                                                                                                    <  2:479642/251‑1 (MQ=255)                                                                                                        gcgcCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAAAGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGa                                                                                                                              >  1:331934/1‑250 (MQ=255)                                                                                                              aCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTGTTACAGGGTGCATTTACGCCTATGAAGTCTGCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGGTGTGTGACGCGCAAACTATCGCCGACCAGGAACGGTTTCGCGCATTAGCACGCGAATATgcg                                                                                                                       <  2:63537/251‑1 (MQ=255)                                                                                                                  aaTCAATTAAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGCGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTAAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGCTTTCGCGCATGATCACGCGAATATGCGCAGt                                                                                                                   >  2:440263/1‑251 (MQ=255)                                                                                                                      aCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCTGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAg                                                                                                               <  2:378464/251‑1 (MQ=255)                                                                                                                             cAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGtt                                                                                                        >  2:422950/1‑251 (MQ=255)                                                                                                                              aGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGttt                                                                                                       <  1:229631/251‑1 (MQ=255)                                                                                                                                  cgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAATGATGTTTcgcg                                                                                                   >  1:431548/1‑251 (MQ=255)                                                                                                                                      cGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGt                                                                                               >  2:192556/1‑251 (MQ=255)                                                                                                                                                      aacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCGATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAAc                                                                               <  2:280449/251‑1 (MQ=255)                                                                                                                                                        cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAATGATGTTTCGCGCTGTTTTACCGACTGGCAACAg                                                                             <  2:431548/251‑1 (MQ=255)                                                                                                                                                                            cAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAAc                                                         >  2:84486/1‑251 (MQ=255)                                                                                                                                                                                                 gTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCgatg                                     >  1:507184/1‑250 (MQ=255)                                                                                                                                                                                                       catTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACAGACTGGCAACAGGTTCAGTAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTg                               >  2:306062/1‑250 (MQ=255)                                                                                                                                                                                                            tCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCt                                                                                                 >  1:53307/1‑179 (MQ=255)                                                                                                                                                                                                            tCTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGGTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCt                                                                                                 <  2:53307/179‑1 (MQ=255)                                                                                                                                                                                                             cTTTTTTTACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCg                        >  1:396437/1‑251 (MQ=255)                                                                                                                                                                                                              tttttttACAGGGTGCATTTACGCCTATGAAGTCTTATATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGAATGGCAACAGGTTCAGGAAGATCTCGAACCCGCACAGAAGATGCTCGATGATCCTGAAATGCgg                       >  2:508005/1‑250 (MQ=255)                                                                                                                                                                                                                   ttACAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTgaga                   <  2:3409/250‑1 (MQ=255)                                                                                                                                                                                                                      cAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGACGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTGAGATGgc               >  2:28756/1‑251 (MQ=255)                                                                                                                                                                                                                      cAGGGTGCATTTACGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGCTGATCCTGAAATGCGTGAGATGgc               >  2:465011/1‑251 (MQ=255)                                                                                                                                                                                                                         ggTGCATTTACGCCTATGAAGACTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGACGTTCAGGCGTTGCTGGGTGACGCGCACACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGTGTGAGATGgcgc             >  2:265699/1‑250 (MQ=255)                                                                                                                                                                                                                                   cGCCTATGAAGTCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTGAGATGGCGCAGGATGAACTg  <  2:524846/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                CGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGATGTTTCGCGCTGTTTTACCGACTGGCAACAGGTTCAGGAAGATATCGAAACCGCACAGATGATGCTCGATGATCCTGAAATGCGTGAGATGGCGCAGGATGAACTG  >  NC_000913/1264781‑1265257

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: ‑