Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,197,576	T→C	D221G (GAC→GGC)	ymfE ←	e14 prophage; putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,197,576	0	T	C	100.0%	185.6 / NA	52	D221G (GAC→GGC)	ymfE	e14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/25);  total (27/25)

AATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAA  >  NC_000913/1197431‑1197720                                                                                                                                                  |                                                                                                                                                 aaTATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCggg                                                                                                                                             <  2:115924/151‑1 (MQ=255)  atatACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGa                                                                                                                                            <  1:439154/151‑1 (MQ=255)     taCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGAttt                                                                                                                                         >  2:236086/1‑151 (MQ=255)         cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGTGATTTTTcc                                                                                                                                     <  1:608097/151‑1 (MQ=255)         cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTcc                                                                                                                                     >  2:537623/1‑151 (MQ=255)                             aTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTgg                                                                                                                 <  2:235847/151‑1 (MQ=255)                               cGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGc                                                                                                               >  1:47615/1‑151 (MQ=255)                                gTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCt                                                                                                              >  1:247949/1‑151 (MQ=255)                                               gcATTCTCTGAGCGAATATTAACTCCGTGCATA‑TATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaa                                                                                              <  2:266413/151‑1 (MQ=255)                                                          gCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAcc                                                                                     >  1:385815/1‑150 (MQ=255)                                                          gCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACca                                                                                    <  1:557784/151‑1 (MQ=255)                                                              atatTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca                                                                                <  2:195109/151‑1 (MQ=255)                                                                atTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcact                                                                              >  1:61610/1‑151 (MQ=255)                                                                    aCTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTgaga                                                                          <  2:131648/151‑1 (MQ=255)                                                                       ccGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATg                                                                       <  1:215758/151‑1 (MQ=255)                                                                         gTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGAt                                                                     >  2:558605/1‑151 (MQ=255)                                                                           gCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATa                                                                    >  1:206868/1‑150 (MQ=255)                                                                            cATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaaa                                                                  <  2:546027/151‑1 (MQ=255)                                                                            cATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaa                                                                   >  2:185490/1‑150 (MQ=255)                                                                             atatTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaaa                                                                  >  1:73360/1‑150 (MQ=255)                                                                                  ataGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTAc                                                             >  2:267755/1‑150 (MQ=255)                                                                                   taGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACa                                                            >  2:516204/1‑150 (MQ=255)                                                                                   taGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACa                                                            >  1:241663/1‑150 (MQ=255)                                                                                    aGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaa                                                           >  1:341555/1‑150 (MQ=255)                                                                                      ttCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaaa                                                         <  2:387860/150‑1 (MQ=255)                                                                                            gCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCGCTGAGAATGATAAAATTACAAAAAAGAATa                                                  <  1:412557/151‑1 (MQ=255)                                                                                            gCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATa                                                  <  1:389957/151‑1 (MQ=255)                                                                                               aGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAg                                               <  2:112675/151‑1 (MQ=255)                                                                                                   ccATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAg                                           <  2:88137/151‑1 (MQ=255)                                                                                                    cATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGc                                          <  2:126924/151‑1 (MQ=255)                                                                                                    cATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGc                                          >  1:10741/1‑151 (MQ=255)                                                                                                    cATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGc                                          <  1:587527/151‑1 (MQ=255)                                                                                                       tCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCaa                                        <  1:311559/150‑1 (MQ=255)                                                                                                        cAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAg                                      <  1:547189/151‑1 (MQ=255)                                                                                                           ttAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAAc                                   >  2:269366/1‑151 (MQ=255)                                                                                                            tAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAAcc                                  >  1:570991/1‑151 (MQ=255)                                                                                                             aaTGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCt                                 >  2:176135/1‑151 (MQ=255)                                                                                                              aTGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCtt                                <  1:355571/151‑1 (MQ=255)                                                                                                                  tATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATGAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAAc                            >  2:600171/1‑151 (MQ=255)                                                                                                                  tATTTTTAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAAc                            >  2:411828/1‑151 (MQ=255)                                                                                                                       ttAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGa                        <  1:185495/150‑1 (MQ=255)                                                                                                                         aGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGaaa                      <  2:378596/150‑1 (MQ=255)                                                                                                                            gAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAat                   >  2:367606/1‑150 (MQ=255)                                                                                                                              aaaaTTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATc                >  1:352032/1‑151 (MQ=255)                                                                                                                                aaTTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATcc               >  1:542823/1‑150 (MQ=255)                                                                                                                                 aTTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCAt             >  1:487365/1‑151 (MQ=255)                                                                                                                                      tttattttatCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCAtttgt         <  1:417300/150‑1 (MQ=255)                                                                                                                                      tttattttatCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCAtttg          <  2:84813/149‑1 (MQ=255)                                                                                                                                      tttattttatCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCAttt           <  2:500475/148‑1 (MQ=255)                                                                                                                                          ttttatCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGc     >  1:16649/1‑150 (MQ=255)                                                                                                                                            ttatCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCaa   >  1:534362/1‑150 (MQ=255)                                                                                                                                             tatCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCaaa  >  1:445670/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                 AATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAA  >  NC_000913/1197431‑1197720

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