Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,981,539	T→G	L218R (CTG→CGG)	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,981,539	0	T	G	100.0%	180.1 / NA	50	L218R (CTG→CGG)	yifK	putative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/25);  total (25/25)

GGCGTGATTTTCTTTGGCTTTGGCAATGGCGGGCAGTCGATTGGTTTTAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCTGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGTGGAATGAAATAGGCAGCA  >  NC_000913/3981394‑3981685                                                                                                                                                  |                                                                                                                                                   ggCGTGATTTTCTTTGGCTTTGGCAATGGCGGGCAGTCGATTGGTTTTAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTg                                                                                                                                               <  2:466818/151‑1 (MQ=255)    gTGATTTTCTTTGGCTTTGGCAATGGCGGGCAGTCGATTGGTTTTAGCAATCTCACAGAGCATGGCGGGTTCTTTGCGGGGGGCTGGAAAGGGGTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTgg                                                                                                                                              >  1:314241/1‑149 (MQ=255)                ggctttggcAATGGCGGGCAGTCGATTGGTTTTAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTg                                                                                                                                >  2:805223/1‑151 (MQ=255)                                    gTCGATTGGTTTTAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGt                                                                                                             <  2:260344/150‑1 (MQ=255)                                         ttGGTTTTAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCt                                                                                                       >  2:776803/1‑151 (MQ=255)                                           ggTTTTAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTgc                                                                                                     <  1:173116/151‑1 (MQ=255)                                             ttttAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTgcgc                                                                                                   >  2:401260/1‑151 (MQ=255)                                               ttAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCa                                                                                                  <  1:482204/150‑1 (MQ=255)                                               ttAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGGGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTgcgc                                                                                                   >  2:577115/1‑149 (MQ=255)                                                tAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGTTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGGGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGGGACGCTGGGCa                                                                                                  >  2:131760/1‑149 (MQ=255)                                                   cAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGcc                                                                                             >  2:648454/1‑151 (MQ=255)                                                    aaTCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCCGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCg                                                                                            >  1:769015/1‑151 (MQ=255)                                                          acaGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCGTTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAgg                                                                                        <  2:728827/149‑1 (MQ=255)                                                          acaGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCa                                                                                      <  2:365978/151‑1 (MQ=255)                                                             gagCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAgg                                                                                   <  1:765948/151‑1 (MQ=255)                                                               gCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTg                                                                                 <  2:825483/151‑1 (MQ=255)                                                                 aTGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCt                                                                               >  2:516893/1‑151 (MQ=255)                                                                   ggCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCtgt                                                                             <  2:352924/151‑1 (MQ=255)                                                                   ggCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCtgt                                                                             >  2:719520/1‑151 (MQ=255)                                                                    gCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTGCCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCtgtg                                                                            <  2:370404/151‑1 (MQ=255)                                                                       gTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCg                                                                         >  2:125742/1‑151 (MQ=255)                                                                           cTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGaac                                                                     >  1:352927/1‑149 (MQ=255)                                                                             ttGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATcc                                                                    <  2:917653/150‑1 (MQ=255)                                                                              tGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCt                                                                   <  1:501123/150‑1 (MQ=255)                                                                                         ggAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACg                                                        <  1:236697/150‑1 (MQ=255)                                                                                           aaaGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACg                                                        >  1:195698/1‑148 (MQ=255)                                                                                            aaGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAgg                                                    >  2:634311/1‑151 (MQ=255)                                                                                               ggTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGAACCTGATTTTCTACGTaag                                                    >  1:931723/1‑146 (MQ=255)                                                                                                gTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGcgcg                                                <  1:683981/151‑1 (MQ=255)                                                                                                      gACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGAttt                                            <  2:142184/149‑1 (MQ=255)                                                                                                             ctGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACTCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATc                                    >  2:665567/1‑150 (MQ=255)                                                                                                               gtgtATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTAt                                     <  2:303571/147‑1 (MQ=255)                                                                                                                 gtATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCAc                               >  1:17225/1‑151 (MQ=255)                                                                                                                 gtATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCAc                               <  1:520064/151‑1 (MQ=255)                                                                                                                  tATTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCAcc                              >  2:869394/1‑151 (MQ=255)                                                                                                                   aTTGTGGTGGCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCAcc                              >  1:695052/1‑150 (MQ=255)                                                                                                                        ggtggCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTc                        <  1:916310/151‑1 (MQ=255)                                                                                                                           ggCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTccc                      <  1:17222/150‑1 (MQ=255)                                                                                                                           ggCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTccc                      <  1:530459/150‑1 (MQ=255)                                                                                                                           ggCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCg                     >  2:787794/1‑151 (MQ=255)                                                                                                                           ggCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCg                     <  1:216613/151‑1 (MQ=255)                                                                                                                            gCGTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGt                    <  2:5202/151‑1 (MQ=255)                                                                                                                              gTCCTACCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGTg                   >  1:756271/1‑150 (MQ=255)                                                                                                                                   aCCAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGTGGAATGa             <  1:805203/151‑1 (MQ=255)                                                                                                                                     cAGGGCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGTGGAATGaa            >  2:384517/1‑150 (MQ=255)                                                                                                                                       gggCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCCCCATCTTCCCGTGGAATGaaa           >  1:365930/1‑149 (MQ=255)                                                                                                                                         gCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGTGGAATGAAATAgg       >  2:313628/1‑151 (MQ=255)                                                                                                                                         gCGTGGAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGTGGAATGAAATAgg       <  1:488262/151‑1 (MQ=255)                                                                                                                                              gAGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGTGGAATGAAATAGgcagca  <  2:765210/151‑1 (MQ=255)                                                                                                                                               aGCGGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGTGGAATGAAATAGgcagca  >  1:161427/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GGCGTGATTTTCTTTGGCTTTGGCAATGGCGGGCAGTCGATTGGTTTTAGCAATCTCACAGAGCATGGCGGTTTCTTTGCGGGTGGCTGGAAAGGGTTCCTGACCGCTCTGTGTATTGTGGTGGCGTCCTACCAGGGCGTGGAGCTGATTGGCATTACTGCCGGTGAAGCGAAGAATCCGCAGGTGACGCTGCGCAGTGCCGTAGGCAAGGTGCTGTGGCGGATCCTGATTTTCTACGTAGGCGCGATTTTCGTTATCGTCACCATCTTCCCGTGGAATGAAATAGGCAGCA  >  NC_000913/3981394‑3981685

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