Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,078,709	G→A	A45A (GCG→GCA)	yoeF →	CP4‑44 prophage; putative uncharacterized protein YoeF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,078,709	0	G	A	100.0%	161.5 / NA	48	A45A (GCG→GCA)	yoeF	CP4‑44 prophage; putative uncharacterized protein YoeF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/20);  total (28/20)

TCCATGGCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCGCTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGG  >  NC_000913/2078572‑2078826                                                                                                                                          |                                                                                                                      tCCATGGCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACt                                                                                                                     >  1:95449/1‑140 (MQ=255)   cATGGCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGGTGGTCTGGCAACAAGAAGGAATCTGAGCACTGa                                                                                                                   >  2:325244/1‑140 (MQ=255)   cATGGCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGa                                                                                                                   >  2:629902/1‑140 (MQ=255)   cATGGCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGa                                                                                                                   >  1:429077/1‑140 (MQ=255)   cATGGCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGa                                                                                                                   >  2:401956/1‑140 (MQ=255)   cATGGCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGAGCCTGCTGTTCTGGCATCAAGCAGGAATCTCAGCACTGa                                                                                                                   >  2:234934/1‑140 (MQ=255)      ggCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGc                                                                                                                <  1:401956/140‑1 (MQ=255)      ggCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGc                                                                                                                >  2:510999/1‑140 (MQ=255)       gCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCa                                                                                                               <  2:466381/140‑1 (MQ=255)             gTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAAtt                                                                                                         <  1:629902/140‑1 (MQ=255)             gTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAAtt                                                                                                         >  1:622327/1‑140 (MQ=255)                cGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTaga                                                                                                      >  1:542007/1‑140 (MQ=255)                      cacGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCGGGAATCTCAGCAATGAAGCAGGAATTAGAGGTaaa                                                                                                >  2:267433/1‑140 (MQ=255)                      cacGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTGCTGGCAACGAGCAGGACTCTCAGCAATGAATCAGGAATTAGAGGTaaa                                                                                                >  2:257329/1‑140 (MQ=255)                            tGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATTTCAGCACTGAAGAAGGAGTTAGAGGTAAAAACGcc                                                                                          >  2:275481/1‑140 (MQ=255)                                gggCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTAc                                                                                      <  1:518433/140‑1 (MQ=255)                                    cATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCggg                                                                                  >  1:564006/1‑140 (MQ=255)                                     aTTCCGGCAGTATGCTCAGAATCTGTCGCGGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGc                                                                                 <  1:234934/140‑1 (MQ=255)                                     aTTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGGAAAAACGCCCTACCGGGc                                                                                 >  1:439139/1‑140 (MQ=255)                                        gagacAGTATGCTCATAATATGCCGCAGATACGCTTCGGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTGAAAACGCCCTACCGGGCCAt                                                                              <  1:296821/136‑1 (MQ=255)                                               tATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCAc                                                                       <  1:564122/140‑1 (MQ=255)                                                 tGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGc                                                                                 >  2:291472/1‑128 (MQ=255)                                                 tGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGc                                                                                 <  1:291472/128‑1 (MQ=255)                                                 tGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGc                                                                                 >  1:349297/1‑128 (MQ=255)                                                 tGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGc                                                                                 <  2:349297/128‑1 (MQ=255)                                                 tGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCAccc                                                                     >  1:47110/1‑140 (MQ=255)                                                         atatGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTg                                                             <  2:429077/140‑1 (MQ=255)                                                         atatGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTg                                                             >  1:370188/1‑140 (MQ=255)                                                         atatGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGAGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTg                                                             >  2:542195/1‑140 (MQ=255)                                                          tatGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGGGAAAACGGCCTACCGGACCATAAATCACCCggctgtctc                                                            >  2:296821/1‑133 (MQ=255)                                                             ggCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGTCCTTGAATCACCCGGTCATTGGagt                                                         >  1:305160/1‑140 (MQ=255)                                                             ggCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGagt                                                         >  1:464531/1‑140 (MQ=255)                                                                  aGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTaa                                                    >  1:514184/1‑140 (MQ=255)                                                                    aTACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTGAAAACGCCTTACGGGGCCATGAATCCCCCGGTCATTGGAGTAGTAAcc                                                  >  2:223947/1‑140 (MQ=255)                                                                        gCTTCCGGATCGATACCGTTCACACCGCACGTGCCTGGTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAg                                              <  1:257329/140‑1 (MQ=255)                                                                                       ccGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCa                               <  1:612340/140‑1 (MQ=255)                                                                                        cGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTAGAGTAGTAACCAACGCAGATTTAGCCAGCAt                              >  2:215444/1‑140 (MQ=255)                                                                                         gTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATg                             <  1:4715/140‑1 (MQ=255)                                                                                         gTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATg                             <  2:95449/140‑1 (MQ=255)                                                                                          ttCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATgg                            <  2:404/140‑1 (MQ=255)                                                                                               aCGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACaa                       >  2:17615/1‑140 (MQ=255)                                                                                                      gTGCCTGCTTTTCTGTCAACAAGCAGGAATCTCAGCAGTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCtt                <  1:267433/140‑1 (MQ=255)                                                                                                      gTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCtt                <  1:510999/140‑1 (MQ=255)                                                                                                      gTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCtt                <  1:428733/140‑1 (MQ=255)                                                                                                         ccTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTgg             >  2:371648/1‑140 (MQ=255)                                                                                                          cTGCTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGt            <  2:370188/140‑1 (MQ=255)                                                                                                             cTTTTCTGGCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGaa         >  2:52403/1‑140 (MQ=255)                                                                                                                    ggCAACAAGCAGGAATCTCAGCACTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTgg  <  1:542195/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                      TCCATGGCAGGAGTTCGTCAACACGGTTGGAGGGCCATTCCGGCAGTATGCTCAGAATATGGCGCAGATACGCTTCCGGATCGATACCGTTCAGACGGCAGGTGCCTGCTTTTCTGGCAACAAGCAGGAATCTCAGCGCTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGG  >  NC_000913/2078572‑2078826

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