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%	137.7 / NA	44	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 (16/28);  total (16/28)

ACCTTCATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATC  >  NC_000913/2476382‑2476874                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      aCCTTCATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGTGTAAAATATTTCGCACATCTAAAAAACACCATGCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTTTGGTTATCCCTAAGACCACATAAAGAcc                                                                                                                                                                                                                                                    <  2:541277/251‑1 (MQ=255)           aTGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGtaaattaaa                                                                                                                                                                                                                                          <  2:558226/251‑1 (MQ=255)                 ggATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGaa                                                                                                                                                                                                                                    <  1:269006/251‑1 (MQ=255)                   agtAAGTCAGGCGACCTTCCTAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCTCGCCAATGGGTAAAAGAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATTCTTGTGTTTATCCCTAAACCCTCATAAAGACCGTAAATTAAATTCGAAca                                                                                                                                                                                                                                  <  2:520241/249‑2 (MQ=255)                               gACCTGCCAAAAGCCAAAGTTATCCACGGGCNAGCAAGATGGTTATATCGGGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTCGCAGGCCGAAt                                                                                                                                                                                                                      <  2:282854/251‑1 (MQ=255)                                    gCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAAACACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAg                                                                                                                                                                                                                   <  2:77537/249‑1 (MQ=255)                                     ccAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGca                                                                                                                                                                                                                <  1:448522/251‑1 (MQ=255)                                        aaaGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCacac                                                                                                                                                                                                             >  1:106451/1‑251 (MQ=255)                                                 gTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTc                                                                                                                                                                                                    <  2:321591/251‑1 (MQ=255)                                                 attatacccaggCGAGCAAGATGGATATATCGTTACCATTGTGAATTCAAAAATAAGCTCGTAAGCCTCGACAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTACAAAACCCCATTCCAATCAATCCCCTGCTTCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAATCCACATAAAGACCGTAAATTAAAGTCGAACTAGCAGACCGAATAAGGctcaacactaactc                                                                                                                                                                                                    <  2:307265/241‑14 (MQ=255)                                                    aTCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTg                                                                                                                                                                                                 >  1:284292/1‑251 (MQ=255)                                                          gggCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGCAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGAAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTCAATTAAATTCGAACTAGCAGGCNGAATAAGGCACACCAGGAACTCGTGGttttt                                                                                                                                                                                           >  2:283310/1‑250 (MQ=255)                                                                   aGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCAGCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATg                                                                                                                                                                                  <  2:52295/251‑1 (MQ=255)                                                                            tatCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTAACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGACCAAGCGGTTTTATTCATTATGAAGTGGttt                                                                                                                                                                         >  2:408350/1‑251 (MQ=255)                                                                             atCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGtttt                                                                                                                                                                        >  1:343172/1‑251 (MQ=255)                                                                                  gACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAACGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCt                                                                                                                                                                   >  2:406725/1‑251 (MQ=255)                                                                                  gACCATTGTGAATTCAAAAATAAGCGCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCTCATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTAcc                                                                                                                                                                    <  2:442181/250‑1 (MQ=255)                                                                                   aCCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTg                                                                                                                                                                  <  1:251053/251‑1 (MQ=255)                                                                                                  aaaaTAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGat                                                                                                                                                   <  2:501662/251‑1 (MQ=255)                                                                                                  aaaaTAAGCACTTAATCCGCGCCACTGGATAAAATAGCGAGCAAAATATTTCTCAAATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATATCGGCAGGGGACACCATGCTTGAGTTTGTACATAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGCGGTTTTATTCCTTATGTAGTGGATTTTACCTTTTCTCTTAACCAGa                                                                                                                                                    >  2:370018/1‑250 (MQ=255)                                                                                                              tAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTTCTTCAATTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGTCCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGTAACTCGTTGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCtt                                                                                                                                        <  1:23598/250‑1 (MQ=255)                                                                                                                         aaTGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGcc                                                                                                                            >  2:444126/1‑251 (MQ=255)                                                                                                                                 aaaTAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGa                                                                                                                                                                                                                                     >  1:151258/1‑138 (MQ=255)                                                                                                                                 aaaTAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGa                                                                                                                                                                                                                                     <  2:151258/138‑1 (MQ=255)                                                                                                                                            aaaaTATTTCTCACATCTAAAAAACACCATACCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGGCCGCAAATTAAATTCAAACTAACAGGCCCAATCAGGAAAACCAGGAACTCGTAGTTTTATTCATTATGAAGTTGTTTTTATCTGTTCTCTTAGCCAGATAATAATCGCTTCACATTTTGGCCACtgcctgccaggaattgg                                                                                                         >  2:210218/1‑242 (MQ=255)                                                                                                                                                    tctcACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCGCATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGTCACTCCAGGAACTCGTGGTTTTATTCATTAGGAAGTGGTTTTTACCTGTTCTCTTAGCTAGATAATAAACGCTTCAATTTTTGTCCACTGCCTGCCCGGTAATGTGTTGATg                                                                                                  <  2:550340/250‑1 (MQ=255)                                                                                                                                                    tctcACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGACAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATg                                                                                                  <  1:318085/250‑1 (MQ=255)                                                                                                                                                     ctctcaTCTAATAAACACCATTCCAATCAATCCCCTGATGCTTCAAGTAGATGTCTGCAGGGGACACAATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAAATCGTGGTTTTCTTCATTATGAAGTGTTTTTTACCTGTTCTCTTAGCCAGTTAATAAACGCTTCAATTTTTGACCACTGCCTGCCCGGTAATGTAGTGCTGt                                                                                                 >  2:305076/5‑250 (MQ=255)                                                                                                                                                       cacaTCTAAAAAACACCATTCCAATCAATCCCCTGCTGCGTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGtaa                                                                                               <  1:457459/250‑1 (MQ=255)                                                                                                                                                            cTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGCCACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGc                                                                                         >  2:469738/1‑251 (MQ=255)                                                                                                                                                                aaaaCACCATTCCAATCAAGCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGAACGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGAt                                                                                     <  1:560005/251‑1 (MQ=255)                                                                                                                                                                                caatcCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCat                                                                     >  1:237987/1‑251 (MQ=255)                                                                                                                                                                                  atcCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCatat                                                                   <  1:194790/251‑1 (MQ=255)                                                                                                                                                                                    cccatgctgcATCAAGTAAATGGCTGCAGGGGACACCATCCTTGCGTTTATCCCTAAAACCACATAAAGAACGTAAATTAAATTCGAACTAGCATGCCGAATAAGTCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGATCTCTCAGCCGGATAATAAACGCTTCAATTTTAGGCCACTGACTGCCCGGTAATGTGGTGTTGTAATAATGCTGATGGCATTTCACCGTCCTATCg                                                                 >  2:488609/5‑251 (MQ=255)                                                                                                                                                                                            gctTCAAGTAGATGTCTGCAGGTGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACgg                                                         <  2:140269/251‑1 (MQ=255)                                                                                                                                                                                                 aaGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGcgcg                                                     >  1:276692/1‑250 (MQ=255)                                                                                                                                                                                                             gCAGGGGACACCATCCTTGTGTTGATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACCAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTAGTCATTATTACGTGGTTTTTCCCTGGTCTCCTAGCCTGATAATAAACTCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGGGATGTAATAATGCTGATGGCATTTCGCCGTCATATCGACAAACGGCGCGACGAGCTccccac                                        <  2:226843/251‑1 (MQ=255)                                                                                                                                                                                                                ggggACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCTAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTTGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTgg                                     <  2:364958/251‑1 (MQ=255)                                                                                                                                                                                                                         aTCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCttt                            >  1:84910/1‑251 (MQ=255)                                                                                                                                                                                                                           ccTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTGTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCttttt                          <  1:222034/251‑1 (MQ=255)                                                                                                                                                                                                                              tgtgTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGTCGAGCTCACCACTGTCAAGCCTTTTTTg                        <  1:465116/250‑1 (MQ=255)                                                                                                                                                                                                                                             aaCCTCATAAAGACCGGAAATGAAATTCGAACTGGCTGGCCGAATTAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTTGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGCCGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttctt        <  2:91098/251‑1 (MQ=255)                                                                                                                                                                                                                                             aaCCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttctt        <  2:119068/251‑1 (MQ=255)                                                                                                                                                                                                                                                   aTAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCAGTATTAAGTGGTTTTTTCCTGTTCTCTTATCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATc  <  2:468135/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      ACCTTCATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATC  >  NC_000913/2476382‑2476874

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