Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,513,003	A→G	D212D (GAT→GAC)	trpS ←	tryptophan‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,513,003	0	A	G	100.0%	365.3 / NA	92	D212D (GAT→GAC)	trpS	tryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (53/39);  total (53/39)

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048                                                |                                              tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             >  2:1103845/1‑50 (MQ=255) tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             >  2:413806/1‑50 (MQ=255) tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             >  2:3916394/1‑50 (MQ=255)  cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  2:2684984/1‑50 (MQ=255)  cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  2:2964637/1‑50 (MQ=255)   ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           >  1:4272649/1‑50 (MQ=255)   ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           >  1:2197826/1‑50 (MQ=255)   ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            <  1:1770255/49‑1 (MQ=255)    gAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          <  1:3219310/50‑1 (MQ=255)    gAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  2:956380/1‑50 (MQ=255)     aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTcc                                         <  1:3628450/50‑1 (MQ=255)     aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:1929698/1‑49 (MQ=255)     aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  2:435610/1‑49 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  2:1456737/1‑50 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:552735/1‑48 (MQ=255)         aGTGACCGCACGTTTGATTTACTTCACTACCGATTTCGGGTCTTCcagc                                      <  2:4214831/49‑1 (MQ=39)            gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  2:1831168/1‑50 (MQ=255)            gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  2:775144/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:2113147/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:2229501/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  2:2423423/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:2622657/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  <  1:3866674/49‑1 (MQ=255)              ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCg                                <  1:1216218/50‑1 (MQ=255)               cGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  2:3488525/1‑50 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              >  2:4198648/1‑49 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  2:1267429/50‑1 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:97386/50‑1 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  2:2535819/50‑1 (MQ=255)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  2:138108/1‑50 (MQ=255)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  2:3561417/50‑1 (MQ=255)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  2:1896313/50‑1 (MQ=255)                   cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           >  2:631317/1‑50 (MQ=255)                   cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           <  1:4041928/50‑1 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:3501042/1‑50 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  2:31030/1‑50 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:503438/1‑50 (MQ=255)                     tttGATTTTCTTTACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:4282549/1‑50 (MQ=255)                     tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:4263493/1‑50 (MQ=255)                     tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  2:2144590/1‑50 (MQ=255)                     tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:2551706/1‑50 (MQ=255)                     tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:72639/1‑50 (MQ=255)                       tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGtt                        <  1:4049081/49‑1 (MQ=255)                        gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:243256/49‑1 (MQ=255)                        gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:3203457/49‑1 (MQ=255)                        gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      >  1:2950083/1‑50 (MQ=255)                        gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      >  2:4060438/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:1201245/49‑1 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:3253955/49‑1 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:2311882/49‑1 (MQ=255)                          ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  1:3318299/50‑1 (MQ=255)                            ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   >  1:3521349/1‑49 (MQ=255)                            ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   <  1:432639/49‑1 (MQ=255)                             tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGa                 >  2:3619767/1‑50 (MQ=255)                              cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATACCGTTATTGCGAt                >  2:2607510/1‑50 (MQ=39)                              cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  2:3686590/1‑50 (MQ=255)                              cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  2:3190395/1‑50 (MQ=255)                                tcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              >  1:3879446/1‑50 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  2:3615167/48‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:2541816/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             >  1:2718558/1‑50 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             <  2:3731435/49‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  2:667373/50‑1 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  2:1047976/1‑50 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  1:1370608/50‑1 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAGTGCGATTATCg           <  2:3374892/50‑1 (MQ=255)                                    tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:726115/1‑50 (MQ=255)                                    tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1693618/1‑50 (MQ=255)                                    tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  2:2568147/1‑50 (MQ=255)                                     aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  2:4098723/1‑49 (MQ=255)                                     aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  2:1839574/1‑49 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          <  1:2223713/48‑1 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCa        <  1:2588678/50‑1 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAg       <  2:4348350/50‑1 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAg       <  2:2726525/50‑1 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  2:3244510/1‑50 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     <  2:453505/50‑1 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:1501593/1‑50 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  2:1573860/1‑50 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:3583038/1‑50 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     <  1:2276360/50‑1 (MQ=255)                                          tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    >  1:1571943/1‑50 (MQ=255)                                          tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    >  2:4510039/1‑50 (MQ=255)                                           ttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   >  2:1270040/1‑50 (MQ=255)                                            tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  1:2461840/49‑1 (MQ=255)                                            tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  2:2173260/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  2:1022101/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:983418/1‑49 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  2:4508332/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  2:1474212/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:2637380/1‑49 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:2240699/1‑49 (MQ=255)                                                |                                              TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048

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