Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,801,156	T→C	9.3%	D472G (GAC→GGC)	thrS ←	threonyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,801,156	0	T	C	9.3%	123.5 / 3.8	43	D472G (GAC→GGC)	thrS	threonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (18/21);  new base C (2/2);  total (20/23)
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.58e-01

TTCCATCGACCCCAGAATTGCGCGGTGAATCATTACCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGA  >  NC_000913/1801014‑1801302                                                                                                                                               |                                                                                                                                                   ttCCATCGACCCCAGAATTGCGCGGTGAATCATTACCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATAc                                                                                                                                               <  6:304161/148‑1 (MQ=255) ttCCATCGACCCCAGAATTGCGCGGTGAATCATTACCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATAc                                                                                                                                               <  5:305608/148‑1 (MQ=255)   ccATCGACCCCAGAATTGCGCGGTGAATCATTAGCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACa                                                                                                                                              <  1:14216/147‑1 (MQ=255)    cATCGACCCCAGAATTGCGCGGTGAATCATTACCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAgg                                                                                                                                            <  4:247898/148‑1 (MQ=255)                    gcgcGGTGAATCATTACCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTc                                                                                                                             <  4:26512/147‑1 (MQ=255)                                 ttACCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAgc                                                                                                               >  5:216192/1‑148 (MQ=255)                                  tACCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAgcg                                                                                                              >  7:139368/1‑148 (MQ=255)                                   aCCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTTCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAgcg                                                                                                              <  8:313882/147‑1 (MQ=255)                                            ttACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAAtttt                                                                                                                              >  3:229755/1‑122 (MQ=255)                                            ttACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAAtttt                                                                                                                              <  4:229755/122‑1 (MQ=255)                                            ttACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGa                                                                                                                          >  1:254340/1‑126 (MQ=255)                                            ttACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGa                                                                                                                          <  2:254340/126‑1 (MQ=255)                                                 ttCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGa                                                                                               <  6:216192/148‑1 (MQ=255)                                                 ttCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGa                                                                                               >  5:285991/1‑148 (MQ=255)                                                 ttCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGCCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGa                                                                                               >  2:160724/1‑148 (MQ=255)                                                    aTTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGCCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGatat                                                                                            >  8:243879/1‑148 (MQ=255)                                                             gCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGCCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACggg                                                                                   <  7:243879/148‑1 (MQ=255)                                                                 aCATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAAtttt                                                                                                                              >  3:273265/1‑101 (MQ=255)                                                                 aCATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAAtttt                                                                                                                              <  4:273265/101‑1 (MQ=255)                                                                   aTAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGCTGATATTCAAACGGGAtgttg                                                                             >  5:316429/1‑148 (MQ=255)                                                                     aagaagCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATgttgtt                                                                           >  1:211639/1‑148 (MQ=255)                                                                            gcTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTcc                                                                    <  8:63778/148‑1 (MQ=255)                                                                                         ggCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCa                                                       <  4:213876/148‑1 (MQ=255)                                                                                         ggCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCa                                                       <  2:211639/148‑1 (MQ=255)                                                                                             aaGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTc                                                   >  1:292185/1‑148 (MQ=255)                                                                                             aaGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTc                                                   >  7:100382/1‑148 (MQ=255)                                                                                              agagAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTcc                                                                    <  8:269340/130‑1 (MQ=255)                                                                                              agagAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTcc                                                                    >  7:269340/1‑130 (MQ=255)                                                                                              agagAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTcc                                                  >  8:283226/1‑148 (MQ=255)                                                                                               gagaAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCg                                                 >  8:210126/1‑148 (MQ=255)                                                                                               gagaAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCg                                                 >  4:78782/1‑148 (MQ=255)                                                                                                      ccAGCTGTACTGTACCGCACTGCCATGCACGAGCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGc                                          <  5:270852/148‑1 (MQ=255)                                                                                                        aGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGa                                                                                  >  1:285421/1‑106 (MQ=255)                                                                                                        aGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGa                                                                                  <  2:285421/106‑1 (MQ=255)                                                                                                             tactgtacCGCACTGCCATGCACGATCGAGGCAGCCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATc                                    <  1:160724/147‑1 (MQ=255)                                                                                                               ctgtacCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATccc                                  <  4:151460/147‑1 (MQ=255)                                                                                                                            ccATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTg                    <  8:139368/148‑1 (MQ=255)                                                                                                                                  aCGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATa              >  2:57124/1‑148 (MQ=255)                                                                                                                                     aTCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAg                                                                  >  4:214553/1‑93 (MQ=255)                                                                                                                                     aTCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAg                                                                  <  3:214553/93‑1 (MQ=255)                                                                                                                                      tCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGtt          <  6:102447/148‑1 (MQ=255)                                                                                                                                              gTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGa  <  4:6770/148‑1 (MQ=255)                                                                                                                                              gTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGa  >  3:6770/1‑148 (MQ=255)                                                                                                                                               tCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGttttt       <  7:267754/142‑1 (MQ=255)                                                                                                                                               tCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGttttt       >  8:267754/1‑142 (MQ=255)                                                                                                                                               |                                                                                                                                                   TTCCATCGACCCCAGAATTGCGCGGTGAATCATTACCGGTACTTTACGTTCATTGTCTTCGCCTACATAAGAAGCGCTCAGACGAGACGGCAAAGAGAAGTCCAGCTGTACTGTACCGCACTGCCATGCACGATCGAGGCAGTCATACAGGGTAAATTCAATTTTCGGACCGTAGAAAGCGCCTTCACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGA  >  NC_000913/1801014‑1801302

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