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%	252.4 / NA	63	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 (38/25);  total (38/25)

CATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCAT  >  NC_000913/2476387‑2476873                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     cATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCTCATAAAGACCGtaa                                                                                                                                                                                                                                               <  2:1050256/250‑1 (MQ=255)                 aaGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGc                                                                                                                                                                                                                              >  2:463253/1‑251 (MQ=255)                      aGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGcc                                                                                                                                                                                                                         <  2:1139306/251‑1 (MQ=255)                                       ccAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAgg                                                                                                                                                                                                        >  1:8209/1‑251 (MQ=255)                                               aTCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTg                                                                                                                                                                                                <  2:452422/251‑1 (MQ=255)                                                  cACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGtt                                                                                                                                                                                             <  1:802903/251‑1 (MQ=255)                                                   aCGGGAGAGAAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATACTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGttt                                                                                                                                                                                            >  1:535931/1‑251 (MQ=255)                                                          aGCAAGATGGTTATATAGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACAAGGAACTAGTGGTTTTATTCAt                                                                                                                                                                                     >  2:177290/1‑251 (MQ=255)                                                           gCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACt                                                                                                                                                                                                    <  2:979227/235‑1 (MQ=255)                                                           gCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGAAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACt                                                                                                                                                                                                    >  1:979227/1‑235 (MQ=255)                                                            cAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCAtt                                                                                                                                                                                    >  2:1094432/1‑250 (MQ=255)                                                                           gTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTAc                                                                                                                                                                    >  1:875825/1‑251 (MQ=255)                                                                            tGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCAtt                                                                                                                                                                                    >  1:681937/1‑234 (MQ=255)                                                                            tGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCAtt                                                                                                                                                                                    <  2:681937/234‑1 (MQ=255)                                                                                           cAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAg                                                                                                                                                    <  1:1008888/251‑1 (MQ=255)                                                                                                         tAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGTGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTATCCAGATAATAAACGCtt                                                                                                                                       <  2:535931/250‑1 (MQ=255)                                                                                                         tAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCt                                                                                                                                        <  1:286714/249‑1 (MQ=255)                                                                                                          aaGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTc                                                                                                                                      >  2:48811/1‑250 (MQ=255)                                                                                                              cgcgCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGaa                                                                                                                                                                                                                      >  2:1144412/1‑166 (MQ=255)                                                                                                              cgcgCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGaa                                                                                                                                                                                                                      <  1:1144412/166‑1 (MQ=255)                                                                                                                   cAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGc                                                                                                                            >  1:610055/1‑251 (MQ=255)                                                                                                                                                  cacaTCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGtaat                                                                                             >  2:1063421/1‑251 (MQ=255)                                                                                                                                                   acaTCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGtaat                                                                                             >  2:566942/1‑250 (MQ=255)                                                                                                                                                       cTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGc                                                                                        >  1:350131/1‑251 (MQ=255)                                                                                                                                                           aaatcacCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGAt                                                                                    >  2:996808/5‑251 (MQ=255)                                                                                                                                                                acCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCAt                                                                               <  2:875825/251‑1 (MQ=255)                                                                                                                                                                    ttCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTc                                                                            >  2:460936/1‑250 (MQ=255)                                                                                                                                                                     tCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCa                                                                           <  1:61778/250‑1 (MQ=255)                                                                                                                                                                     tCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCa                                                                           <  1:1125702/250‑1 (MQ=255)                                                                                                                                                                     tCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCAc                                                                          >  1:515312/1‑251 (MQ=255)                                                                                                                                                                     tCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCAc                                                                          >  2:403791/1‑251 (MQ=255)                                                                                                                                                                        aatcaatcCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGt                                                                       >  2:886238/1‑251 (MQ=255)                                                                                                                                                                          tcaatcCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACCCCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGCACTCGTGTTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGAGGTTATGTAATAATGCTGCTGGCATTTCACCGTc                                                                      >  2:132564/1‑250 (MQ=255)                                                                                                                                                                           caatcCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCat                                                                    >  2:969526/1‑251 (MQ=255)                                                                                                                                                                           caatcCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCat                                                                    <  1:566942/251‑1 (MQ=255)                                                                                                                                                                             atcCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCatat                                                                  <  1:463253/251‑1 (MQ=255)                                                                                                                                                                             atcCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCatat                                                                  <  1:460936/251‑1 (MQ=255)                                                                                                                                                                               ccccTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATc                                                                 <  2:123603/250‑1 (MQ=255)                                                                                                                                                                                      tgctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGCTGGCATTTCACCGTCATATCGCCAAACg                                                         >  2:4887/1‑251 (MQ=255)                                                                                                                                                                                       gctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGAAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACgg                                                        >  2:1041289/1‑251 (MQ=255)                                                                                                                                                                                       gctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAAGTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAAc                                                          <  2:8209/249‑1 (MQ=255)                                                                                                                                                                                        ctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAAt                                                                                                          >  1:189070/1‑200 (MQ=255)                                                                                                                                                                                        ctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAAt                                                                                                          <  2:189070/200‑1 (MQ=255)                                                                                                                                                                                        ctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGgc                                                       >  2:829327/1‑251 (MQ=255)                                                                                                                                                                                         ttCAAGGAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGgcg                                                      >  2:1114703/1‑251 (MQ=255)                                                                                                                                                                                                         cAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACt                                      >  2:100386/1‑251 (MQ=255)                                                                                                                                                                                                         cAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACt                                      >  2:871544/1‑251 (MQ=255)                                                                                                                                                                                                            gggACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGc                                   <  2:610055/251‑1 (MQ=255)                                                                                                                                                                                                             ggACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGCATAAGGAACACCAGGAACTCGTGGCTTTATTCATTATGAAGTGGTATTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGACTGCCCGGTAATGTGGTGATGTATTAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGc                                   >  1:513306/1‑250 (MQ=255)                                                                                                                                                                                                             ggACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGc                                   >  2:320716/1‑250 (MQ=255)                                                                                                                                                                                                              gACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGgcg                                                      >  2:1080206/1‑230 (MQ=255)                                                                                                                                                                                                              gACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGgcg                                                      <  1:1080206/230‑1 (MQ=255)                                                                                                                                                                                                               acacCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGc                                                                                                                                                       >  2:21529/1‑132 (MQ=255)                                                                                                                                                                                                                      ccTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCtttt                          <  2:350131/250‑1 (MQ=255)                                                                                                                                                                                                                             tttATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTAcc                   >  2:691890/1‑250 (MQ=255)                                                                                                                                                                                                                                   ccTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGtt            >  1:1066231/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaCCACATAAAGACCGTAAATTAAATTCGAGCTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGtttt          <  2:72250/248‑1 (MQ=255)                                                                                                                                                                                                                                        aaCCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGataat                                                                                                                                               >  2:640158/1‑115 (MQ=255)                                                                                                                                                                                                                                        aaCCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGataat                                                                                                                                               <  1:640158/115‑1 (MQ=255)                                                                                                                                                                                                                                        aaCCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGtttt          >  1:72250/1‑248 (MQ=255)                                                                                                                                                                                                                                        aaCCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttctt       <  1:1041289/251‑1 (MQ=255)                                                                                                                                                                                                                                            acaTAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTccc    >  1:854126/1‑250 (MQ=255)                                                                                                                                                                                                                                             caTAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCAt  >  2:728699/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     CATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCAT  >  NC_000913/2476387‑2476873

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