Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,820,792	T→C	D407G (GAC→GGC)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,820,792	0	T	C	86.0%	127.9 / 10.2	50	D407G (GAC→GGC)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (21/22);  total (25/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 = 1.00e+00

ATGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTT  >  NC_000913/2820646‑2820924                                                                                                                                                   |                                                                                                                                     aTGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCa                                                                                                                                    >  4:184663/1‑149 (MQ=255)  tGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCaa                                                                                                                                   >  5:66303/1‑149 (MQ=255)            acacGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGcc                                                                                                                         <  5:106232/149‑1 (MQ=255)              acGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCa                                                                                                                        >  3:43972/1‑148 (MQ=255)                ggATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     <  6:129760/149‑1 (MQ=255)                     aTTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTc                                                                                                                >  5:51285/1‑149 (MQ=255)                      ttGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCg                                                                                                               <  6:66303/149‑1 (MQ=255)                         cGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTa                                                                                                             >  7:82792/1‑148 (MQ=255)                          gTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAAc                                                                                                                                  <  6:36717/126‑1 (MQ=255)                          gTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAAc                                                                                                                                  >  5:36717/1‑126 (MQ=255)                                   ggCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGaaa                                                                                                   <  7:86492/148‑1 (MQ=255)                                     cACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCa                                                                                                <  4:43972/149‑1 (MQ=255)                                         aaaGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGtt                                                                                             >  6:24803/1‑148 (MQ=255)                                         aaaGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGtt                                                                                             <  5:24803/148‑1 (MQ=255)                                           aGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCa                                                                                          <  3:30917/149‑1 (MQ=255)                                             ccGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     <  2:210693/120‑1 (MQ=255)                                             ccGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     <  2:207491/120‑1 (MQ=255)                                             ccGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     >  1:210693/1‑120 (MQ=255)                                             ccGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     >  1:207491/1‑120 (MQ=255)                                                     ttCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGc                                                                                 >  2:113977/1‑148 (MQ=255)                                                       cgcgcgcgcgACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGc                                                                                 <  1:113977/146‑1 (MQ=255)                                                        gcgcgcgcgACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTa                                                                                                             >  1:14036/1‑117 (MQ=255)                                                        gcgcgcgcgACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTa                                                                                                             <  2:14036/117‑1 (MQ=255)                                                           cgcgcgACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCAc                                                                          <  2:212572/149‑1 (MQ=255)                                                                      gctgctCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGCCGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGtttt                                                               >  3:138473/1‑149 (MQ=255)                                                                      gctgctCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGtttt                                                               <  8:82792/149‑1 (MQ=255)                                                                         gctCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGc                                                             >  2:62344/1‑148 (MQ=255)                                                                             ttCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAg                                                          <  1:62344/147‑1 (MQ=255)                                                                              tCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGctc                                                       >  2:163400/1‑149 (MQ=255)                                                                                 aTTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCtt                                                     >  8:42227/1‑148 (MQ=255)                                                                                        cTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCATATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGc                                             <  3:184663/149‑1 (MQ=255)                                                                                         ttCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGtttt                                                               >  6:76933/1‑130 (MQ=255)                                                                                          tCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCaa                                           >  5:146191/1‑149 (MQ=255)                                                                                             aaaCCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGc                                        <  1:163400/149‑1 (MQ=255)                                                                                                    cTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCa                                                                                          >  6:35872/1‑92 (MQ=255)                                                                                                    cTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCa                                                                                          <  5:35872/92‑1 (MQ=255)                                                                                                         tCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTc                            <  6:146191/149‑1 (MQ=255)                                                                                                         tCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTc                            >  4:6739/1‑149 (MQ=255)                                                                                                             cTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCaga                        <  3:6739/149‑1 (MQ=255)                                                                                                                        cgcTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCaaa             <  7:42227/149‑1 (MQ=255)                                                                                                                         gcTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCaaa             >  4:120374/1‑148 (MQ=255)                                                                                                                         gcTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCaaa             <  3:120374/148‑1 (MQ=255)                                                                                                                         gcTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAAc            <  6:51285/149‑1 (MQ=255)                                                                                                                           tCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGAACAGACGCAAAGAAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTg          >  7:54281/1‑149 (MQ=255)                                                                                                                                    aCATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTcttctt  >  5:42912/1‑148 (MQ=255)                                                                                                                                     cATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTcttctt  >  6:93927/1‑147 (MQ=255)                                                                                                                                      aTCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGa                                    <  7:22192/112‑1 (MQ=255)                                                                                                                                      aTCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGa                                    >  8:22192/1‑112 (MQ=255)                                                                                                                                         aGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCa                          <  2:64807/119‑1 (MQ=255)                                                                                                                                         aGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCa                          >  1:64807/1‑119 (MQ=255)                                                                                                                                                   |                                                                                                                                     ATGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTT  >  NC_000913/2820646‑2820924

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