Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,980,957	(G)6→7	coding (71/1386 nt)	yifK →	putative APC family amino acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,980,951	1	.	G	100.0%	209.6 / NA	59	L22R (CTG→CGG)	yifK	putative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (28/31);  total (28/31)

TTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAAT  >  NC_000913/3980815‑3981083                                                                                                                                          |                                                                                                                                      tttGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gg                                                                                                                                    >  2:22171/1‑139 (MQ=255)                  aCCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGt                                                                                                                    >  2:1168910/1‑139 (MQ=255)                     cATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCgg                                                                                                                 <  2:253214/139‑1 (MQ=255)                             tAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGGCACCATTGGCGTCGGCCt                                                                                                              >  1:33301/1‑135 (MQ=255)                             tAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGGCACCATTGGCGTCGGCCt                                                                                                              <  2:33301/135‑1 (MQ=255)                             tAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGc                                                                                                                >  2:1466904/1‑132 (MQ=255)                             tAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGc                                                                                                                <  1:1466904/132‑1 (MQ=255)                                  ggTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                    <  1:301729/139‑1 (MQ=255)                                  ggTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                    >  2:301729/1‑139 (MQ=255)                                   gTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATggggg                                                                                                   >  1:823518/1‑139 (MQ=255)                                    tAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                                                                  <  1:22171/139‑1 (MQ=255)                                          cGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATggggg                                                                                                   >  2:872072/1‑132 (MQ=255)                                          cGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATggggg                                                                                                   <  1:872072/132‑1 (MQ=255)                                           ggAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                    >  2:524901/1‑130 (MQ=255)                                           ggAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                    <  1:524901/130‑1 (MQ=255)                                            gAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa                                                                                          >  2:780288/1‑139 (MQ=255)                                            gAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa                                                                                          >  1:432265/1‑139 (MQ=255)                                              aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTg                                                                                                             <  1:396036/118‑1 (MQ=255)                                              aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTg                                                                                                             >  2:396036/1‑118 (MQ=255)                                                   tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaa                                                                                   <  2:269816/139‑1 (MQ=255)                                                     cAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                                                                 <  2:432265/139‑1 (MQ=255)                                                       aacacaacccGAGGAATTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGGAGGTCGACATATCGAACTCAACGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgg                                                                               <  2:494785/139‑1 (MQ=255)                                                           caacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATg                                                                                                       >  2:411147/1‑111 (MQ=255)                                                           caacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATg                                                                                                       <  1:411147/111‑1 (MQ=255)                                                              cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                        >  1:494785/1‑139 (MQ=255)                                                               acGAGGATCTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggcc                                                                       <  2:1288934/139‑1 (MQ=255)                                                                        tATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAtt                                                              >  2:999880/1‑139 (MQ=255)                                                                         aTGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttg                                                             <  2:422434/139‑1 (MQ=255)                                                                               gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                                                             <  2:667210/101‑1 (MQ=255)                                                                               gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                                                             >  1:667210/1‑101 (MQ=255)                                                                               gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc                                                       <  1:999880/139‑1 (MQ=255)                                                                               gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc                                                       >  1:1345746/1‑139 (MQ=255)                                                                                 tAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc                                                       >  2:603106/1‑137 (MQ=255)                                                                                 tAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc                                                       <  1:603106/137‑1 (MQ=255)                                                                                    cAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTatat                                                  <  2:229099/139‑1 (MQ=255)                                                                                            agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCg                                          >  2:1386834/1‑139 (MQ=255)                                                                                                  aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGt                                    <  1:1168910/139‑1 (MQ=255)                                                                                                      tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGt                                <  2:823518/139‑1 (MQ=255)                                                                                                      tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGt                                >  2:1498641/1‑139 (MQ=255)                                                                                                      tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggg                                                                                                                                   >  2:125397/1‑39 (MQ=255)                                                                                                      tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggg                                                                                                                                   <  1:125397/39‑1 (MQ=255)                                                                                                        ggCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCt                              <  1:1498641/139‑1 (MQ=255)                                                                                                         gCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGcc                                                                                                               >  1:1061696/1‑57 (MQ=255)                                                                                                         gCTGGAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGcc                                                                                                               <  2:1061696/57‑1 (MQ=255)                                                                                                             gAAGCTCGACATATCGATCTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt                                                           <  2:1466053/105‑1 (MQ=255)                                                                                                             gAAGCTCGACATATCGATCTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt                                                           >  1:1466053/1‑105 (MQ=255)                                                                                                             gAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATc                                                                    >  1:498795/1‑96 (MQ=255)                                                                                                             gAAGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATc                                                                    <  2:498795/96‑1 (MQ=255)                                                                                                              aaGCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTtca                        >  1:1121291/1‑139 (MQ=255)                                                                                                                gCTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTtcatc                      <  2:1345746/139‑1 (MQ=255)                                                                                                                 cTCGACATATCGAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTtcatca                     <  1:1386834/139‑1 (MQ=255)                                                                                                                            gAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt                                                            <  1:136569/89‑1 (MQ=255)                                                                                                                            gAACTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt                                                            >  2:136569/1‑89 (MQ=255)                                                                                                                              aCTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                    <  1:69761/47‑1 (MQ=255)                                                                                                                              aCTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                    >  2:69761/1‑47 (MQ=255)                                                                                                                              aCTCATCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATggg        <  2:389520/139‑1 (MQ=255)                                                                                                                                  aTCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                                                             <  1:632636/50‑1 (MQ=255)                                                                                                                                  aTCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                                                             >  2:632636/1‑50 (MQ=255)                                                                                                                                   tCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggc                                                                             >  1:879655/1‑65 (MQ=255)                                                                                                                                   tCGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggc                                                                             <  2:879655/65‑1 (MQ=255)                                                                                                                                    cGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAAt  >  2:1196359/1‑139 (MQ=255)                                                                                                                                    cGCCCTG‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAAt  >  2:1198147/1‑139 (MQ=255)                                                                                                                                          |                                                                                                                                      TTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAAT  >  NC_000913/3980815‑3981083

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