Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,398,576	T→C	D16G (GAC→GGC)	uspE ←	universal stress protein with a role cellular motility

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,398,576	0	T	C	100.0%	212.0 / NA	58	D16G (GAC→GGC)	uspE	universal stress protein with a role cellular motility
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/34);  total (24/34)

TAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGTGATGAACCTTCTGC  >  NC_000913/1398445‑1398711                                                                                                                                    |                                                                                                                                        tAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGTCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTa                                                                                                                                 <  2:1148680/140‑1 (MQ=255) tAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTa                                                                                                                                 >  2:1269914/1‑140 (MQ=255)  aGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAg                                                                                                                                >  1:837817/1‑140 (MQ=255)     ggTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCATCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGAt                                                                                                                             <  2:55468/140‑1 (MQ=255)       tACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCg                                                                                                                           >  2:563423/1‑140 (MQ=255)         cGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGAt                                                                                                                         >  1:367766/1‑140 (MQ=255)          gTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATa                                                                                                                        >  2:480786/1‑140 (MQ=255)              gTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATaacaa                                                                                                                    >  2:837629/1‑140 (MQ=255)                   gggAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGc                                                                                                               >  1:140729/1‑140 (MQ=255)                    ggAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCa                                                                                                              <  2:837817/140‑1 (MQ=255)                    ggAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCa                                                                                                              <  2:367766/140‑1 (MQ=255)                     gagaGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCAt                                                                                                             <  2:532007/140‑1 (MQ=255)                      agagCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATg                                                                                                            <  2:1049403/140‑1 (MQ=255)                      agagCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATg                                                                                                            <  1:94806/140‑1 (MQ=255)                      agagCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATg                                                                                                            <  1:944005/140‑1 (MQ=255)                          cAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCt                                                                                                        <  1:169426/140‑1 (MQ=255)                              gtggtCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATa                                                                                                    <  2:1221402/140‑1 (MQ=255)                              gtggtCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATa                                                                                                    >  1:1221402/1‑140 (MQ=255)                              gtggtCATTTCGTATGAGAAGTCATAGAGCGGCAAAAAGGCTTTAATTTTGCCGCCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATa                                                                                                    <  2:151779/140‑1 (MQ=255)                                  tCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGc                                                                                                               <  2:758406/125‑1 (MQ=255)                                  tCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGc                                                                                                               >  1:758406/1‑125 (MQ=255)                                  tCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATa                                                                                                <  2:300233/140‑1 (MQ=255)                                  tCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCAATAACAACGAGCATGTTCTGATACATa                                                                                                <  2:1172849/140‑1 (MQ=255)                                       tCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCAt                                                                                           >  2:1061322/1‑140 (MQ=255)                                         gTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATAc                                                                                         <  1:563423/140‑1 (MQ=255)                                           aTGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATa                                                                                                >  2:825491/1‑131 (MQ=255)                                           aTGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATa                                                                                                <  1:825491/131‑1 (MQ=255)                                              agaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATGACAACGAGCATGTTCTGATACATAGCCATACAGGGt                                                                                    <  1:1269914/140‑1 (MQ=255)                                                 aGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtctc                                                                                 <  1:364204/140‑1 (MQ=255)                                                          tCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTacaaca                                                                        <  2:1076719/140‑1 (MQ=255)                                                             gCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACt                                                                     >  1:1261550/1‑140 (MQ=255)                                                             gCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACt                                                                     >  1:468337/1‑140 (MQ=255)                                                             gCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTAGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACt                                                                     >  2:1126656/1‑140 (MQ=255)                                                                    ggCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACg                                                              <  1:389694/140‑1 (MQ=255)                                                                     gCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGc                                                             <  1:837629/140‑1 (MQ=255)                                                                               tGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAAt                                                   >  2:1090526/1‑140 (MQ=255)                                                                                gCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAAtt                                                  >  2:1193632/1‑140 (MQ=255)                                                                                  caccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTaa                                                <  2:213840/140‑1 (MQ=255)                                                                                  caccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTaa                                                <  1:1090526/140‑1 (MQ=255)                                                                                       aTCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGAt                                           <  1:1193632/140‑1 (MQ=255)                                                                                       aTCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGAt                                           <  1:480786/140‑1 (MQ=255)                                                                                       aTCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGAt                                           <  2:468337/140‑1 (MQ=255)                                                                                        tCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGAtt                                          <  2:965678/140‑1 (MQ=255)                                                                                        tCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGAtt                                          <  2:140729/140‑1 (MQ=255)                                                                                                         aCAGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGtgaa                         <  1:802429/140‑1 (MQ=255)                                                                                                           aGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGtgaa                         >  2:216892/1‑138 (MQ=255)                                                                                                           aGCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGtgaa                         <  1:216892/138‑1 (MQ=255)                                                                                                            gCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGtgaatga                      >  1:785973/1‑140 (MQ=255)                                                                                                            gCTCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGtgaatga                      >  2:904132/1‑140 (MQ=255)                                                                                                              tCGCCGCAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGaaa                    <  1:904132/140‑1 (MQ=255)                                                                                                                gccgcAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGaaa                    >  1:1255305/1‑138 (MQ=255)                                                                                                                gccgcAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGaaa                    <  2:1255305/138‑1 (MQ=255)                                                                                                                  cgcAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGaaa                    >  2:284244/1‑136 (MQ=255)                                                                                                                  cgcAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGaaa                    <  1:284244/136‑1 (MQ=255)                                                                                                                  cgcAATGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGt                >  2:770832/1‑140 (MQ=255)                                                                                                                      aTGCTGGTTGGTCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGtgatg            <  1:929323/140‑1 (MQ=255)                                                                                                                            gttggtCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGTGATGAACCtt      >  2:1227866/1‑140 (MQ=255)                                                                                                                               ggtCGCCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGTGATGAACCTTCTg   >  1:313715/1‑140 (MQ=255)                                                                                                                                gtcgccCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGTGATGAACCTTCTGc  <  1:1061322/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        TAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGTGATGAACCTTCTGC  >  NC_000913/1398445‑1398711

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