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%	197.0 / NA	57	S14R (AGC→AGA)	dusA	tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (36/21);  total (36/21)

TGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261577‑4261844                                                                                                                                      |                                                                                                                                       tGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATg                                                                                                                                  >  2:81768/1‑140 (MQ=255)    tataGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCt                                                                                                                               >  1:200771/1‑140 (MQ=255)        gTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGaaa                                                                                                                           >  1:269648/1‑140 (MQ=255)           ccccGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGaaaaaa                                                                                                                        >  2:256368/1‑140 (MQ=255)             ccGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACt                                                                                                                      <  1:239680/140‑1 (MQ=255)              cGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTg                                                                                                                     >  1:412792/1‑140 (MQ=255)              cGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTg                                                                                                                     >  1:315272/1‑140 (MQ=255)              cGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTg                                                                                                                     >  1:300105/1‑140 (MQ=255)                cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   >  2:270292/1‑140 (MQ=255)                cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   <  1:175943/140‑1 (MQ=255)                     gAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCa                                                                                                              <  2:243298/140‑1 (MQ=255)                      aGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCAc                                                                                                             >  2:300553/1‑140 (MQ=255)                          tcctGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGa                                                                                                         >  1:23460/1‑140 (MQ=255)                           cctGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAg                                                                                                        >  2:254252/1‑140 (MQ=255)                           cctGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAg                                                                                                        >  1:413902/1‑140 (MQ=255)                                  aGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCgg                                                                                                 >  1:117856/1‑140 (MQ=255)                                    tGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGtt                                                                                               >  2:115906/1‑140 (MQ=255)                                      aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  1:290841/1‑140 (MQ=255)                                      aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  2:29214/1‑140 (MQ=255)                                       tGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCACAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAg                                                                                            <  1:284434/140‑1 (MQ=255)                                         cTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCg                                                                                          >  1:452348/1‑140 (MQ=255)                                              taacataacCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCa                                                                                     >  1:123014/1‑140 (MQ=255)                                               aacataacCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCAc                                                                                    <  2:427765/140‑1 (MQ=255)                                                      ccccATGATATATCGATAAAATAATCTCTACATTTTAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  2:229669/1‑140 (MQ=255)                                                      ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  2:53711/1‑140 (MQ=255)                                                      ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  2:220044/1‑140 (MQ=255)                                                            gATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCAc                                                                                    <  1:337440/127‑1 (MQ=255)                                                            gATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCAc                                                                                    >  2:337440/1‑127 (MQ=255)                                                              tatatCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTgga                                                                     <  1:290834/140‑1 (MQ=255)                                                                tatCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacg                                                                   >  2:379957/1‑140 (MQ=255)                                                                 atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacgg                                                                  >  2:308883/1‑140 (MQ=255)                                                                 atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTgcacgg                                                                  >  1:409163/1‑140 (MQ=255)                                                                          ataatCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACAtt                                                         >  2:356359/1‑140 (MQ=255)                                                                             atCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGcc                                                      <  2:452348/140‑1 (MQ=255)                                                                              tctcTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCg                                                     >  2:378026/1‑140 (MQ=255)                                                                                    cATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTAttt                                               >  1:154954/1‑140 (MQ=255)                                                                                    cATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTAttt                                               >  1:416406/1‑140 (MQ=255)                                                                                       ttGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCtt                                            >  2:195445/1‑140 (MQ=255)                                                                                         gAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGc                                          >  2:77916/1‑140 (MQ=255)                                                                                          aaaaTGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCg                                         <  2:154954/140‑1 (MQ=255)                                                                                          aaaaTGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCg                                         <  2:23460/140‑1 (MQ=255)                                                                                                 aCGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCtt                                  <  1:308883/140‑1 (MQ=255)                                                                                                      aaTTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCg                             <  1:53711/140‑1 (MQ=255)                                                                                                        ttCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCa                           >  1:272658/1‑140 (MQ=255)                                                                                                           tGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATa                        <  1:81768/140‑1 (MQ=255)                                                                                                           tGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATa                        <  1:378026/140‑1 (MQ=255)                                                                                                             aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTc                                >  1:8566/1‑130 (MQ=255)                                                                                                             aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTc                                <  2:8566/130‑1 (MQ=255)                                                                                                             aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACg                      <  1:115906/140‑1 (MQ=255)                                                                                                             aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACg                      <  1:29214/140‑1 (MQ=255)                                                                                                               aTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGtt                    >  1:143502/1‑140 (MQ=255)                                                                                                                tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  1:125367/140‑1 (MQ=255)                                                                                                                 gCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   >  2:265296/1‑139 (MQ=255)                                                                                                                 gCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  1:265296/139‑1 (MQ=255)                                                                                                                  cAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGGGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCt                 >  2:362258/1‑140 (MQ=255)                                                                                                                         tCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATAcc          <  2:200771/140‑1 (MQ=255)                                                                                                                                 accaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGt  <  2:168164/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       TGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261577‑4261844

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: ‑