Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,476,629	A→G	intergenic (+119/+65)	torI → / ← dsdC	prophage CPS‑53 recombination directionality factor and response regulator inhibitor/DNA‑binding transcriptional dual regulator DsdC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,476,629	0	A	G	100.0%	121.3 / NA	32	intergenic (+119/+65)	torI/dsdC	prophage CPS‑53 recombination directionality factor and response regulator inhibitor/DNA‑binding transcriptional dual regulator DsdC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/13);  total (19/13)

GATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGC  >  NC_000913/2476399‑2476860                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        gATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAAc                                                                                                                                                                                                                     <  2:310428/251‑1 (MQ=255)                       aaaGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGACGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCacac                                                                                                                                                                                               >  2:294981/1‑251 (MQ=255)                                gTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTc                                                                                                                                                                                      <  1:471579/251‑1 (MQ=255)                                            cGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTc                                                                                                                                                                          <  1:272505/251‑1 (MQ=255)                                             gAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTc                                                                                                                                                                          >  2:361013/1‑250 (MQ=255)                                                 aaGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTAt                                                                                                                                                                     >  2:481931/1‑251 (MQ=255)                                                         tataTCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTgg                                                                                                                                                              <  1:407200/250‑1 (MQ=255)                                                            atCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAATTGGttt                                                                                                                                                           <  2:367043/250‑1 (MQ=255)                                                                 gACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCt                                                                                                                                                     >  1:71524/1‑251 (MQ=255)                                                                   ccATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGt                                                                                                                                                   >  2:220817/1‑251 (MQ=255)                                                                                 aaaaTAAGCTCTTAAGCCGCGCCCATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTATCTTAGcc                                                                                                                                         >  2:170030/1‑247 (MQ=255)                                                                                 aaaaTAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGcc                                                                                                                                         >  2:378251/1‑247 (MQ=255)                                                                                                           gggTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCAct                                                                                                           <  1:432463/251‑1 (MQ=255)                                                                                                           gggTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCAAct                                                                                                           >  2:392630/1‑251 (MQ=255)                                                                                                                           aaaaTATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAAtg                                                                                            <  1:466856/250‑1 (MQ=255)                                                                                                                           aaaaTATTTCTCACATCTAAAAAACACCATTCCAATCAATCCACTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGAGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAAtgt                                                                                           >  2:230530/1‑251 (MQ=255)                                                                                                                             aaTATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATgtgg                                                                                         >  1:251592/1‑251 (MQ=255)                                                                                                                                 tttCTCACATCTAAAAAACACCATTCCAATCAATCCACTGCTGCTTCAAGTAGATGTCAGCAGGGGACACCATCCTTGTGTTTATCCCTAAAAACACATAAAGACAGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCACTTTTTGGCCACTGCCTGCCCGGTAATGTGGTGAt                                                                                     >  2:277483/1‑251 (MQ=255)                                                                                                                                      cacaTCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACGCGTGGTTTTATTCAt                                                                                                                                                                        <  2:383667/163‑1 (MQ=255)                                                                                                                                      cacaTCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGCAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAACACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCAt                                                                                                                                                                        >  1:383667/1‑163 (MQ=255)                                                                                                                                       acatctcaaatccacCATTCCAATCAATCACCAGCTGCTTCAAGTAGATGTCTGCAGGGGCCCACATTCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGATTTATTCATTATGAAGTGGTTTTTACCTATTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTTCCCGGTAATGTGGTGAttaaat                                                                                >  2:442194/13‑245 (MQ=255)                                                                                                                                               aaaaCACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGAt                                                                       >  1:30701/1‑251 (MQ=255)                                                                                                                                                 aaCACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGTTAATGTGGTGATGTAATAATGCTGATgg                                                                     <  2:310815/251‑1 (MQ=255)                                                                                                                                                   cacCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGc                                                                    >  1:385823/1‑250 (MQ=255)                                                                                                                                                                      ctgctgCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGcc                                                 >  1:342158/1‑250 (MQ=255)                                                                                                                                                                        gctgctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCa                                                <  1:400101/249‑1 (MQ=255)                                                                                                                                                                           gctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAAc                                             <  1:220817/249‑1 (MQ=255)                                                                                                                                                                             ttCAAGCAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGgcg                                         >  1:296343/1‑251 (MQ=255)                                                                                                                                                                                                     acCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGtaa                                                                                 >  2:227657/1‑187 (MQ=255)                                                                                                                                                                                                     acCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGtaa                                                                                 <  1:227657/187‑1 (MQ=255)                                                                                                                                                                                                       cATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCtt               >  1:240192/1‑251 (MQ=255)                                                                                                                                                                                                                    aTCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGc  <  1:277483/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        GATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGC  >  NC_000913/2476399‑2476860

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