Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,762,713	T→C	T276T (ACA→ACG)	abpB ←	CP4‑57 prophage; putative helicase YfjK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,762,713	0	T	C	100.0%	166.0 / NA	46	T276T (ACA→ACG)	abpB	CP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/22);  total (24/22)

CGCAACAGTTGGTGACTTGCAATAAATGATCGTTGCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACTGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTGGAGACCTTAAGCAGTTTGCTTAACGCAATG  >  NC_000913/2762584‑2762839                                                                                                                                  |                                                                                                                               cccAACAGTTTGTGACTTGCAATGAATGTTCGTTGCGTCGATACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTa                                                                                                                             <  2:297863/131‑1 (MQ=255)         ttGGTGACTTGCAATAAATGATCGTTGCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTc                                                                                                                     >  1:614903/1‑133 (MQ=255)         ttGGTGACTTGCAATAAATGATCGTTGCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTc                                                                                                                     >  1:507381/1‑133 (MQ=255)           ggTGACTTGCAATAAATGATCGTTGCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCgg                                                                                                                   >  2:502313/1‑133 (MQ=255)                              tCGTTGCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGa                                                                                                          <  2:380456/123‑1 (MQ=255)                              tCGTTGCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGa                                                                                                          >  1:380456/1‑123 (MQ=255)                              tCGTTGCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCgg                                                                                                >  1:404742/1‑133 (MQ=255)                              tCGTTGCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCgg                                                                                                >  2:538027/1‑133 (MQ=255)                                   gCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATc                                                                                           >  1:526357/1‑133 (MQ=255)                                    cGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTATGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATcc                                                                                          >  2:391445/1‑133 (MQ=255)                                            caTGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAAAAACGAAAGTAGGCGGATATCAAAGttttt                                                                                  >  1:300171/1‑133 (MQ=255)                                             aTGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTc                                                                                 <  1:594264/133‑1 (MQ=255)                                                cgtgcgCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCgag                                                                              <  2:404742/133‑1 (MQ=255)                                                cgtgcgCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCgag                                                                              <  1:500941/133‑1 (MQ=255)                                                    cgcgATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACcg                                                                          >  2:13764/1‑133 (MQ=255)                                                    cgcgATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACcg                                                                          >  1:97874/1‑133 (MQ=255)                                                    cgcgATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACcg                                                                          >  2:408429/1‑133 (MQ=255)                                                    cgcgATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACcg                                                                          >  1:112334/1‑133 (MQ=255)                                                      cgATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACcgcg                                                                        <  2:97874/133‑1 (MQ=255)                                                                    gcTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTaacga                                                          <  1:87986/133‑1 (MQ=255)                                                                          aGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTc                                                    <  1:502313/133‑1 (MQ=255)                                                                            cTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTT‑ACGAACGGTCCCg                                                 <  1:53410/133‑1 (MQ=255)                                                                             ttGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCg                                                 <  1:408429/133‑1 (MQ=255)                                                                             tcggCTTTGTCGTCAGTCGCTTTGGTTCGGAAAGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCg                                                 <  2:300171/131‑1 (MQ=255)                                                                                       cgtcATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTaaaa                                       <  2:507381/133‑1 (MQ=255)                                                                                        gtcATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAAc                                      <  1:22010/133‑1 (MQ=255)                                                                                         tcATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTTGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACt                                     <  1:547099/133‑1 (MQ=255)                                                                                          cATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGGTAATGAACGGTCTCGTCGGGTAAAACTg                                    >  2:247866/1‑133 (MQ=255)                                                                                           aTTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGc                                   >  2:126752/1‑133 (MQ=255)                                                                                           aTTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGc                                   >  2:142592/1‑133 (MQ=255)                                                                                               gCTTTGATTCCGAACCTCATAATATCTAGCCCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTgg                               <  1:237483/133‑1 (MQ=255)                                                                                               gCTTTGATGCTGAACTTCTTAACACCTATCGCCACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTgg                               <  2:244350/133‑1 (MQ=255)                                                                                               gCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTgg                               <  1:601763/133‑1 (MQ=255)                                                                                               gCTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTgg                               <  1:538027/133‑1 (MQ=255)                                                                                                cTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGTATGCTGTTTACGAACGCTCCCGTCAGGTAAAACTGCCTGga                              >  2:236932/1‑133 (MQ=255)                                                                                                cTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCTCGTCAGGTaaa                                        >  2:513072/1‑123 (MQ=255)                                                                                                cTTTGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCTCGTCAGGTaaa                                        <  1:513072/123‑1 (MQ=255)                                                                                                   tGATGCCGAACGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATCCTGTTAACAAACTGTCGCGTCAGGTTAAACCGCCTGGAGAc                           >  2:237483/1‑133 (MQ=255)                                                                                                           aaCGTCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCTTCAGGTAAAACAGCATGGAGATCTTAAGCa                   >  1:272418/1‑133 (MQ=255)                                                                                                               tCTTAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTGGAGACCTTAAGCAGttt               >  1:625327/1‑133 (MQ=255)                                                                                                                 ttAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTGGAGACCTTAAGCAGTTTGc             <  1:407233/133‑1 (MQ=255)                                                                                                                 ttAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTGAAGGCCTTAAACAGTTTGc             >  1:298486/1‑133 (MQ=255)                                                                                                                  tAACATCTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTGGAGACCTTAAGCAGTTTGCt            <  2:526357/133‑1 (MQ=255)                                                                                                                        cTAGCGCGACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTGGAGACCTTAAGCAGTTTGCTTAACGc      >  1:617842/1‑133 (MQ=255)                                                                                                                            cgcgACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTGGAGACCTTAAGCAGTTTGCTTAACGCAATg  >  1:566021/1‑133 (MQ=255)                                                                                                                            cgcgACCGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTGGAGACCTTAAGCAGTTTGCTTAACGCAATg  <  1:142592/133‑1 (MQ=255)                                                                                                                                  |                                                                                                                               CGCAACAGTTGGTGACTTGCAATAAATGATCGTTGCGTCGACACATGCGTGCGCGATTTCCCCCAGCGCTTTGAGCTTGGCTTTGTCGTCATTCGCTTTGATGCCGAACGTCTTAACATCTAGCGCGACTGTATTGAAGTCGGTAGAAACGAAAGTATGCGGATATCCAAGTTTTTCGAGACCGCGAATGCTGTTAACGAACGGTCCCGTCAGGTAAAACTGCCTGGAGACCTTAAGCAGTTTGCTTAACGCAATG  >  NC_000913/2762584‑2762839

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