Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,410,733	G→A	E152K (GAA→AAA)	dusB →	tRNA‑dihydrouridine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,410,733	0	G	A	100.0%	194.5 / NA	56	E152K (GAA→AAA)	dusB	tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/28);  total (28/28)

TCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGG  >  NC_000913/3410597‑3410868                                                                                                                                         |                                                                                                                                        tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAag                                                                                                                                      >  2:333730/1‑140 (MQ=255) tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAag                                                                                                                                      >  1:626275/1‑140 (MQ=255)            cAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCaa                                                                                                                           >  1:321358/1‑140 (MQ=255)            cAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCaa                                                                                                                           >  1:376046/1‑140 (MQ=255)             aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAAc                                                                                                                          <  2:321358/140‑1 (MQ=255)             aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAAc                                                                                                                          <  2:587506/140‑1 (MQ=255)             aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAAc                                                                                                                          <  2:376046/140‑1 (MQ=255)             aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAAc                                                                                                                          <  1:333730/140‑1 (MQ=255)             aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAAc                                                                                                                          <  2:376504/140‑1 (MQ=255)             aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAAc                                                                                                                          <  2:109323/140‑1 (MQ=255)                          ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACtg                                                                                                             >  1:683828/1‑140 (MQ=255)                               gTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAAc                                                                                                                          <  1:64183/122‑1 (MQ=255)                               gTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAAc                                                                                                                          >  2:64183/1‑122 (MQ=255)                                 aCCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                                                      <  1:951286/140‑1 (MQ=255)                                 aCCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                                                      <  1:77963/140‑1 (MQ=255)                                  cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTc                                                                                                     <  2:596286/140‑1 (MQ=255)                                     ggATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGAAAAGAGAGTGCCCAACTGGCTGAAGACTGTGGCAGTCAgg                                                                                                  >  2:1082686/1‑140 (MQ=255)                                      gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGc                                                                                                 <  1:74052/140‑1 (MQ=255)                                           cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGa                                                                                            <  2:683828/140‑1 (MQ=255)                                                 aTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                                                      <  2:1066045/124‑1 (MQ=255)                                                 aTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                                                      >  1:1066045/1‑124 (MQ=255)                                                 aTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTc                                                                                      >  2:865545/1‑140 (MQ=255)                                                     aTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATgg                                                                                  >  2:756732/1‑140 (MQ=255)                                                      tCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGc                                                                                 >  2:276907/1‑140 (MQ=255)                                                          tACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTa                                                                             <  2:626275/140‑1 (MQ=255)                                                                ggTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAcgcg                                                                       <  1:740806/140‑1 (MQ=255)                                                                ggTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAcgcg                                                                       >  1:665357/1‑140 (MQ=255)                                                                ggTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAcgcg                                                                       >  1:36158/1‑140 (MQ=255)                                                                ggTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAcgcg                                                                       >  2:137305/1‑140 (MQ=255)                                                                 gtcgtcAATGCAGTGGACGTTGCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACgcgc                                                                      >  2:869526/1‑140 (MQ=255)                                                                 gtcgtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACgcgc                                                                      >  2:674504/1‑140 (MQ=255)                                                                    gtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtg                                                                   >  2:497507/1‑140 (MQ=255)                                                                     tcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgt                                                                  <  1:251684/140‑1 (MQ=255)                                                                     tcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgt                                                                  <  2:175988/140‑1 (MQ=255)                                                                      cAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtt                                                                 >  1:757338/1‑140 (MQ=255)                                                                        aTGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttg                                                               <  1:263589/140‑1 (MQ=255)                                                                             gtagACGTTCCATTTACCCTGAGGAGTCGCTCCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAACACTGTCGCATTCAGGCTCTGACCATTCATGCCCGTACACGCGCCTGTTTGTTCaa                                                          <  2:274413/137‑1 (MQ=255)                                                                                     tCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAg                                                  <  1:756732/140‑1 (MQ=255)                                                                                       cTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCt                                                <  1:723409/140‑1 (MQ=255)                                                                                          ttACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAg                                                  <  2:987935/135‑1 (MQ=255)                                                                                          ttACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAg                                                  >  1:987935/1‑135 (MQ=255)                                                                                           tACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGt                                            <  1:674504/140‑1 (MQ=255)                                                                                            aCCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                                           <  1:865545/140‑1 (MQ=255)                                                                                               cTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGa                                        <  2:757338/140‑1 (MQ=255)                                                                                               cTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGa                                        <  1:911803/140‑1 (MQ=255)                                                                                               cTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGa                                        >  1:704555/1‑140 (MQ=255)                                                                                                     aTTCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAt                                  >  1:831075/1‑140 (MQ=255)                                                                                                       tCGCACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTc                                >  1:404033/1‑140 (MQ=255)                                                                                                          cACCGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCggg                             <  2:1196366/140‑1 (MQ=255)                                                                                                            ccGGCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCa                           >  1:197951/1‑140 (MQ=255)                                                                                                               gCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGtt                        >  2:8281/1‑140 (MQ=255)                                                                                                               gCTGGGCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGtt                        <  1:1082686/140‑1 (MQ=255)                                                                                                                   ggCACCGGAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGc                    <  1:137305/140‑1 (MQ=255)                                                                                                                          gAACACCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTATGACAGTATTCGGGCAGTTAAGCAGAAAGt             >  2:791006/1‑140 (MQ=255)                                                                                                                             cacCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTc          >  2:1127136/1‑140 (MQ=255)                                                                                                                              acCGTAACTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTcc         >  1:890814/1‑140 (MQ=255)                                                                                                                                     cTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCgg  >  2:1113739/1‑140 (MQ=255)                                                                                                                                     cTGCAAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCgg  >  2:242465/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        TCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGG  >  NC_000913/3410597‑3410868

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