Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,223,567	A→G	S100P (TCG→CCG)	dusC ←	tRNA dihydrouridine(16) synthase DusC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,223,567	0	A	G	100.0%	131.2 / NA	38	S100P (TCG→CCG)	dusC	tRNA dihydrouridine(16) synthase DusC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/24);  total (14/24)

GCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGACGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCG  >  NZ_CP009273/2223436‑2223700                                                                                                                                    |                                                                                                                                      gcccagACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAg                                                                                                                                >  2:504466/1‑139 (MQ=255)    cagACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCg                                                                                                                             <  1:171065/139‑1 (MQ=255)      gACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGc                                                                                                                            >  1:207356/1‑138 (MQ=255)      gACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGc                                                                                                                            >  1:422336/1‑138 (MQ=255)         gCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAAtt                                                                                                                        <  1:137684/139‑1 (MQ=255)                 aCGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATcc                                                                                                                 <  2:419840/138‑1 (MQ=255)                 aCGCGGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGGCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCa                                                                                                                <  2:578004/139‑1 (MQ=255)                  cGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCAc                                                                                                               <  2:427670/139‑1 (MQ=255)                       aCGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCTTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGccc                                                                                                           >  2:93285/1‑138 (MQ=255)                        cGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCa                                                                                                         <  1:546714/139‑1 (MQ=255)                         gggCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCCCCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAg                                                                                                        <  2:382859/139‑1 (MQ=255)                                       gTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCAc                                                                                          >  1:510099/1‑139 (MQ=255)                                                cccGCATCGCTTTTGCCCCCTGGTAGATGAGTCCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACggg                                                                                  <  2:285478/136‑1 (MQ=255)                                                 aCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGggcg                                                                                <  1:531570/139‑1 (MQ=255)                                                  cGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCCCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGCTCCACGCCCCAGGAACCTAACTCCACCGCACGggcgg                                                                               >  2:9829/1‑139 (MQ=255)                                                      tCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCGCGGGCGGCGtt                                                                           <  2:3005/139‑1 (MQ=255)                                                               cACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTgccag                                                                  <  2:422336/139‑1 (MQ=255)                                                                    tGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCAtt                                                             >  2:386495/1‑139 (MQ=255)                                                                           tGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGaa                                                      <  2:510099/139‑1 (MQ=255)                                                                           tGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGaa                                                      <  1:13225/139‑1 (MQ=255)                                                                           tGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGaa                                                      >  1:324405/1‑139 (MQ=255)                                                                           tGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGaa                                                      <  1:386495/139‑1 (MQ=255)                                                                             aGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAAc                                                     <  1:11232/138‑1 (MQ=255)                                                                             aGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACt                                                    <  2:207356/139‑1 (MQ=255)                                                                                tCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGAc                                                 >  2:286310/1‑139 (MQ=255)                                                                                 cAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGAcc                                                >  1:523788/1‑139 (MQ=255)                                                                                 cAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGAc                                                 <  1:286310/138‑1 (MQ=255)                                                                                        tttGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTaacaac                                         >  1:344178/1‑139 (MQ=255)                                                                                         ttGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACt                                                    >  1:243848/1‑127 (MQ=255)                                                                                         ttGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACt                                                    <  2:243848/127‑1 (MQ=255)                                                                                          tGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACt                                        <  1:93285/138‑1 (MQ=255)                                                                                              tAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGaa                                                      <  1:469026/120‑1 (MQ=255)                                                                                              tAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGaa                                                      >  2:469026/1‑120 (MQ=255)                                                                                                aCGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCAcgcg                                 <  1:18124/139‑1 (MQ=255)                                                                                                 cGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGcccc                                                                                                          <  1:63036/65‑1 (MQ=255)                                                                                                 cGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGcccc                                                                                                          >  2:63036/1‑65 (MQ=255)                                                                                                       cccGCCGCCGCTACCGTTAACCGTTTTCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAAATGCACGCGCACCAGc                          <  1:504466/139‑1 (MQ=255)                                                                                                                                ttCGGCGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCg  <  2:523788/138‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      GCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGACGGGCAGCCGCAATTGAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCG  >  NZ_CP009273/2223436‑2223700

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