Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,070,778	T→C	intergenic (+304/‑8)	yoeA → / → yeeP	CP4‑44 prophage; TonB‑dependent receptor plug domain‑containing protein YoeA/CP4‑44 prophage; putative GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,070,778	0	T	C	100.0%	135.1 / NA	40	intergenic (+304/‑8)	yoeA/yeeP	CP4‑44 prophage; TonB‑dependent receptor plug domain‑containing protein YoeA/CP4‑44 prophage; putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/26);  total (14/26)

TCTCGCGTAAGCGGGTACCGTGACATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGT  >  NC_000913/2070645‑2070913                                                                                                                                      |                                                                                                                                        tctcGCGTAAGCGGGTACCGTGACATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGtt                                                                                                                                   <  1:315779/140‑1 (MQ=255)  ctcGCGTAAGCGGGTACCGTGACATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGttt                                                                                                                                  <  1:660469/140‑1 (MQ=255)        tAAGCGGGTACCGTGACATTCTGCCTGAACTCGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTgtggtg                                                                                                                            >  1:791060/1‑140 (MQ=255)           gCGGGTACCGTGACATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACg                                                                                                                         <  1:752190/140‑1 (MQ=255)                       aCATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACaaa                                                                                                             >  1:481009/1‑140 (MQ=255)                             tGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAg                                                                                                       <  2:481009/140‑1 (MQ=255)                             tGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAg                                                                                                       <  1:31117/140‑1 (MQ=255)                               ccTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGcc                                                                                                     >  2:203124/1‑140 (MQ=255)                               ccTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACCCAGGCCGACAAAACGGAGcc                                                                                                     >  1:393741/1‑140 (MQ=255)                                 tGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCt                                                                                                   <  2:160103/140‑1 (MQ=255)                                    aCTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGcc                                                                                                <  1:203124/140‑1 (MQ=255)                                    aCTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGcc                                                                                                <  1:411487/140‑1 (MQ=255)                                    aCTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGcc                                                                                                <  1:656167/140‑1 (MQ=255)                                        gACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGa                                                                                            <  2:791060/140‑1 (MQ=255)                                         aCCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                                                                           >  1:747663/1‑140 (MQ=255)                                         aCCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                                                                           <  1:459726/140‑1 (MQ=255)                                         aCCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                                                                           >  2:364004/1‑140 (MQ=255)                                             ggTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTTGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCTTGCCATGAATggg                                                                                       >  1:92046/1‑140 (MQ=255)                                             ggTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATggg                                                                                       <  2:227255/140‑1 (MQ=255)                                               tACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGat                                                                                     <  1:600075/140‑1 (MQ=255)                                                aCTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGata                                                                                    >  1:865923/1‑140 (MQ=255)                                                    tgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTTGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCTTGCCATGAATGGGATATGGc                                                                                <  2:92046/140‑1 (MQ=255)                                                                gCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGcc                                                                    <  1:145878/140‑1 (MQ=255)                                                                    aTGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTtct                                                                >  1:465027/1‑140 (MQ=255)                                                                          gTGCCGTGTCTGTGGATGAGTATTTCTGCCGACACACCCGGCAGAGCGGACCTTAACCTGCTCTGTTTGGGCTGCCGCGGGCCGACAAACCGGACCGCAGCCATGAATGAGATATTGCTTTTCTTCAGCCTTCTCCTgct                                                          <  2:757780/140‑2 (MQ=255)                                                                          gTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgca                                                          <  1:729281/140‑1 (MQ=255)                                                                          gTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgca                                                          <  2:662813/140‑1 (MQ=255)                                                                           tGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcac                                                         <  2:740865/140‑1 (MQ=255)                                                                           tGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcac                                                         >  1:298387/1‑140 (MQ=255)                                                                                            aGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTcgc                                        >  2:74851/1‑140 (MQ=255)                                                                                             gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTcgcg                                       >  2:203953/1‑140 (MQ=255)                                                                                                     ggCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACg                               >  1:757004/1‑140 (MQ=255)                                                                                                      gCGACACATCCTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACgg                              <  2:359483/140‑1 (MQ=255)                                                                                                               ccTGCAGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTAt                     <  1:74851/140‑1 (MQ=255)                                                                                                                    aGTGCGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGt                <  2:471207/140‑1 (MQ=255)                                                                                                                        cGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTAt                     >  1:184456/1‑131 (MQ=255)                                                                                                                        cGGACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTAt                     <  2:184456/131‑1 (MQ=255)                                                                                                                         ggACATCAGCAGGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGtt           <  2:747663/140‑1 (MQ=255)                                                                                                                               cagcagGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGcc     <  1:364004/140‑1 (MQ=255)                                                                                                                               cagcagGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGcc     <  2:298387/140‑1 (MQ=255)                                                                                                                                  cagGCGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGt  <  2:757004/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        TCTCGCGTAAGCGGGTACCGTGACATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGT  >  NC_000913/2070645‑2070913

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: ‑