Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,411,139	T→C	L287P (CTC→CCC)	dusB →	tRNA‑dihydrouridine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,411,139	0	T	C	100.0%	182.9 / NA	49	L287P (CTC→CCC)	dusB	tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/21);  total (28/21)

ACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCAGAGGTTAAGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCTCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATAAAGAGCTGACAGAAC  >  NC_000913/3411004‑3411269                                                                                                                                        |                                                                                                                                   acacTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCAGAGGTTAAGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCa                                                                                                                                 >  1:797643/1‑139 (MQ=255)     tGGGGAGTTGCTGCCCCCGCTGCCTTTGGCAGAGGTTAAGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGaa                                                                                                                             >  1:1098954/1‑139 (MQ=255)                       gCTGCCTTTGGCAGAGGTTAAGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAg                                                                                                           >  1:1221579/1‑139 (MQ=255)                        cTGCCTTTGGCAGAGGTTAAGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCa                                                                                                                                 >  1:804370/1‑116 (MQ=255)                        cTGCCTTTGGCAGAGGTTAAGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGCATCCCCa                                                                                                                                 <  2:804370/116‑1 (MQ=255)                                        ttAAGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCaa                                                                                          <  2:1098954/139‑1 (MQ=255)                                          aaGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACg                                                                                        >  2:981811/1‑139 (MQ=255)                                           aGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGca                                                                                                 <  2:928499/129‑1 (MQ=255)                                           aGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATAGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGca                                                                                                 >  1:928499/1‑129 (MQ=255)                                               cttgcttTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAac                                                                                                                            >  1:1302581/1‑98 (MQ=255)                                               cttgcttTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAac                                                                                                                            <  2:1302581/98‑1 (MQ=255)                                                 tgcttTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCaaa                                                                                                                  >  1:77067/1‑106 (MQ=255)                                                 tgcttTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCaaa                                                                                                                  <  2:77067/106‑1 (MQ=255)                                                 tgcttTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGa                                                                                 >  2:1136281/1‑139 (MQ=255)                                                                ttCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACa                                                                  >  1:676007/1‑139 (MQ=255)                                                                    ggAACTGCATGACTTTTATGGTCCGGCAAAAGGGTTCCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGGGGATGCCAGCGAACAGCTg                                                              <  2:877589/139‑1 (MQ=255)                                                                     gAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTgg                                                             <  2:797643/139‑1 (MQ=255)                                                                       aCTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAg                                                           <  1:996528/139‑1 (MQ=255)                                                                           caggtcTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGt                                                       >  1:764729/6‑139 (MQ=255)                                                                           cATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGca                                                                                                 >  1:686292/1‑97 (MQ=255)                                                                           cATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGca                                                                                                 <  2:686292/97‑1 (MQ=255)                                                                           cATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGt                                                       >  1:526591/1‑139 (MQ=255)                                                                           cATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGt                                                       >  2:158349/1‑139 (MQ=255)                                                                           cATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGt                                                       >  2:956017/1‑139 (MQ=255)                                                                           cATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGt                                                       >  1:1185838/1‑139 (MQ=255)                                                                             tGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTg                                                     >  1:905970/1‑139 (MQ=255)                                                                                   ttATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCa                                               <  2:840679/139‑1 (MQ=255)                                                                                     aTGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATa                                             >  1:494305/1‑139 (MQ=255)                                                                                          ccGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCg                                        >  2:851574/1‑139 (MQ=255)                                                                                                       aCCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTa                           <  1:956017/139‑1 (MQ=255)                                                                                                       aCCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTa                           <  2:233070/139‑1 (MQ=255)                                                                                                       aCCGAATTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTa                           <  1:1136281/139‑1 (MQ=255)                                                                                                            aTTGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATAGAGGATGCCAGTGAACAGCTGGAGGCGTTGGATGCATACTTCGAAAATTTTGCGTAAACAg                      <  1:851574/139‑1 (MQ=255)                                                                                                             ttGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGa                     >  1:903143/1‑139 (MQ=255)                                                                                                              tGCACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGATAATTTTGCGTAAACAGaa                    >  2:128008/1‑139 (MQ=255)                                                                                                                cACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCaa                                                                                          <  2:1023247/67‑1 (MQ=255)                                                                                                                cACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCaa                                                                                          >  1:1023247/1‑67 (MQ=255)                                                                                                                cACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCaa                                                                                          <  1:776625/67‑1 (MQ=255)                                                                                                                cACGTAAACACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCaa                                                                                          >  2:776625/1‑67 (MQ=255)                                                                                                                  cGTAAGCACGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATaa                >  2:994213/1‑139 (MQ=255)                                                                                                                  cGTAAACACGTTTCCTGGTATCCCCAGGAGCACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATaa                >  2:638359/1‑139 (MQ=255)                                                                                                                       acacGTTTCCTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATAAAGAGc           <  2:905970/139‑1 (MQ=255)                                                                                                                               ccTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATAAAGAGCTGa        >  2:680340/1‑134 (MQ=255)                                                                                                                               ccTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATAAAGAGCTGa        <  1:680340/134‑1 (MQ=255)                                                                                                                               ccTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATAAAGAGCTGACAGaa   >  2:161171/1‑139 (MQ=255)                                                                                                                                cTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATAAAGAGCTGACAGAAc  <  2:75282/139‑1 (MQ=255)                                                                                                                                cTGGTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATAAAGAGCTGACAGAAc  <  2:399591/139‑1 (MQ=255)                                                                                                                                   gTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGt                                                       >  2:674894/1‑83 (MQ=255)                                                                                                                                   gTATCCCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGt                                                       <  1:674894/83‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   ACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCAGAGGTTAAGCGCTTGCTTTGCGCGCACGTTCGGGAACTGCATGACTTTTATGGTCCGGCAAAAGGGTACCGAATTGCACGTAAACACGTTTCCTGGTATCTCCAGGAACACGCTCCAAATGACCAGTTTCGGCGCACATTCAACGCCATTGAGGATGCCAGCGAACAGCTGGAGGCGTTGGAGGCATACTTCGAAAATTTTGCGTAAACAGAAATAAAGAGCTGACAGAAC  >  NC_000913/3411004‑3411269

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