Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,716,441	A→T	I261N (ATC→AAC)	tyrS ←	tyrosine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,716,441	0	A	T	96.1%	164.6 / ‑1.9	51	I261N (ATC→AAC)	tyrS	tyrosine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (25/24);  minor base C (1/1);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

TCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGAT  >  NC_000913/1716309‑1716576                                                                                                                                     |                                                                                                                                        tCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGaa                                                                                                                                  <  2:492285/140‑1 (MQ=255)   cGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACt                                                                                                                                >  2:440812/1‑140 (MQ=255)              aCTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAATGTTGTTCCAGAACTGGTAGAATTTg                                                                                                                     >  2:712308/1‑140 (MQ=255)                  ggCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTgg                                                                                                                              <  1:117734/127‑1 (MQ=255)                  ggCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTgg                                                                                                                              >  2:117734/1‑127 (MQ=255)                  ggCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACg                                                                                                                 >  2:393122/1‑140 (MQ=255)                  ggCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACg                                                                                                                 >  2:275119/1‑140 (MQ=255)                   gcgcgcGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTgg                                                                                                                              >  2:780014/1‑126 (MQ=255)                   gcgcgcGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTgg                                                                                                                              <  1:780014/126‑1 (MQ=255)                   gcgcgcGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACgg                                                                                                                >  2:71444/1‑140 (MQ=255)                   gcgcgcGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACg                                                                                                                 <  1:172259/139‑1 (MQ=255)                   gcgcgcGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACg                                                                                                                 >  2:172259/1‑139 (MQ=255)                          gTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGtt                                                                                                         >  2:155801/1‑140 (MQ=255)                          gTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGtt                                                                                                         <  1:712308/140‑1 (MQ=255)                                   cGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGa                                                                                                <  1:275119/140‑1 (MQ=255)                                     cTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATc                                                                                              >  1:113283/1‑140 (MQ=255)                                        tttttATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCaa                                                                                           >  2:801240/1‑140 (MQ=255)                                         ttttATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAAc                                                                                          <  2:632912/140‑1 (MQ=255)                                             aTCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGa                                                                                                <  2:586572/130‑1 (MQ=255)                                             aTCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGa                                                                                                >  1:586572/1‑130 (MQ=255)                                                   ttcttcCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTgcgc                                                                                <  1:393122/140‑1 (MQ=255)                                                         cAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCtt                                                                          <  2:774973/140‑1 (MQ=255)                                                           gggCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCa                                                                        >  1:588467/1‑140 (MQ=255)                                                               gTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGttt                                                                    >  2:783786/1‑140 (MQ=255)                                                                  gATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTAc                                                                 <  1:155801/140‑1 (MQ=255)                                                                             tGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTg                                                      >  1:159457/1‑140 (MQ=255)                                                                                tCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCa                                                   <  1:440812/140‑1 (MQ=255)                                                                                 cATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCAt                                                  <  1:125099/140‑1 (MQ=255)                                                                                 cATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCAt                                                  <  1:276308/140‑1 (MQ=255)                                                                                  aTAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATc                                                 >  2:578674/1‑140 (MQ=255)                                                                                  aTAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATc                                                 >  2:784391/1‑140 (MQ=255)                                                                                   tAAAGGTGAAGAACTCCAGG‑AGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTg                                               <  1:754672/140‑1 (MQ=255)                                                                                    aaaGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTg                                               <  1:784391/140‑1 (MQ=255)                                                                                    aaaGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTg                                               <  1:783786/140‑1 (MQ=255)                                                                                    aaaGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTg                                               <  2:588467/140‑1 (MQ=255)                                                                                       ggTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCtt                                            <  1:578674/140‑1 (MQ=255)                                                                                            agaaCTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCTGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTg                                       >  1:54788/1‑140 (MQ=255)                                                                                            agaaCTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTg                                       >  1:596094/1‑140 (MQ=255)                                                                                            agaaCTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTg                                       >  2:148471/1‑140 (MQ=255)                                                                                              aaCTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGAt                                     <  1:71444/140‑1 (MQ=255)                                                                                                   cAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCg                                >  1:725269/1‑140 (MQ=255)                                                                                                       aaGCGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAAc                            <  2:553189/140‑1 (MQ=255)                                                                                                        aGCGGTAAACGTCGACATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACg                           >  2:586478/1‑140 (MQ=255)                                                                                                          cGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGt                         <  2:159457/140‑1 (MQ=255)                                                                                                          cGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGt                         >  2:151602/1‑140 (MQ=255)                                                                                                          cGGTAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGt                         >  1:745845/1‑140 (MQ=255)                                                                                                          cGGTAAACGTCGGCATCCGCAGTGTTGCTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGaa                             >  2:211572/1‑136 (MQ=255)                                                                                                          cGGTAAACGTCGGCATCCGCAGTGTTGCTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGaa                             <  1:211572/136‑1 (MQ=255)                                                                                                             tAAACGTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAg                      >  2:135626/1‑140 (MQ=255)                                                                                                                  gTCGGCATCCGCAGTGTTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCaa                 <  1:135626/140‑1 (MQ=255)                                                                                                                               gtgtTGTTCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTg    <  2:113283/140‑1 (MQ=255)                                                                                                                                 gttgttCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCTAACACCTGATTCTGAt  >  2:577520/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                        TCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGAT  >  NC_000913/1716309‑1716576

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