Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,200,881	G→A	A51V (GCC→GTC)	xisE ←	e14 prophage, putative excisionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,200,881	0	G	A	100.0%	166.9 / NA	53	A51V (GCC→GTC)	xisE	e14 prophage, putative excisionase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/25);  total (28/25)

AGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAG  >  NC_000913/1200754‑1201014                                                                                                                                |                                                                                                                                      aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGc                                                                                                                            >  2:411096/1‑139 (MQ=255) aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGc                                                                                                                            >  1:577386/1‑139 (MQ=255)  gCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                           <  2:1140569/139‑1 (MQ=255)   cGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCAc                                                                                                                          >  2:1065475/1‑139 (MQ=255)     gATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTc                                                                                                                        >  1:67284/1‑139 (MQ=255)        ataAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGcc                                                                                                                     >  1:820748/1‑139 (MQ=255)                  gACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGcc                                                                                                           >  2:133235/1‑139 (MQ=255)                  gACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGcc                                                                                                           <  1:133235/139‑1 (MQ=255)                     tGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                           <  1:25786/120‑1 (MQ=255)                     tGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCAAGCa                                                                                                                           >  2:25786/1‑120 (MQ=255)                      gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGt                                                                                                       >  1:654103/1‑139 (MQ=255)                       aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGtt                                                                                                      >  2:728488/1‑139 (MQ=255)                           attattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCAt                                                                                                  <  1:904721/139‑1 (MQ=255)                                    aCGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTATCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTgg                                                                                         <  1:178866/139‑1 (MQ=255)                                    aCGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                             >  2:737312/1‑103 (MQ=255)                                    aCGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                             <  1:737312/103‑1 (MQ=255)                                    aCGCGGTCTGGCTGCAATTACTCAATATCCGTTGCAAAGGTAATGAGTCATTTTTCTTGCTTACAGGTGATACCAATTCCCCAATTAAATAGACGACCCCACGCACTCGCGATAACCTGCCTTCTTTCACGGCCGGtcg                                                                                         >  1:513010/1‑137 (MQ=255)                                      gcgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                       <  2:654103/139‑1 (MQ=255)                                         gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                                    >  2:477080/1‑139 (MQ=255)                                         gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                                    >  2:1081335/1‑139 (MQ=255)                                          tCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                       >  1:72167/1‑135 (MQ=255)                                          tCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                       <  2:72167/135‑1 (MQ=255)                                          tCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaaa                                                                                   >  1:54717/1‑139 (MQ=255)                                           cTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaat                                                                                  <  1:993644/139‑1 (MQ=255)                                                tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAGCCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                             <  1:1398561/139‑1 (MQ=255)                                                tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                             <  1:728488/139‑1 (MQ=255)                                                    aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                           <  1:411878/89‑1 (MQ=255)                                                    aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                           >  2:411878/1‑89 (MQ=255)                                                     ttACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCt                                                                        <  1:411096/139‑1 (MQ=255)                                                      tACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCtt                                                                       >  2:1373149/1‑139 (MQ=255)                                                             aaTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAg                                                                     <  2:789818/134‑1 (MQ=255)                                                             aaTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAg                                                                     >  1:789818/1‑134 (MQ=255)                                                                   ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                          <  2:555828/139‑1 (MQ=255)                                                                             aaTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGaca                                                <  2:891790/139‑1 (MQ=255)                                                                                    cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgg                                         <  1:416692/139‑1 (MQ=255)                                                                                     aTTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgga                                        >  2:1164373/1‑139 (MQ=255)                                                                                                    aGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCAt                                                                                                  <  1:61531/66‑1 (MQ=255)                                                                                                    aGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCAt                                                                                                  >  2:61531/1‑66 (MQ=255)                                                                                                     ggTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAtt                        <  1:477080/139‑1 (MQ=255)                                                                                                        gtTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAttttt                     >  1:1222607/1‑139 (MQ=255)                                                                                                         ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTc                    >  2:287503/1‑139 (MQ=255)                                                                                                           aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag                  >  2:1163956/1‑139 (MQ=255)                                                                                                                cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                       >  1:556657/1‑65 (MQ=255)                                                                                                                cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                       <  2:556657/65‑1 (MQ=255)                                                                                                                cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                           <  2:241385/93‑1 (MQ=255)                                                                                                                cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                           >  1:241385/1‑93 (MQ=255)                                                                                                                cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc             <  2:590710/139‑1 (MQ=255)                                                                                                                cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc             >  2:1074418/1‑139 (MQ=255)                                                                                                                cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc             >  2:400364/1‑139 (MQ=255)                                                                                                                  ccccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCa           <  2:67284/139‑1 (MQ=255)                                                                                                                  ccccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCa           <  1:1164373/139‑1 (MQ=255)                                                                                                                  ccccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCa           <  1:1081335/139‑1 (MQ=255)                                                                                                                           actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAag  >  2:1216419/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                      AGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAG  >  NC_000913/1200754‑1201014

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