Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	281,581	G→A	R135W (CGG→TGG)	yagA ←	CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	281,581	0	G	A	100.0%	179.5 / NA	53	R135W (CGG→TGG)	yagA	CP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/27);  total (26/27)

TCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCC  >  NC_000913/281484‑281674                                                                                                  |                                                                                              tcCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTGCGGCGCGTCGTGTTCGAACCAGc                                                                                             <  1:201015/100‑1 (MQ=255)  cccgcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGAGGCCCATAAAATCCATCTGCCAGAGGCGGTTCGGCGCCTCCTGCTCCAACCAGcc                                                                                            >  2:470816/4‑100 (MQ=255)   cagcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGccc                                                                                           >  2:600947/1‑100 (MQ=255)     gcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGt                                                                                         >  1:187054/1‑100 (MQ=255)      cagGGTGGGGGGATGGCAGCGTCCACCGCCGGAGGGAAAGTGGCCCTTAGAAGCCATGTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTg                                                                                        <  2:445450/100‑1 (MQ=255)        gggTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGc                                                                                      >  2:482286/1‑100 (MQ=255)         ggTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCg                                                                                     <  1:292983/100‑1 (MQ=255)                ggATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGc                                                                              >  2:13671/1‑100 (MQ=255)                   tGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCg                                                                           >  1:551404/1‑100 (MQ=255)                     gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCGAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGt                                                                         >  2:875674/1‑100 (MQ=255)                      cAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAGTGGCCGGTg                                                                        >  2:707658/1‑100 (MQ=255)                      cAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTg                                                                        <  2:931744/100‑1 (MQ=255)                         cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAg                                                                     <  1:613714/100‑1 (MQ=255)                         cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAg                                                                     <  1:1412129/100‑1 (MQ=255)                            ccACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAgcgc                                                                  >  1:917406/1‑100 (MQ=255)                               ccGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCg                                                               <  1:133773/100‑1 (MQ=255)                                 gCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgc                                                             >  1:1123018/1‑100 (MQ=255)                                 gCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgc                                                             >  1:1159406/1‑100 (MQ=255)                                  ccAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgca                                                            <  1:1170424/100‑1 (MQ=255)                                   cAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcag                                                           <  2:204424/100‑1 (MQ=255)                                    aaaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcagc                                                          >  2:615635/1‑100 (MQ=255)                                        ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGc                                                      >  2:573386/1‑100 (MQ=255)                                        ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGc                                                      <  2:1042075/100‑1 (MQ=255)                                          aaaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCa                                                    >  2:524694/1‑100 (MQ=255)                                            aGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATg                                                  <  2:1279560/100‑1 (MQ=255)                                               ggCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCg                                               >  1:1131882/1‑100 (MQ=255)                                                 cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                             <  1:831079/100‑1 (MQ=255)                                                 cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                             >  2:211733/1‑100 (MQ=255)                                                  ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                            <  1:1293822/100‑1 (MQ=255)                                                  ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                            >  2:188349/1‑100 (MQ=255)                                                   cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGcc                                           >  1:821049/1‑100 (MQ=255)                                                    ttAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCa                                          <  1:805808/100‑1 (MQ=255)                                                    ttAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAACGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCa                                          <  2:974057/100‑1 (MQ=255)                                                       aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCa                                       <  2:1234468/100‑1 (MQ=255)                                                            cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGtt                                   <  2:1355807/99‑1 (MQ=255)                                                            cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTa                                  >  1:206389/1‑100 (MQ=255)                                                             aTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTAt                                 >  2:880810/1‑100 (MQ=255)                                                               cTGCCAGAGGCGGTTGGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATgg                               <  2:776287/100‑1 (MQ=255)                                                               cTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATgg                               <  2:953175/100‑1 (MQ=255)                                                               cTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATgg                               <  2:907334/100‑1 (MQ=255)                                                                tGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGa                              >  1:653204/1‑100 (MQ=255)                                                                 gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGTCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                             >  2:1065599/1‑100 (MQ=255)                                                                               ggCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg               <  2:1122636/100‑1 (MQ=255)                                                                               ggCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg               <  1:898520/100‑1 (MQ=255)                                                                                      cgtGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgt        <  2:482295/100‑1 (MQ=255)                                                                                       gtgtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg       <  2:1413136/100‑1 (MQ=255)                                                                                       gtgtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg       >  2:603506/1‑100 (MQ=255)                                                                                        tgtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtgt      <  1:524695/100‑1 (MQ=255)                                                                                         gtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg     >  1:1349110/1‑100 (MQ=255)                                                                                         gtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCAGGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg     <  1:1021283/100‑1 (MQ=255)                                                                                           tCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcc   >  2:1210746/1‑100 (MQ=255)                                                                                            cGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGccc  <  1:386137/100‑1 (MQ=255)                                                                                            cGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGccc  >  2:1505693/1‑100 (MQ=255)                                                                                                  |                                                                                              TCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCC  >  NC_000913/281484‑281674

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