Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,261,710	C→A	S14R (AGC→AGA)	dusA →	tRNA‑dihydrouridine synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,710	0	C	A	100.0%	220.8 / NA	62	S14R (AGC→AGA)	dusA	tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/30);  total (32/30)

ATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261581‑4261844                                                                                                                                  |                                                                                                                                       ataGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTg                                                                                                                              >  1:262704/1‑140 (MQ=255)  taGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGa                                                                                                                             >  1:453416/1‑140 (MQ=255)    gTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGaaa                                                                                                                           >  1:23277/1‑140 (MQ=255)      aCCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGaaaaa                                                                                                                         >  2:490750/1‑140 (MQ=255)       ccccGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGaaaaaa                                                                                                                        >  1:552641/1‑140 (MQ=255)          cGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTg                                                                                                                     >  1:615542/1‑140 (MQ=255)          cGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTg                                                                                                                     >  2:663154/1‑140 (MQ=255)            cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATTGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   <  1:189648/140‑1 (MQ=255)            cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   <  2:181470/140‑1 (MQ=255)            cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   <  1:245123/140‑1 (MQ=255)            cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   <  1:16336/140‑1 (MQ=255)            cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   >  1:123683/1‑140 (MQ=255)            cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   <  1:594061/140‑1 (MQ=255)             tATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACg                                                                                                                  >  1:173658/1‑140 (MQ=255)               ttGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGtt                                                                                                                <  1:171197/140‑1 (MQ=255)                    cctcctGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTg                                                                                                           <  2:557773/140‑1 (MQ=255)                      tcctGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGa                                                                                                         <  1:266369/140‑1 (MQ=255)                          gAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTgg                                                                                                     <  2:262704/140‑1 (MQ=255)                           aaTAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCATGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGt                                                                                                    <  2:514992/140‑1 (MQ=255)                                 gATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGttt                                                                                              <  2:611306/140‑1 (MQ=255)                                 gATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGttt                                                                                              <  2:401725/140‑1 (MQ=255)                                  aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  2:100522/1‑140 (MQ=255)                                  aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  1:339919/1‑140 (MQ=255)                                  aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  1:450241/1‑140 (MQ=255)                                  aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  1:476922/1‑140 (MQ=255)                                  aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  2:419162/1‑140 (MQ=255)                                   tGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAg                                                                                            <  1:252466/140‑1 (MQ=255)                                          taacataacCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAg                                                                                            <  1:508418/133‑1 (MQ=255)                                          taacataacCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAg                                                                                            >  2:508418/1‑133 (MQ=255)                                                  ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  2:104761/1‑140 (MQ=255)                                                  ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  2:555566/1‑140 (MQ=255)                                                  ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  1:651126/1‑140 (MQ=255)                                                    ccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCg                                                                           >  2:69978/1‑140 (MQ=255)                                                     cATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGa                                                                          >  2:306855/1‑140 (MQ=255)                                                               cGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacggac                                                                >  1:307764/1‑140 (MQ=255)                                                               cGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacggac                                                                >  1:289981/1‑140 (MQ=255)                                                               cGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacggac                                                                >  2:360580/1‑140 (MQ=255)                                                                 aTAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTg                                                                       <  2:123024/131‑1 (MQ=255)                                                                 aTAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTg                                                                       >  1:123024/1‑131 (MQ=255)                                                                 aTAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacag                                                              >  1:194065/1‑140 (MQ=255)                                                                 aTAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacag                                                              >  2:68647/1‑140 (MQ=255)                                                                 aTAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacag                                                              >  2:370332/1‑140 (MQ=255)                                                                   aaaaTAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacagac                                                            >  1:426177/1‑140 (MQ=255)                                                                           ctctACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGc                                                    <  1:663154/140‑1 (MQ=255)                                                                            tctACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCt                                                   >  1:234871/1‑140 (MQ=255)                                                                                  tttGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCt                                             <  2:217702/140‑1 (MQ=255)                                                                                          tGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTg                                     >  1:249476/1‑140 (MQ=255)                                                                                                    ttCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCa                           >  1:321044/1‑140 (MQ=255)                                                                                                     tCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCaa                          <  1:89585/140‑1 (MQ=255)                                                                                                         aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACg                      <  1:100522/140‑1 (MQ=255)                                                                                                         aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACg                      <  1:490750/140‑1 (MQ=255)                                                                                                           aTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGtt                    <  1:117719/140‑1 (MQ=255)                                                                                                           aTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGtt                    >  2:227126/1‑140 (MQ=255)                                                                                                            tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  2:307764/140‑1 (MQ=255)                                                                                                            tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  2:234871/140‑1 (MQ=255)                                                                                                            tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  2:23277/140‑1 (MQ=255)                                                                                                                 aaaaTCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGta              <  2:354671/140‑1 (MQ=255)                                                                                                                        accaaaccaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGaa       <  1:104761/140‑1 (MQ=255)                                                                                                                        accaaaccaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGaa       <  1:68647/140‑1 (MQ=255)                                                                                                                             accaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGt  <  2:615542/140‑1 (MQ=255)                                                                                                                             accaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGt  <  2:37319/140‑1 (MQ=255)                                                                                                                             accaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGt  <  2:357417/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                       ATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261581‑4261844

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