Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,950,367	A→T	D52E (GAT→GAA)	aspS ←	aspartyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,950,367	0	A	T	100.0%	218.8 / NA	66	D52E (GAT→GAA)	aspS	aspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (21/45);  total (21/45)

AATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACC  >  NC_000913/1950274‑1950459                                                                                              |                                                                                             aaTATTTTTTTTGTCATGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGTCCAGCTTTAACGCGTCCGCACGTTCCGGa                                                                                        >  1:1127588/1‑100 (MQ=255) aaTATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGa                                                                                        <  1:1135273/100‑1 (MQ=255)   tatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGTCCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATc                                                                                      <  1:1526658/100‑1 (MQ=255)   tatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATc                                                                                      <  2:1136971/100‑1 (MQ=255)   tatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATc                                                                                      >  2:1504282/1‑100 (MQ=255)    atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCTGATCg                                                                                     >  2:885110/1‑100 (MQ=255)     tttttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGa                                                                                    <  1:768422/100‑1 (MQ=255)       tttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGaaa                                                                                  <  1:157283/100‑1 (MQ=255)       tttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGaaa                                                                                  <  2:1593191/100‑1 (MQ=255)       tttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGaaa                                                                                  <  2:1099319/100‑1 (MQ=255)         tttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGaaaaa                                                                                <  2:1796798/100‑1 (MQ=255)          ttCGTCATGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAAt                                                                               <  2:1139793/100‑1 (MQ=255)          ttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAAt                                                                               <  1:62401/100‑1 (MQ=255)          ttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAAt                                                                               <  1:1351738/100‑1 (MQ=255)           tCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGATCTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATa                                                                              <  1:497999/100‑1 (MQ=255)             gTCACGCGCACGTACGTTGCCCGTTACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATAcc                                                                            <  2:824945/100‑1 (MQ=255)             gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCCGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATAc                                                                             >  1:1140457/1‑99 (MQ=255)             gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATAcc                                                                            >  2:742820/1‑100 (MQ=255)             gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATAcc                                                                            <  2:1208435/100‑1 (MQ=255)              tCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCt                                                                           >  2:257950/1‑100 (MQ=255)               cACGCGCACGTACGGTGCCCGTGATCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTg                                                                          <  2:1528786/100‑1 (MQ=255)                 cgcgCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCa                                                                        <  2:1505302/100‑1 (MQ=255)                 cgcgCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCa                                                                        >  2:378597/1‑100 (MQ=255)                  gcgcACGTACGGTGCCCGTGACCTGAATGAAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCAc                                                                       <  1:362527/100‑1 (MQ=255)                      aCGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATa                                                                   <  2:1007409/100‑1 (MQ=255)                      aCGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATa                                                                   <  2:1425864/100‑1 (MQ=255)                        gTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATAcc                                                                 <  2:169269/100‑1 (MQ=255)                         tACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCt                                                                >  1:729961/1‑100 (MQ=255)                          aCGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCtt                                                               <  1:175867/100‑1 (MQ=255)                            ggTGCCCGGGACCTGAAGGCAGAAGGTAGTACGCAGTTCAGAGGGCAGCTTTAGCGCGTCCGAACGTTCCGAATCGAAAAATACCTGCACGATACCTTcg                                                             <  1:977976/100‑1 (MQ=255)                             gTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTcgc                                                            <  1:837111/100‑1 (MQ=255)                               gCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTcgcg                                                           <  1:515519/99‑1 (MQ=255)                                 ccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTc                                                        >  2:314240/1‑100 (MQ=255)                                   gTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTAAAAAAACAAAAAATACCTGCACGATACCTTCGCGGTcgc                                                      <  1:257949/100‑1 (MQ=255)                                    ctaCCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTcgcg                                                     <  2:615176/98‑1 (MQ=255)                                      aCCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTcgcgag                                                   >  1:615172/1‑98 (MQ=255)                                        cTGAATGCAGAACGCATTACGCAGATCAGAAGCCAGATATAACGCTTCAGCACGTTCCGGATCGAAAAATACCTGCACGATAGCTTTGCGGTCGGGGata                                                 >  1:1119845/1‑100 (MQ=255)                                         tGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCatat                                                >  1:175405/1‑100 (MQ=255)                                            aTGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCgn                                             <  2:1467772/100‑2 (MQ=255)                                              gCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATg                                           >  1:895451/1‑100 (MQ=255)                                                  aaCTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGa                                       <  2:1507199/100‑1 (MQ=255)                                                   aCTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGAt                                      >  2:1629968/1‑100 (MQ=255)                                                     tCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCa                                    <  2:605168/100‑1 (MQ=255)                                                      cATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAg                                   >  2:1062505/1‑100 (MQ=255)                                                      cATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAg                                   <  1:475331/100‑1 (MQ=255)                                                        ttACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc                                 <  2:313805/100‑1 (MQ=255)                                                        ttACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc                                 <  2:1626003/100‑1 (MQ=255)                                                           cGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTAc                              >  2:1220693/1‑100 (MQ=255)                                                           cGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTAc                              >  2:1211470/1‑100 (MQ=255)                                                            gCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTAcc                             >  2:1526542/1‑100 (MQ=255)                                                               gTTCAGAGGACAGCTTTAACGCGTCCGCACGTTCTGGATCGAAAAATACCTGCACGATACCTTCGCTGTCGCGCATATCGATGAAGATCAGGCTACCAAg                          <  1:175056/100‑1 (MQ=255)                                                                ttCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGa                         <  1:1050535/100‑1 (MQ=255)                                                                ttCAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGa                         <  1:208801/100‑1 (MQ=255)                                                                  cAGAGGCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATc                       <  2:1127599/100‑1 (MQ=255)                                                                      ggCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCacga                   >  1:582645/1‑100 (MQ=255)                                                                       gCCAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCacgac                  <  1:73611/100‑1 (MQ=255)                                                                         cAGCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCacgacga                >  1:521041/1‑100 (MQ=255)                                                                           gCTTTAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCacgacgacg              <  2:359/100‑1 (MQ=255)                                                                              ttAACGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGtt           <  1:1812045/100‑1 (MQ=255)                                                                                 aCGCGTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCTATGAAGATCAGGCTACCAAGATCACGACGACGGTTGAc        <  1:757223/100‑1 (MQ=255)                                                                                  cgcgTCTGTACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGAcc       <  2:765808/100‑1 (MQ=255)                                                                                     gTCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCaa    <  1:976435/100‑1 (MQ=255)                                                                                      tCCGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAAc   >  2:923018/1‑100 (MQ=255)                                                                                       ccGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAAcc  <  1:611910/100‑1 (MQ=255)                                                                                       ccGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAAcc  <  2:1254354/100‑1 (MQ=255)                                                                                       ccGCACGTTCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAAcc  <  1:885104/100‑1 (MQ=255)                                                                                              |                                                                                             AATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACC  >  NC_000913/1950274‑1950459

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