Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	24,992	G→A	6.8%	A868T (GCT→ACT)	ileS →	isoleucyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	24,992	0	G	A	6.8%	179.8 / 2.7	59	A868T (GCT→ACT)	ileS	isoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base A (2/2);  total (28/31)
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.60e-01

AACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGC  >  NC_000913/24848‑25134                                                                                                                                                 |                                                                                                                                               aaCAAAGTCATTGAGCAAGTGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTAc                                                                                                                                            >  1:40879/1‑149 (MQ=255) aaCAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAATCTGACCTCGCGGGGCGATGAATTACGATTTGTCCAGTTGACCTCCGGCGCTAc                                                                                                                                            >  3:14422/1‑149 (MQ=255) aaCAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTAc                                                                                                                                            >  1:245604/1‑149 (MQ=255)    aaaGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGt                                                                                                                                         >  2:12281/1‑149 (MQ=255)          aTTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGa                                                                                                                                   >  1:269828/1‑149 (MQ=255)                 aaGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAac                                                                                                                            <  2:245604/149‑1 (MQ=255)                       gTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATaa                                                                                                                             >  1:230411/1‑142 (MQ=255)                       gTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATaa                                                                                                                             <  2:230411/142‑1 (MQ=255)                           agACATGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTAGGCAGAACCGGAAATGTCGGCGAACCTGCCCGCGCTGGGCGATGAATTACTCTTTGTCCTTTGGACCTTCTGCGCTACCGTTTCAGACTATAACTCCGCACatg                                                                                                                  >  4:273569/2‑149 (MQ=255)                               aaGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGa                                                                                                              >  5:165556/1‑149 (MQ=255)                                aGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAt                                                                                                             <  5:85255/149‑1 (MQ=255)                                 gAAAGTGGGTGGGTCGCTGGACTCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATg                                                                                                            <  4:14422/149‑1 (MQ=255)                                 gAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATg                                                                                                            >  7:259274/1‑149 (MQ=255)                                        ggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagca                                                                                                     <  2:40879/149‑1 (MQ=255)                                                gCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGaa                                                                                              <  3:59433/148‑1 (MQ=255)                                                    gAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACt                                                                                         >  7:139631/1‑149 (MQ=255)                                                     aaGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTc                                                                                        <  8:139631/149‑1 (MQ=255)                                                       gCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCa                                                                                       >  2:255602/1‑148 (MQ=255)                                                           cAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAggg                                                                                  >  8:78148/1‑149 (MQ=255)                                                                   ttGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGc                                                                                                 <  4:82689/126‑1 (MQ=255)                                                                   ttGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGc                                                                                                 >  3:82689/1‑126 (MQ=255)                                                                        tGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTAc                                                                                          >  6:43630/1‑128 (MQ=255)                                                                        tGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTAc                                                                                          <  5:43630/128‑1 (MQ=255)                                                                          cAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGtt                                                                    <  3:64332/148‑1 (MQ=255)                                                                           aGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTTTAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGa                                                                  <  8:258673/149‑1 (MQ=255)                                                                           aGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGa                                                                  <  1:12281/149‑1 (MQ=255)                                                                           aGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGa                                                                  <  8:142968/149‑1 (MQ=255)                                                                            gAACCTGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAg                                                                 <  1:255602/149‑1 (MQ=255)                                                                                cGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTaaa                                                             <  8:259274/149‑1 (MQ=255)                                                                                cGGAACTGTCGGCGAAACTGACCGCGCTGGGCGAGGAATTCCGATTTGGCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTaaa                                                             <  6:51780/149‑1 (MQ=255)                                                                                  gAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCt                                                                                                          >  3:1208/1‑102 (MQ=255)                                                                                  gAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCt                                                                                                          <  4:1208/102‑1 (MQ=255)                                                                                        tCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAgg                                                     <  7:78148/149‑1 (MQ=255)                                                                                        tCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAg                                                      <  6:165556/148‑1 (MQ=255)                                                                                          ggCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTc                                                                         <  2:382628/127‑1 (MQ=255)                                                                                          ggCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTc                                                                         >  1:382628/1‑127 (MQ=255)                                                                                           gCGAAACTTACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTACGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCATAGCGAAGTACTCAAAGGGCTGAAATTCGCGttta                                                                  >  5:5428/1‑131 (MQ=255)                                                                                           gCGAAACTGACCGCGCTTGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGATCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGAGAAGTACTCAAAGGGCTGAAAGTCGCGTTGa                                                                  <  6:5428/133‑1 (MQ=255)                                                                                           gCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTga                                                  <  2:269828/149‑1 (MQ=255)                                                                                             gAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCACTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTgag                                                 >  5:202324/1‑148 (MQ=255)                                                                                             gAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCACTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTgag                                                 >  5:197151/1‑148 (MQ=255)                                                                                              aaaCTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCACTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGaa                                               <  6:202324/149‑1 (MQ=255)                                                                                              aaaCTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCACTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGaa                                               <  6:197151/149‑1 (MQ=255)                                                                                                           tGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGtt                                                                    <  2:3790/115‑1 (MQ=255)                                                                                                               cGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTgg                               >  7:363386/1‑148 (MQ=255)                                                                                                                gATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCa                             >  5:135609/1‑149 (MQ=255)                                                                                                                  tGAATTACGATTTGTCCTGGTGCCCTCCGGCGCTACCGGTGCAGACTTTCCCGACGCACCGGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGa                                                                  <  7:267693/110‑1 (MQ=255)                                                                                                                    aaTTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAg                                                      <  4:186249/120‑1 (MQ=255)                                                                                                                    aaTTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAg                                                      >  3:186249/1‑120 (MQ=255)                                                                                                                     aTTACGATTTGTCCTGGTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTAATCAAAGGCCTTAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTTCCCACGCTGCTGGCACTaca                        <  3:273569/149‑1 (MQ=255)                                                                                                                      ttACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAgg                                                                                   >  7:40005/1‑89 (MQ=255)                                                                                                                      ttACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAgg                                                                                   <  8:40005/89‑1 (MQ=255)                                                                                                                      ttACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTc                                                                         >  7:133653/1‑99 (MQ=255)                                                                                                                            tttGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGa                 >  7:382944/1‑149 (MQ=255)                                                                                                                                    gTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGc          >  2:55486/1‑148 (MQ=255)                                                                                                                                        aCCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTaaa                                                             >  7:265551/1‑93 (MQ=255)                                                                                                                                        aCCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTaaa                                                             <  8:265551/93‑1 (MQ=255)                                                                                                                                         ccTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTg    <  8:382944/149‑1 (MQ=255)                                                                                                                                           tCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGc  >  1:149502/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                               AACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGC  >  NC_000913/24848‑25134

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