Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,802,309	G→C	P88A (CCG→GCG)	thrS ←	threonine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,802,309	0	G	C	100.0%	211.6 / NA	54	P88A (CCG→GCG)	thrS	threonine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (33/21);  total (33/21)

TCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTC  >  NC_000913/1802173‑1802444                                                                                                                                         |                                                                                                                                        tcgtAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCa                                                                                                                                      >  1:1033482/1‑140 (MQ=255)  cgtAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCaa                                                                                                                                     >  2:738609/1‑140 (MQ=255)               aGCAAGCTCATGTATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATcg                                                                                                                        >  2:499783/1‑140 (MQ=255)               aGCAAGCTCATGTATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATcg                                                                                                                        <  1:499783/140‑1 (MQ=255)               aGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATcg                                                                                                                        <  2:605517/140‑1 (MQ=255)               aGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATcg                                                                                                                        >  1:757221/1‑140 (MQ=255)                  aaGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGt                                                                                                                     <  1:1125552/140‑1 (MQ=255)                              ccGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGt                                                                                                         <  2:509557/140‑1 (MQ=255)                              ccGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGt                                                                                                         <  2:1067201/140‑1 (MQ=255)                              ccGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGt                                                                                                         <  2:714134/140‑1 (MQ=255)                               cGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTg                                                                                                        >  1:1013427/1‑140 (MQ=255)                                 cTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTgcg                                                                                                      <  2:757221/140‑1 (MQ=255)                                   tctcGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGca                                                                                                    <  1:738609/140‑1 (MQ=255)                                   tctcGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGca                                                                                                    <  1:684221/140‑1 (MQ=255)                                               cGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGGgcg                                                                                                      <  2:111398/126‑1 (MQ=255)                                               cGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGGgcg                                                                                                      >  1:111398/1‑126 (MQ=255)                                                 aCATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAAt                                                                                      <  2:1033482/140‑1 (MQ=255)                                                         cTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCa                                                                              >  2:1163847/1‑140 (MQ=255)                                                           ggggTAACGTGCGGTCAAGATCAGCGTAGTAAGAAAAACCGTTGTCGATAGCCGGTCCGATTGCCCTTTTGGTACGCGCCCAAAGTAGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTGCAGa                                                                            <  2:1337495/140‑1 (MQ=255)                                                           gggTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGa                                                                            <  1:597942/140‑1 (MQ=255)                                                             gTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGAcc                                                                          >  2:248968/1‑140 (MQ=255)                                                                 aCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACcttct                                                                      >  1:200616/1‑140 (MQ=255)                                                                 aCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACcttct                                                                      >  1:67492/1‑140 (MQ=255)                                                                 aCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACcttct                                                                      >  2:760369/1‑140 (MQ=255)                                                                      tggTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTc                                                                 >  2:812872/2‑140 (MQ=255)                                                                       ggTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCt                                                                >  1:223740/1‑140 (MQ=255)                                                                        gTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCtt                                                               >  1:1087370/1‑140 (MQ=255)                                                                           aaGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTgg                                                            <  1:1163847/140‑1 (MQ=255)                                                                                   cgtcgtAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATg                                                    >  1:492326/1‑140 (MQ=255)                                                                                   cgtcgtAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATg                                                    >  1:1111769/1‑140 (MQ=255)                                                                                   cgtcgtAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATg                                                    >  1:144822/1‑140 (MQ=255)                                                                                     tcgtAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGAt                                                  <  2:223740/140‑1 (MQ=255)                                                                                        taataaAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGa                                               >  1:461387/1‑140 (MQ=255)                                                                                        taataaAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGa                                               >  2:417746/1‑140 (MQ=255)                                                                                        taataaAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGa                                               >  2:519060/1‑140 (MQ=255)                                                                                                   gTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGcgt                                    >  2:728399/1‑140 (MQ=255)                                                                                                   gTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGcgt                                    >  2:483066/1‑140 (MQ=255)                                                                                                    ttGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGcgtc                                   >  2:1334909/1‑140 (MQ=255)                                                                                                     tGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGcgtcg                                  >  2:278728/1‑140 (MQ=255)                                                                                                     tGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGcgtcg                                  >  1:1222881/1‑140 (MQ=255)                                                                                                       tCAATAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGtt                                <  2:144822/140‑1 (MQ=255)                                                                                                         aaTAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGtttt                              >  2:832303/1‑140 (MQ=255)                                                                                                         aaTAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGtttt                              >  2:900216/1‑140 (MQ=255)                                                                                                           tAACCGGGCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCa                            >  2:734744/1‑140 (MQ=255)                                                                                                                 ggCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGa                      >  1:902178/1‑140 (MQ=255)                                                                                                                 ggCCGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGa                      >  1:1386483/1‑140 (MQ=255)                                                                                                                   ccGATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATc                    >  1:713043/1‑140 (MQ=255)                                                                                                                     gATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGc                  <  2:1222881/140‑1 (MQ=255)                                                                                                                     gATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGc                  <  1:429809/140‑1 (MQ=255)                                                                                                                     gATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGc                  <  1:344055/140‑1 (MQ=255)                                                                                                                     gATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGc                  <  1:1014794/140‑1 (MQ=255)                                                                                                                     gATTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGAGCGTCGTTTTCAATCAGATCGc                  <  1:1149240/140‑1 (MQ=255)                                                                                                                      aTTGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCa                 >  2:781484/1‑140 (MQ=255)                                                                                                                       ttGCCATTTTGGTATGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCaa                <  1:900216/140‑1 (MQ=255)                                                                                                                                     tGCGCCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTc  >  2:177377/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        TCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTC  >  NC_000913/1802173‑1802444

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