Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,201,688	A→C	T64T (ACA→ACC)	ymfI →	e14 prophage; uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,201,688	0	A	C	100.0%	191.2 / NA	53	T64T (ACA→ACC)	ymfI	e14 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/26);  total (27/26)

AACAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGATT  >  NC_000913/1201543‑1201833                                                                                                                                                  |                                                                                                                                                  aaCAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTc                                                                                                                                              >  2:533461/1‑151 (MQ=255)  aCAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCa                                                                                                                                             <  2:718585/151‑1 (MQ=255)     gACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAAtt                                                                                                                                          <  2:901195/151‑1 (MQ=255)     gACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAAtt                                                                                                                                          >  2:789804/1‑151 (MQ=255)         aaaaaaTTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTa                                                                                                                                      <  1:378609/151‑1 (MQ=255)          aaaaaTTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTa                                                                                                                                      >  2:914346/1‑150 (MQ=255)             aaTTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAAc                                                                                                                                    >  2:364386/1‑149 (MQ=255)             aaTTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAg                                                                                                                                  <  1:178820/151‑1 (MQ=255)                  gaagaaTTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCg                                                                                                                             >  2:138473/1‑151 (MQ=255)                    agaaTTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGgt                                                                                                                           >  2:857922/1‑151 (MQ=255)                    agaaTTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGgt                                                                                                                           <  1:247826/151‑1 (MQ=255)                     gaaTTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGgtg                                                                                                                          <  2:379526/151‑1 (MQ=255)                      aaTTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGgtgt                                                                                                                         >  2:63066/1‑151 (MQ=255)                      aaTTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGATGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGgagt                                                                                                                         >  1:77703/1‑151 (MQ=255)                          aaTGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGa                                                                                                                      <  2:930253/150‑1 (MQ=255)                           aTGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACg                                                                                                                    <  1:78217/151‑1 (MQ=255)                                 aaCTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAat                                                                                                               >  1:412744/1‑150 (MQ=255)                                 aaCTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAat                                                                                                               <  1:770213/150‑1 (MQ=255)                                  aCTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAatat                                                                                                             <  2:429513/151‑1 (MQ=255)                                    tttCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAatatat                                                                                                           >  2:724410/1‑151 (MQ=255)                                      tCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATAttt                                                                                                         >  1:221467/1‑151 (MQ=255)                                       cTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTg                                                                                                        >  2:349557/1‑151 (MQ=255)                                        tGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGt                                                                                                       >  1:709809/1‑151 (MQ=255)                                            ccAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGt                                                                                                   >  1:232789/1‑151 (MQ=255)                                            ccAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGt                                                                                                   >  1:615188/1‑151 (MQ=255)                                                aaaaCAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAAt                                                                                               <  2:492892/151‑1 (MQ=255)                                                   aCAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTg                                                                                            <  1:164006/151‑1 (MQ=255)                                                            gagTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAgcg                                                                                   <  2:401996/151‑1 (MQ=255)                                                             agTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAgcgc                                                                                  <  2:857488/151‑1 (MQ=255)                                                             agTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAgcgc                                                                                  >  1:471327/1‑151 (MQ=255)                                                             agTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAgcgc                                                                                  <  1:419332/151‑1 (MQ=255)                                                                    cGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTa                                                                           <  2:402469/151‑1 (MQ=255)                                                                      aTAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACt                                                                         <  2:402993/151‑1 (MQ=255)                                                                      aTAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACt                                                                         <  2:777398/151‑1 (MQ=255)                                                                                     tGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCtt                                                          >  2:553518/1‑151 (MQ=255)                                                                                         tCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGt                                                       <  1:150163/150‑1 (MQ=255)                                                                                            agagGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGt                                                   >  2:787346/1‑151 (MQ=255)                                                                                               gggATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAAc                                                <  2:501887/151‑1 (MQ=255)                                                                                               gggATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAAc                                                <  2:534271/151‑1 (MQ=255)                                                                                                     tatCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCCTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTc                                          >  1:516402/1‑151 (MQ=255)                                                                                                       tCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGt                                        >  1:352/1‑151 (MQ=255)                                                                                                        ccTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGtt                                       >  2:752289/1‑151 (MQ=255)                                                                                                              aaGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTAc                                 >  1:478954/1‑151 (MQ=255)                                                                                                                              caacaaATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAg                  <  2:583275/150‑1 (MQ=255)                                                                                                                              caacaaATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGc                 >  2:720282/1‑151 (MQ=255)                                                                                                                               aacaaATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGcc                >  1:404754/1‑151 (MQ=255)                                                                                                                                     tCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCa           >  2:47157/1‑150 (MQ=255)                                                                                                                                      ccTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAg         <  2:759913/151‑1 (MQ=255)                                                                                                                                        ttGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTACGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCt       >  1:217347/1‑151 (MQ=255)                                                                                                                                            tCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGAt   >  1:223364/1‑151 (MQ=255)                                                                                                                                             cAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGAtt  <  2:209772/151‑1 (MQ=255)                                                                                                                                             cAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGAt   <  1:464834/150‑1 (MQ=255)                                                                                                                                             cAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGAt   <  1:533444/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  AACAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGATT  >  NC_000913/1201543‑1201833

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