Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,868,399	C→T	C304Y (TGT→TAT)	yidR ←	DUF3748 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,868,399	0	C	T	100.0%	167.5 / NA	47	C304Y (TGT→TAT)	yidR	DUF3748 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/22);  total (25/22)

GACTGAATATCCGTTTTGTTATGGGTTAACTGGCGCGGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTACAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGCCACGCGGTGGCGCGGGCAGGGTTGTTTCCGTTCCGCTTAACGGCGCATCAC  >  NC_000913/3868273‑3868531                                                                                                                               |                                                                                                                                     gACTGAATATCCGTTTTGTTATGGGTTAACTGGCGCGGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGt                                                                                                                          <  2:574498/139‑1 (MQ=255)              ttttGTTATGGGTTAACTGGCGCGGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGtt                                                                                                             >  2:507808/1‑139 (MQ=255)                  gTTATGGGTTAACTGGCGCGGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACt                                                                                                         <  1:268788/139‑1 (MQ=255)                            aaCTGGCGCGGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGAt                                                                                               <  1:957482/139‑1 (MQ=255)                            aaCTGGCGCGGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGAt                                                                                               <  1:549140/139‑1 (MQ=255)                                 gcgcGGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGAACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGcc                                                                                          <  2:1025517/139‑1 (MQ=255)                                  cgcgGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGccc                                                                                         >  2:793968/1‑139 (MQ=255)                                   gcgGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGc                                                                                           >  1:1112807/1‑136 (MQ=255)                                   gcgGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGc                                                                                           <  2:1112807/136‑1 (MQ=255)                                       cTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAAcc                                                                                                     >  2:176009/1‑122 (MQ=255)                                       cTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAAcc                                                                                                     <  1:176009/122‑1 (MQ=255)                                                      cGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTaaa                                                                     <  2:1120606/139‑1 (MQ=255)                                                             aGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCt                                                              >  2:164710/1‑139 (MQ=255)                                                             aGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCt                                                              >  1:895789/1‑139 (MQ=255)                                                                   aGTTGCACAATGCCGTTTTCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCacga                                                        >  2:641933/1‑139 (MQ=255)                                                                   aGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCg                                                                <  2:94927/131‑1 (MQ=255)                                                                   aGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCg                                                                >  1:94927/1‑131 (MQ=255)                                                                       gCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGcc                                                    >  2:160000/1‑139 (MQ=255)                                                                        cacaATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAAcc                                                                                                     <  1:204074/89‑1 (MQ=255)                                                                        cacaATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAAcc                                                                                                     >  2:204074/1‑89 (MQ=255)                                                                         acaATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGgcgc                                                                                                                     >  1:1059922/1‑72 (MQ=255)                                                                         acaATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGgcgc                                                                                                                     <  2:1059922/72‑1 (MQ=255)                                                                         acaATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTaaa                                                                     <  1:827047/120‑1 (MQ=255)                                                                         acaATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTaaa                                                                     >  2:827047/1‑120 (MQ=255)                                                                               ccGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGcgcg                                                                                                                    >  1:48904/1‑67 (MQ=255)                                                                               ccGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGcgcg                                                                                                                    <  2:48904/67‑1 (MQ=255)                                                                               ccGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGCCACGCGGTg                                            >  2:443922/1‑139 (MQ=255)                                                                                    atcatcTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGcgcg                                                                                                                    >  2:466851/1‑62 (MQ=255)                                                                                    atcatcTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGcgcg                                                                                                                    <  1:466851/62‑1 (MQ=255)                                                                                             cATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGcc                                                    >  1:935057/1‑117 (MQ=255)                                                                                             cATTAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGcc                                                    <  2:935057/117‑1 (MQ=255)                                                                                               ttAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGt                                                                                                                          >  2:399645/1‑45 (MQ=255)                                                                                               ttAAAAACGCGATTTGTGTACCCTGCGGATTATAGCGCACCCAGt                                                                                                                          <  1:399645/45‑1 (MQ=255)                                                                                                        cgATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCa                                                           >  2:554584/1‑99 (MQ=255)                                                                                                        cgATTTGTGTACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCa                                                           <  1:554584/99‑1 (MQ=255)                                                                                                               tgtACCCTGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGCCACGCGGTGGCGCGGGCAGGGTTGTTTCCGTTCCGCTTAAc            <  1:641933/139‑1 (MQ=255)                                                                                                                      tGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGAc                                                                                                          <  1:467368/38‑1 (MQ=255)                                                                                                                      tGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGAc                                                                                                          >  2:467368/1‑38 (MQ=255)                                                                                                                      tGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCg                                                                                        >  2:587409/1‑56 (MQ=255)                                                                                                                      tGCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCg                                                                                        <  1:587409/56‑1 (MQ=255)                                                                                                                       gCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGCCACGCGGTGGCGCGGGCAGGGTTGTTTCCGTTCCGCTTAACGGCGCATc    >  2:354167/1‑139 (MQ=255)                                                                                                                       gCGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGCCACGCGGTGGCGCGGGCAGGGTTGTTTCCGTTCCGCTTAACGGCGCATc    >  2:582190/1‑139 (MQ=255)                                                                                                                        cGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGAt                                                                                               >  2:875637/1‑47 (MQ=255)                                                                                                                        cGGATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGAt                                                                                               <  1:875637/47‑1 (MQ=255)                                                                                                                         ggATTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGCCACGCGGTGGCGCGGGCAGGGTTGTTTCCGTTCCGCTTAACGGCGCATCAc  >  1:1031857/1‑139 (MQ=255)                                                                                                                           aTTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCg                                                                                        <  2:140767/51‑1 (MQ=255)                                                                                                                           aTTATAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCg                                                                                        >  1:140767/1‑51 (MQ=255)                                                                                                                               |                                                                                                                                     GACTGAATATCCGTTTTGTTATGGGTTAACTGGCGCGGCTCGCCGCCCTGTGGCGAGATAAGCCACAGTTGCACAATGCCGTTATCATCTCGCATTAAAAACGCGATTTGTGTACCCTGCGGATTACAGCGCACCCAGTGGCGCGGGACGTTGACTAACCCCGGATAAGCCCGATGGTGGGTAAAGGTTAAACGTCGCTGCACGACGCCACGCGGTGGCGCGGGCAGGGTTGTTTCCGTTCCGCTTAACGGCGCATCAC  >  NC_000913/3868273‑3868531

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