Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	940,380	T→C	8.3%	L318S (TTG→TCG)	serS →	seryl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	940,380	0	T	C	8.3%	140.4 / 4.1	48	L318S (TTG→TCG)	serS	seryl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base C (2/2);  total (23/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 = 6.85e-01

CTGAAGCCGGTTCATATGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCA  >  NC_000913/940233‑940524                                                                                                                                                    |                                                                                                                                                 cTGAAGCCGGTTCATATGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTg                                                                                                                                                 >  1:115910/1‑149 (MQ=255)  tGAAGCCGGTTCATATGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGc                                                                                                                                                >  1:205001/1‑149 (MQ=255)      gCCGGTTCATATGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                                                                                                                            >  2:44588/1‑149 (MQ=255)      gCCGGTTCATATGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                                                                                                                            >  1:6388/1‑149 (MQ=255)      gCCGGTTCATATGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTgg                                                                                                                                             >  3:357600/1‑148 (MQ=255)         ggTTCATATGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGcc                                                                                                                                          <  4:357600/148‑1 (MQ=255)                atGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgta                                                                                                                                   <  4:314935/148‑1 (MQ=255)                 tGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtacc                                                                                                                                 <  2:380539/149‑1 (MQ=255)                   gTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCACTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtaccgt                                                                                                                               >  1:291120/1‑149 (MQ=255)                      gTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTg                                                                                                                                              >  1:194264/1‑131 (MQ=255)                      gTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTg                                                                                                                                              <  2:194264/131‑1 (MQ=255)                           acCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                                                                                                                        <  3:211192/148‑1 (MQ=255)                             ccGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatcatc                                                                                                                      >  4:40217/1‑148 (MQ=255)                             ccGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATcc                                                                                                                     >  2:235990/1‑149 (MQ=255)                                tGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCtt                                                                                                                   <  3:228168/148‑1 (MQ=255)                                    cTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCa                                                                                                               >  3:69412/1‑148 (MQ=255)                                             aTGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACAt                                                                                                     <  1:156782/149‑1 (MQ=255)                                             aTGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACAt                                                                                                     <  1:186240/149‑1 (MQ=255)                                                  ccAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt                                                                                                >  2:114937/1‑149 (MQ=255)                                                    aGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCttt                                                                                              <  1:235990/149‑1 (MQ=255)                                                             aaGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTg                                                                                      >  3:279672/1‑148 (MQ=255)                                                                ttGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCaa                                                                                   >  3:265326/1‑148 (MQ=255)                                                                 tGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCaaaa                                                                                 <  2:389936/149‑1 (MQ=255)                                                                     aTGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTa                                                                             <  2:6388/149‑1 (MQ=255)                                                                      tGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTAc                                                                            <  1:114937/149‑1 (MQ=255)                                                                            aGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAATAAGACCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGACTTTGGCGCTTGCAAAACTTACGACCTg                                                                      <  2:291120/149‑1 (MQ=255)                                                                                 gTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGt                                                                 >  2:187410/1‑149 (MQ=255)                                                                                   gcgcCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTa                                                                <  4:69412/148‑1 (MQ=255)                                                                                       ccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGAt                                                           >  1:289634/1‑149 (MQ=255)                                                                                          gaagaCTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCAGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATccc                                                        <  1:297564/149‑1 (MQ=255)                                                                                                     ggcggcGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAac                                              <  3:284767/148‑1 (MQ=255)                                                                                                                      tGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATctct                            >  2:189400/1‑149 (MQ=255)                                                                                                                      tGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATctct                            <  1:44588/149‑1 (MQ=255)                                                                                                                        aCTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTc                          <  2:205001/149‑1 (MQ=255)                                                                                                                        aCTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTc                          <  1:187410/149‑1 (MQ=255)                                                                                                                                   aGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTCCTGCTCCAACg               >  2:371056/1‑149 (MQ=255)                                                                                                                                    gAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCgg                                                      <  3:285718/109‑1 (MQ=255)                                                                                                                                    gAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCgg                                                      >  4:285718/1‑109 (MQ=255)                                                                                                                                     aaaaaGTCCTGCAGTCGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTcc                         <  3:253780/137‑1 (MQ=255)                                                                                                                                     aaaaaGTCCTGCAGTCGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTcc                         <  3:148627/137‑1 (MQ=255)                                                                                                                                     aaaaaGTCCTGCAGTCGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTcc                         >  4:148627/1‑137 (MQ=255)                                                                                                                                     aaaaaGTCCTGCAGTCGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTcc                         >  4:253780/1‑137 (MQ=255)                                                                                                                                             cTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAacac                                            <  1:317637/110‑1 (MQ=255)                                                                                                                                             cTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAacac                                            >  2:317637/1‑110 (MQ=255)                                                                                                                                               gCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTcc   >  1:99738/1‑149 (MQ=255)                                                                                                                                                cAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTCCTGCTCCAACGTTTg           >  2:212892/1‑140 (MQ=255)                                                                                                                                                cAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTCCTGCTCCAACGTTTg           <  1:212892/140‑1 (MQ=255)                                                                                                                                                cAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCa  <  1:189400/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CTGAAGCCGGTTCATATGGTCGTGACACCCGTGGTCTGATCCGTATGCACCAGTTCGACAAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAAACTTACGACCTGGAAGTATGGATCCCGGCACAGAACACCTACCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCA  >  NC_000913/940233‑940524

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