Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,418,386	G→A	intergenic (‑157/‑285)	kilR ← / → sieB	killing protein, Rac prophage, FtsZ inhibitor protein/phage superinfection exclusion protein, Rac prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,418,386	0	G	A	100.0%	114.7 / NA	36	intergenic (‑157/‑285)	kilR/sieB	killing protein, Rac prophage, FtsZ inhibitor protein/phage superinfection exclusion protein, Rac prophage
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/16);  total (20/16)

TTTGTTGATGCGGTGTCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCGCATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCCCACCAT  >  NC_000913/1418256‑1418505                                                                                                                                   |                                                                                                                        tttGTTGATGCGGTGTCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACaa                                                                                                                 >  1:248524/1‑139 (MQ=255)     ttgATGCGGTGTCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATcaca                                                                                                                  <  1:708642/134‑1 (MQ=255)     ttgATGCGGTGTCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATcaca                                                                                                                  >  2:708642/1‑134 (MQ=255)                       ccTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGa                                                                                           >  2:873624/1‑139 (MQ=255)                       ccTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGa                                                                                           >  1:413002/1‑139 (MQ=255)                       ccTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGa                                                                                           >  1:361394/1‑139 (MQ=255)                        cTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGAc                                                                                          <  2:989156/139‑1 (MQ=255)                          ccAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACAt                                                                                        <  2:413002/139‑1 (MQ=255)                            aGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAAtt                                                                                                             >  1:470222/1‑116 (MQ=255)                            aGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAAtt                                                                                                             <  2:470222/116‑1 (MQ=255)                                   tGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCg                                                                                  >  2:536809/1‑136 (MQ=255)                                   tGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCg                                                                                  <  1:536809/136‑1 (MQ=255)                                   tGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAAc                                                                               >  1:567077/1‑139 (MQ=255)                                                    tACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGa                                                              >  2:222650/1‑139 (MQ=255)                                                         tCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccg                                                         >  1:298373/1‑139 (MQ=255)                                                         tCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccg                                                         >  1:406953/1‑139 (MQ=255)                                                               ttCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGa                                                   <  2:248524/139‑1 (MQ=255)                                                                         aaCCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAg                                         <  1:873624/139‑1 (MQ=255)                                                                              aTAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATg                                    <  1:354811/139‑1 (MQ=255)                                                                               tAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGt                                   >  1:183941/1‑139 (MQ=255)                                                                                 aGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCg                                 <  2:567077/139‑1 (MQ=255)                                                                                         cATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTAc                         >  2:187092/1‑139 (MQ=255)                                                                                            gACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAAc                      <  2:183941/139‑1 (MQ=255)                                                                                             aCTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGa                                                   >  1:1129814/1‑109 (MQ=255)                                                                                             aCTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGa                                                   <  2:1129814/109‑1 (MQ=255)                                                                                             aCTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAAcc                     >  1:766239/1‑139 (MQ=255)                                                                                              cTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCg                    <  1:359537/139‑1 (MQ=255)                                                                                                 ttAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAg                 >  1:1299959/1‑139 (MQ=255)                                                                                                        tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAggg                                                                                               <  2:1427415/54‑1 (MQ=255)                                                                                                        tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAggg                                                                                               >  1:1427415/1‑54 (MQ=255)                                                                                                        tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAACACTGATGCCGCCAGGa                                                   >  1:868165/1‑98 (MQ=255)                                                                                                        tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGa                                                   <  2:868165/98‑1 (MQ=255)                                                                                                        tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCgg                   >  2:1180905/1‑130 (MQ=255)                                                                                                        tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCgg                   <  1:1180905/130‑1 (MQ=255)                                                                                                               tGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCccacca   >  1:1111826/1‑139 (MQ=255)                                                                                                                gcgcTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCCCACCat  <  2:1299959/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                        TTTGTTGATGCGGTGTCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCGCATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCCCACCAT  >  NC_000913/1418256‑1418505

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