Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,998,412	T→C	W144R (TGG→CGG)	uvrD →	DNA‑dependent ATPase I and helicase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,998,412	0	T	C	100.0%	141.9 / NA	40	W144R (TGG→CGG)	uvrD	DNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/19);  total (21/19)

GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAG  >  NC_000913/3998269‑3998555                                                                                                                                                |                                                                                                                                                gTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACAt                                                                                                                                          >  1:517831/1‑151 (MQ=255)     gctgcGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAAc                                                                                                                                      >  2:688973/1‑151 (MQ=255)        gcgtgcgCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACa                                                                                                                                     <  1:430150/149‑1 (MQ=255)          gtgcgCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCa                                                                                                                                 <  2:376698/151‑1 (MQ=255)              gcACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGaa                                                                                                                              >  2:432072/1‑150 (MQ=255)                  cATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATg                                                                                                                         <  2:207018/151‑1 (MQ=255)                    tatGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGa                                                                                                                        >  1:207680/1‑150 (MQ=255)                         aCGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGcc                                                                                                                   <  1:656362/150‑1 (MQ=255)                              aaTCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTc                                                                                                             <  2:150195/151‑1 (MQ=255)                                tCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCg                                                                                                           >  1:741099/1‑151 (MQ=255)                                 cTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGc                                                                                                          >  1:124504/1‑151 (MQ=255)                                    ccGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGcatc                                                                                                       <  2:479760/151‑1 (MQ=255)                                      gCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCATAGCCATGAACCTCGACTAGAAGCAGTGGTCGCCGCGGCACGCAATGCGGTACATCAACAGCCAGCAAGATGAAGGCCTGCCTCCGcatcat                                                                                                     >  1:744669/1‑151 (MQ=255)                                                   aTCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACg                                                                                        >  1:492615/1‑151 (MQ=255)                                                      cTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACgg                                                                                       >  1:256666/1‑149 (MQ=255)                                                                  gACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGc                                                                         <  2:409547/151‑1 (MQ=255)                                                                   aCCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGc                                                                         >  2:625885/1‑150 (MQ=255)                                                                    ccAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGACGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGc                                                                         >  1:444007/1‑149 (MQ=255)                                                                    ccAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAg                                                                       <  2:573356/151‑1 (MQ=255)                                                                          aCGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCt                                                                   <  1:599516/149‑1 (MQ=255)                                                                          aCGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTgg                                                                 <  2:114518/151‑1 (MQ=255)                                                                                gCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGaa                                                            <  2:560452/150‑1 (MQ=255)                                                                                   tAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGgtg                                                        >  2:236170/1‑151 (MQ=255)                                                                                    aaGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGgtgt                                                       <  2:637466/151‑1 (MQ=255)                                                                                    aaGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGgtgt                                                       <  2:741173/151‑1 (MQ=255)                                                                                           tGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGc                                                >  2:341027/1‑151 (MQ=255)                                                                                               cAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTAt                                            >  1:108419/1‑151 (MQ=255)                                                                                                aaaGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATc                                           >  1:265257/1‑151 (MQ=255)                                                                                                     cATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGaa                                      <  1:93114/151‑1 (MQ=255)                                                                                                      aTGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGa                                       >  1:144448/1‑149 (MQ=255)                                                                                                           ccTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCgtgt                                >  2:716074/1‑151 (MQ=255)                                                                                                            cTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGtgtg                               >  2:308837/1‑151 (MQ=255)                                                                                                                acgaGAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACcg                           <  2:317484/151‑1 (MQ=255)                                                                                                                  gagaAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACcgcg                         >  2:14163/1‑151 (MQ=255)                                                                                                                  gagaAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGAATCAGGAAGCGTGTGACcg                           >  1:328989/1‑149 (MQ=255)                                                                                                                     aaGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCgg                      >  1:774833/1‑151 (MQ=255)                                                                                                                         aGTGGCCGCCGCGGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCt                  <  2:619254/151‑1 (MQ=255)                                                                                                                                  cgcgGCAGGCAATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACtt         <  2:620731/151‑1 (MQ=255)                                                                                                                                       caggcaATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCg    <  1:624465/151‑1 (MQ=255)                                                                                                                                         ggcaATGCGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAg  <  1:487838/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTGGTACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAG  >  NC_000913/3998269‑3998555

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