Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,806,796	A→C	T263P (ACC→CCC)	proW →	glycine betaine/proline ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,806,796	0	A	C	100.0%	187.9 / NA	50	T263P (ACC→CCC)	proW	glycine betaine/proline ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/30);  total (20/30)

GACGATCATCTTTGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTCGCGGTATCGGTC‑GTCTGGATATGGGGCTTGCCACCGTTGGCGGCG  >  NC_000913/2806648‑2806940                                                                                                                                                     |                                                                                                                                                   gACGATCATCTTTGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGccc                                                                                                                                                  <  1:636195/151‑1 (MQ=255)       catcTTTGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATg                                                                                                                                            <  2:899509/151‑1 (MQ=255)            ttGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATggcggg                                                                                                                                       >  1:38292/1‑151 (MQ=255)             tGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATggcggg                                                                                                                                       >  2:462891/1‑150 (MQ=255)              gCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATggcgggc                                                                                                                                      <  1:899890/150‑1 (MQ=255)                      gccgATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAg                                                                                                                             <  1:603675/151‑1 (MQ=255)                            tATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCt                                                                                                                        <  1:487542/150‑1 (MQ=255)                               ccGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATg                                                                                                                    >  1:683249/1‑151 (MQ=255)                               ccGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAAACAGACGc                                                                                                                         >  1:914856/1‑146 (MQ=255)                                  tCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTg                                                                                                                 >  1:409539/1‑151 (MQ=255)                                   cTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGAGTTAAACAGACGCTGCTGCTgg                                                                                                                >  1:466508/1‑151 (MQ=255)                                       ccATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGGTAACCAGAAGCTGATGCTGGccc                                                                                                             >  2:470821/1‑150 (MQ=255)                                           tCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTc                                                                                                         <  1:514298/150‑1 (MQ=255)                                            cTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCCTTCt                                                                                                        >  2:871623/1‑150 (MQ=255)                                                ggATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGGTAACCAGACGCTGATGCTGGCCCCTTCTAtgg                                                                                                    >  1:732541/1‑150 (MQ=255)                                                 gATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTAtggtg                                                                                                  <  2:380857/151‑1 (MQ=255)                                                  aTTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTAtggggg                                                                                                 >  1:635996/1‑151 (MQ=255)                                                          ggTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGcc                                                                                         <  1:703971/151‑1 (MQ=255)                                                             tCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTc                                                                                       <  2:329146/150‑1 (MQ=255)                                                                ggcggATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCgatg                                                                                   <  1:277951/151‑1 (MQ=255)                                                                ggcggATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCgat                                                                                    >  1:102261/1‑150 (MQ=255)                                                                   ggATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGGTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGAtt                                                                                >  2:880801/1‑151 (MQ=255)                                                                   ggATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCAGCCCATTATGGCGGGCGGTAACCAGAAACGGATGCTGGCCCCTTCTCTGGTGGGAACAGCCTCcatgat                                                                                 >  1:269348/1‑150 (MQ=255)                                                                       cTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGAAGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGcc                                                                             >  1:895844/1‑150 (MQ=255)                                                                         gATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTg                                                                               <  2:757905/146‑1 (MQ=255)                                                                          aTTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTcg                                                                         <  1:279315/151‑1 (MQ=255)                                                                          aTTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTcg                                                                         <  1:91163/151‑1 (MQ=255)                                                                          aTTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTc                                                                          <  1:725480/150‑1 (MQ=255)                                                                             gAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCTCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTcggc                                                                       <  2:667611/150‑1 (MQ=255)                                                                              aaGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTcggc‑gg                                                                    <  2:242631/151‑1 (MQ=255)                                                                                 ccTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGtt                                                                 <  2:31023/151‑1 (MQ=255)                                                                                   tCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGGCGGC‑gggg                                                                  >  2:836760/1‑147 (MQ=255)                                                                                    cgcgCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTgg                                                               <  2:522685/150‑1 (MQ=255)                                                                                             tCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATgg                                                      <  1:606117/150‑1 (MQ=255)                                                                                             tCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGGTGGGGCCGATgg                                                      >  1:709249/1‑150 (MQ=255)                                                                                              cGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGt                                                     <  2:58828/150‑1 (MQ=255)                                                                                              cGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTa                                                    <  1:181462/151‑1 (MQ=255)                                                                                                gTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACt                                                  <  1:138210/151‑1 (MQ=255)                                                                                                  gCCAGCCCGCGCCAGATGCTGTTCAAAGTTTAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACtt                                                 <  1:777840/150‑1 (MQ=255)                                                                                                   ccAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTcg                                               <  2:287382/151‑1 (MQ=255)                                                                                                      gCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTCGCgg                                            <  2:635981/151‑1 (MQ=255)                                                                                                            gcCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTCGCGGTATCgg                                      <  2:527156/151‑1 (MQ=255)                                                                                                                gATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTCGCGGTATCGgtc‑g                                  <  2:754629/150‑1 (MQ=255)                                                                                                                   gCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTCGCGGTATCGGTC‑GTCt                               <  1:882297/150‑1 (MQ=255)                                                                                                                   gCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTCGCGGTATCGGTC‑GTCTg                              <  2:363936/151‑1 (MQ=255)                                                                                                                        tCAAAGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGCTGATTGCCGTCGGC‑GGGGTGGGGCAGAATGTACTTCGCGGGATCGGTT‑GGCTggg                            >  1:228986/1‑147 (MQ=255)                                                                                                                          aaaGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTCGCGGTATCGGTC‑GTCTGGATATgg                       <  1:545859/151‑1 (MQ=255)                                                                                                                           aaGTTCAGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTCGCGGTATCGGTC‑GTCTGGATATggg                      <  1:701621/151‑1 (MQ=255)                                                                                                                                 aGTTACCGCTGGCGATGCCGCCCATTATGGCGGGCGTTTACCAGACGCTGATGCTGGCCCCTTCTATGGGGGGCGTCGCCTCGATGATTGCCGTCGGCGGGGTTGGGT‑AGATGGGACTTTGCGGGAACGGTCTGCCGGGGTATGGGGCt                   >  2:831959/1‑149 (MQ=255)                                                                                                                                           ggCGATGCCGCCCATTATGGCGGGCGGTAACCACACGCTGATGCTGGCCCTTTCTATGGGGGTCATCCCCCCGATGATTTCCGGCGGC‑GGGGTGGGGCCGATGGGACTTTGCGGGATCCGTC‑GCCTGGGTATGGGGCGTGCcccccgtggg      >  2:771724/1‑145 (MQ=255)                                                                                                                                                  ccGCCCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGGTGGGTCAGATGGGACTTCGCGGGATCCGGC‑GTCTGGGTATGGGGGTTGCCCCCGCTGgcggcg  >  1:114371/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GACGATCATCTTTGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGC‑GGGTTGGGTCAGATGGTACTTCGCGGTATCGGTC‑GTCTGGATATGGGGCTTGCCACCGTTGGCGGCG  >  NC_000913/2806648‑2806940

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