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%	206.2 / NA	58	S43P (TCT→CCT)	dusB	tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/27);  total (31/27)

TCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTA  >  NC_000913/3410271‑3410532                                                                                                                                        |                                                                                                                               tCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGcctt                                                                                                                             >  1:1127732/1‑139 (MQ=255) tCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGcctt                                                                                                                             >  2:277561/1‑139 (MQ=255)          aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                    >  2:151204/1‑139 (MQ=255)          aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                    <  2:1311480/139‑1 (MQ=255)          aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                    >  2:385961/1‑139 (MQ=255)          aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                    >  2:815625/1‑139 (MQ=255)          aTGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                    >  1:170534/1‑139 (MQ=255)              gcATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGt                                                                                                                <  1:1360504/139‑1 (MQ=255)               cATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGtt                                                                                                               >  1:1611097/1‑139 (MQ=255)                              gCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                    <  1:59895/119‑1 (MQ=255)                              gCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCcac                                                                                                                    >  2:59895/1‑119 (MQ=255)                               cTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAt                                                                                               >  1:1271906/1‑139 (MQ=255)                               cTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAt                                                                                               <  2:749927/139‑1 (MQ=255)                                tCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtc                                                                                              <  1:815625/139‑1 (MQ=255)                                 cAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                             >  1:582742/1‑139 (MQ=255)                                          ccTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACgg                                                                                    <  1:909602/139‑1 (MQ=255)                                           cTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                   <  2:1611097/139‑1 (MQ=255)                                               tCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGt                                                                               <  2:304969/139‑1 (MQ=255)                                                  cAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGca                                                                            <  2:582742/139‑1 (MQ=255)                                                     cgcCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCca                                                                                                                     <  1:1515565/95‑1 (MQ=255)                                                     cgcCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCca                                                                                                                     >  2:1515565/1‑95 (MQ=255)                                                          aTGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATtgatg                                                                    <  1:480977/139‑1 (MQ=255)                                                             gctggcATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCAGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAc                                                                 <  1:85181/139‑1 (MQ=255)                                                               tggcATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAAccc                                                               <  2:473974/139‑1 (MQ=255)                                                                 gcATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGATTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCgg                                                             >  2:530861/1‑139 (MQ=255)                                                                          acagacCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCAc                                                    <  1:151204/139‑1 (MQ=255)                                                                              acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGccttc                                                                                                                            <  1:537089/63‑1 (MQ=255)                                                                              acCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGccttc                                                                                                                            >  2:537089/1‑63 (MQ=255)                                                                                 ttttCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACgg                                                                                    <  1:1452862/100‑1 (MQ=255)                                                                                 ttttCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACgg                                                                                    >  2:1452862/1‑100 (MQ=255)                                                                                  tttCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCca                                                                                                                     >  1:1476096/1‑66 (MQ=255)                                                                                  tttCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCca                                                                                                                     <  2:1476096/66‑1 (MQ=255)                                                                                      ggACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGc                                        <  1:1154188/139‑1 (MQ=255)                                                                                             gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCca                                                                                                                     <  1:612224/55‑1 (MQ=255)                                                                                             gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCca                                                                                                                     >  2:612224/1‑55 (MQ=255)                                                                                             gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGc                                 >  1:1392840/1‑139 (MQ=255)                                                                                             gtgCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGc                                 >  2:897821/1‑139 (MQ=255)                                                                                                      gaTGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaaa                        >  1:1242415/1‑139 (MQ=255)                                                                                                        tGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                   >  1:865298/1‑78 (MQ=255)                                                                                                        tGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGa                                                                                   <  2:865298/78‑1 (MQ=255)                                                                                                        tGGGAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaaaga                      >  1:1104139/1‑139 (MQ=255)                                                                                                          ggAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGca                                                                            >  1:539087/1‑83 (MQ=255)                                                                                                          ggAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGca                                                                            <  2:874688/83‑1 (MQ=255)                                                                                                          ggAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGca                                                                            <  2:539087/83‑1 (MQ=255)                                                                                                          ggAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGca                                                                            >  1:874688/1‑83 (MQ=255)                                                                                                          ggAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaaagaaa                    >  2:1332948/1‑139 (MQ=255)                                                                                                          ggAGCCGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCgaaagaaa                    >  2:1275312/1‑139 (MQ=255)                                                                                                              ccGGATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGc                >  1:382284/1‑139 (MQ=255)                                                                                                                 gATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                             >  1:331557/1‑59 (MQ=255)                                                                                                                 gATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                             <  2:331557/59‑1 (MQ=255)                                                                                                                 gATTGACAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGa             <  1:385961/139‑1 (MQ=255)                                                                                                                      aCAGTATCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATgcag        >  2:557531/1‑139 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                             >  1:776257/1‑48 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAAtct                                                                                             <  2:776257/48‑1 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACgg                                                                                    >  2:512991/1‑57 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACgg                                                                                    <  1:512991/57‑1 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGt                                                                               >  1:12826/1‑62 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGt                                                                               <  2:12826/62‑1 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTg                                                                              >  1:1078486/1‑63 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTg                                                                              <  2:1078486/63‑1 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTa  >  2:276361/1‑139 (MQ=255)                                                                                                                            tCCGAGATGATGCCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACAGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCCCGTa  >  2:1148576/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                               TCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGATTGACAGTATCCGAGATGATGTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTA  >  NC_000913/3410271‑3410532

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