Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,410,406	T→C	S43P (TCT→CCT)	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,406	0	T	C	100.0%	235.6 / NA	65	S43P (TCT→CCT)	dusB	tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/28);  total (37/28)

ATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACG  >  NC_000913/3410280‑3410538                                                                                                                               |                                                                                                                                     aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                          <  2:87752/139‑1 (MQ=255) aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                          >  1:87752/1‑139 (MQ=255) aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                          >  1:1392493/1‑139 (MQ=255) aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                          <  2:1392493/139‑1 (MQ=255)                     gCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAAccc                                                                                                                            <  1:435545/117‑1 (MQ=255)                     gCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAAccc                                                                                                                            >  2:435545/1‑117 (MQ=255)                      cTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAAccc                                                                                                                            >  2:1409333/1‑116 (MQ=255)                      cTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAAccc                                                                                                                            <  1:1409333/116‑1 (MQ=255)                      cTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAt                                                                                                     >  1:1326362/1‑139 (MQ=255)                        cAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                                   >  2:130917/1‑139 (MQ=255)                          gAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCg                                                                                                 >  2:1395209/1‑139 (MQ=255)                           aaaTCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGt                                                                                                >  2:1645922/1‑139 (MQ=255)                                gCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGccttct                                                                                                                                 >  1:1079916/1‑101 (MQ=255)                                gCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGccttct                                                                                                                                 <  2:1079916/101‑1 (MQ=255)                                  cTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                         <  2:1326362/139‑1 (MQ=255)                                  cTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                         <  1:210403/139‑1 (MQ=255)                                    gATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATg                                                                                       <  1:1325061/139‑1 (MQ=255)                                    gATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATg                                                                                       <  1:761538/139‑1 (MQ=255)                                     aTCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAt                                                                                                     <  1:1649697/124‑1 (MQ=255)                                     aTCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGCTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAt                                                                                                     >  2:1649697/1‑124 (MQ=255)                                                     ctggcATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAcc                                                                      <  1:1395209/139‑1 (MQ=255)                                                      tggcATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                     <  1:1178978/139‑1 (MQ=255)                                                                   agacCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACA‑TGATGAACCCGGTATTCGCACCGt                                                       >  2:952824/1‑139 (MQ=255)                                                                     acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                          >  1:1212935/1‑71 (MQ=255)                                                                     acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                          <  2:1212935/71‑1 (MQ=255)                                                                      ccTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGc                                                     >  1:1402020/1‑139 (MQ=255)                                                                          ttCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAAt                                                 >  2:1210318/1‑139 (MQ=255)                                                                             ggACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGc                                              <  2:459559/139‑1 (MQ=255)                                                                                   tgtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAAGCCACAGGTTTGGGAAAGCGCCAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                        <  1:1537534/139‑1 (MQ=255)                                                                                   tgtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                        >  1:956139/1‑139 (MQ=255)                                                                                    gtgtTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGc                                       >  1:1509339/1‑139 (MQ=255)                                                                                    gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCg                                                                                                 <  1:1620037/81‑1 (MQ=255)                                                                                    gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCg                                                                                                 >  2:1620037/1‑81 (MQ=255)                                                                                    gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                     >  2:1537440/1‑109 (MQ=255)                                                                                    gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                     <  1:1537440/109‑1 (MQ=255)                                                                                    gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                        <  2:1525761/138‑1 (MQ=255)                                                                                    gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                        >  1:1525761/1‑138 (MQ=255)                                                                                       cTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGAt                                    <  1:1645922/139‑1 (MQ=255)                                                                                        tACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATc                                   >  1:494108/1‑139 (MQ=255)                                                                                           gagaTGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCga                                >  2:162902/1‑139 (MQ=255)                                                                                            agaTGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCaaa                                                  <  1:139468/120‑1 (MQ=255)                                                                                            agaTGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCaaa                                                  >  2:139468/1‑120 (MQ=255)                                                                                            agaTGGGAGCCAGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaa                               >  1:1693808/1‑139 (MQ=255)                                                                                             gaTGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaaa                              >  2:245294/1‑139 (MQ=255)                                                                                               tGGGAGCCGGATTGACAGTATCCGAGATGATGccttc                                                                                                                                  <  1:360471/37‑1 (MQ=255)                                                                                               tGGGAGCCGGATTGACAGTATCCGAGATGATGccttc                                                                                                                                  >  2:360471/1‑37 (MQ=255)                                                                                               tGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaaaga                            >  2:1235862/1‑139 (MQ=255)                                                                                               tGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCCCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaaaga                            >  2:438020/1‑139 (MQ=255)                                                                                                  gAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                     <  2:63625/95‑1 (MQ=255)                                                                                                  gAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                     >  1:63625/1‑95 (MQ=255)                                                                                                     ccGGATTGACAGTATCCGAGATGATGCCTTCTAACCca                                                                                                                           >  2:1408721/1‑38 (MQ=255)                                                                                                     ccGGATTGACAGTATCCGAGATGATGCCTTCTAACCca                                                                                                                           <  1:1408721/38‑1 (MQ=255)                                                                                                        gATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGa                   >  2:1716639/1‑139 (MQ=255)                                                                                                             aCAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                        >  1:184448/1‑113 (MQ=255)                                                                                                             aCAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                        <  2:184448/113‑1 (MQ=255)                                                                                                             aCAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATgcag              <  1:1210318/139‑1 (MQ=255)                                                                                                                   tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                     >  1:663009/1‑78 (MQ=255)                                                                                                                   tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                                     <  2:663009/78‑1 (MQ=255)                                                                                                                   tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                        <  2:1636163/107‑1 (MQ=255)                                                                                                                   tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAg                                        >  1:1636163/1‑107 (MQ=255)                                                                                                                   tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATgcagc             <  1:462228/134‑1 (MQ=255)                                                                                                                   tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATgaagc             >  2:462228/1‑134 (MQ=255)                                                                                                                        gatgatGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGa                   <  1:798771/123‑1 (MQ=255)                                                                                                                        gatgatGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGa                   >  2:798771/1‑123 (MQ=255)                                                                                                                        gatgatGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAAc   >  2:423784/1‑139 (MQ=255)                                                                                                                         atgatgCCTTCTAACCCACAGGTTTGGGAAAGCGa                                                                                                          >  2:1619021/1‑35 (MQ=255)                                                                                                                         atgatgCCTTCTAACCCACAGGTTTGGGAAAGCGa                                                                                                          <  1:1619021/35‑1 (MQ=255)                                                                                                                         atgatgCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACg  >  2:1616247/1‑139 (MQ=255)                                                                                                                            atgCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATgca               >  2:114937/1‑123 (MQ=255)                                                                                                                               |                                                                                                                                     ATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACG  >  NC_000913/3410280‑3410538

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