Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,410,526	G→A	A83T (GCA→ACA)	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,526	0	G	A	100.0%	177.7 / NA	51	A83T (GCA→ACA)	dusB	tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/25);  total (26/25)

AGATGATGTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGT  >  NC_000913/3410398‑3410662                                                                                                                                 |                                                                                                                                         agaTGATGTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAAc                                                                                                                               <  1:591944/140‑1 (MQ=255)       tgTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGaa                                                                                                                         >  1:733905/1‑140 (MQ=255)        gTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGaaa                                                                                                                        >  1:393682/1‑140 (MQ=255)          cttctAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGc                                                                                                                      <  1:563214/140‑1 (MQ=255)           ttctAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCg                                                                                                                     <  2:581873/140‑1 (MQ=255)           ttctAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCg                                                                                                                     <  2:171937/140‑1 (MQ=255)            tctAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCgg                                                                                                                    >  1:48893/1‑140 (MQ=255)              tAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTg                                                                                                                  >  1:189359/1‑140 (MQ=255)                aCCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGcc                                                                                                                >  1:323726/1‑140 (MQ=255)                   cacaGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGTGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGc                                                                                                                 >  2:204178/1‑136 (MQ=255)                   cacaGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGTGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGc                                                                                                                 <  1:204178/136‑1 (MQ=255)                   cacaGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGcc                                                                                                                >  2:683246/1‑137 (MQ=255)                   cacaGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGcc                                                                                                                <  1:683246/137‑1 (MQ=255)                    acaGGCTTGGGAAAGCGACAAATCTCCTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGa                                                                                                            <  2:733905/140‑1 (MQ=255)                     caGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGc                                                                                                                      <  2:30858/129‑1 (MQ=255)                     caGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGc                                                                                                                      >  1:30858/1‑129 (MQ=255)                       ggTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGatta                                                                                                         >  1:80916/1‑140 (MQ=255)                            gggAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAATCGGTGCCCAGATTATTGat                                                                                                    <  2:540631/140‑1 (MQ=255)                            gggAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGat                                                                                                    <  2:492273/140‑1 (MQ=255)                            gggAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGat                                                                                                    <  2:323726/140‑1 (MQ=255)                                  gCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGat                                                                                                    <  1:421057/134‑1 (MQ=255)                                  gCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGat                                                                                                    >  2:421057/1‑134 (MQ=255)                                  gCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAat                                                                                              >  2:542941/1‑140 (MQ=255)                                             cGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGccc                                                                                   <  2:189359/140‑1 (MQ=255)                                                    ggATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTaaa                                                                            >  2:116582/1‑140 (MQ=255)                                                                 tgatgaACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCa                                                               >  1:371952/1‑140 (MQ=255)                                                                   atgaACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAg                                                             >  2:523633/1‑140 (MQ=255)                                                                        cccGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGc                                                        <  1:116582/140‑1 (MQ=255)                                                                           ggTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAgg                                                     >  2:52169/1‑140 (MQ=255)                                                                           ggTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAgg                                                     >  2:108623/1‑140 (MQ=255)                                                                             tATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGtt                                                   >  1:384436/1‑140 (MQ=255)                                                                                tCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAg                                                <  1:503060/140‑1 (MQ=255)                                                                                tCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAg                                                <  2:393682/140‑1 (MQ=255)                                                                                   cACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGccc                                             >  2:393380/1‑140 (MQ=255)                                                                                   cACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGccc                                             >  1:263067/1‑140 (MQ=255)                                                                                        tGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTg                                        <  1:393380/140‑1 (MQ=255)                                                                                         gCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGc                                       >  1:145438/1‑140 (MQ=255)                                                                                                   ggTTGCGATCCGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGa                             <  2:499359/140‑1 (MQ=255)                                                                                                             cGAAAGAAATGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTaaa                   <  2:80916/140‑1 (MQ=255)                                                                                                              gaaagaaaTGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGa                             <  2:328172/129‑1 (MQ=255)                                                                                                              gaaagaaaTGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGa                             >  1:328172/1‑129 (MQ=255)                                                                                                              gaaagaaaTGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAAt                  >  1:80917/1‑140 (MQ=255)                                                                                                              gaaagaaaTGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGACCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAAt                  <  2:38603/140‑1 (MQ=255)                                                                                                              gaaagaaaTGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGACCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAAt                  >  1:38603/1‑140 (MQ=255)                                                                                                                aagaaaTGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGa                             >  2:354422/1‑127 (MQ=255)                                                                                                                aagaaaTGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGa                             <  1:354422/127‑1 (MQ=255)                                                                                                                     aTGGCAGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCt           >  1:89233/1‑140 (MQ=255)                                                                                                                          aGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCg      <  2:145438/140‑1 (MQ=255)                                                                                                                          aGATGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCg      <  1:52169/140‑1 (MQ=255)                                                                                                                            aTGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAg    <  2:263067/140‑1 (MQ=255)                                                                                                                            aTGCAACACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAg    >  1:746390/1‑140 (MQ=255)                                                                                                                              gcaacaCGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGgt  <  2:746390/140‑1 (MQ=255)                                                                                                                              gcaacaCGTATTAACGTGGAAAGCGGTGCCCAGATTAATGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGgt  <  1:779801/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                         AGATGATGTCTTCTAACCCACAGGTTTGGGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGT  >  NC_000913/3410398‑3410662

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