Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,518,793	C→T	R89R (CGC→CGT)	insA7 →	IS1 family repressor protein InsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,518,793	0	C	T	100.0%	122.4 / NA	41	R89R (CGC→CGT)	insA7	IS1 family repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/17);  total (24/17)

ATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTA  >  NC_000913/4518547‑4519037                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      aTGTCTTTCCTGCCCCGCAACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGAAACAACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGTCGTCGGATGTCGCGCCAGTGCTCGCATTATGGGCGTTGGCTTCAACACGGTTGTACGTCACTTTAAAAACTCAGGCCGTAGTc                                                                                                                                                                                                                                                   <  2:229472/251‑1 (MQ=17)                    acggAAGGCGTGGTGCGTAACGGCAAAAGCCCTGCCGGACATCAGCGCTACCTCTGCTCTAATTGCCGTAAAACATGACAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACaa                                                                                                                                                                                                                                <  2:173355/248‑1 (MQ=17)                    aCTGAAGGCGTGGGGCGTAACGGCAAAACCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTGCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCTGTAACCTCTCGCATACa                                                                                                                                                                                                                                 <  2:129003/250‑1 (MQ=17)                     cTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAAc                                                                                                                                                                                                                               >  2:84395/1‑251 (MQ=255)                        aaGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCgg                                                                                                                                                                                                                            >  1:128522/1‑251 (MQ=255)                                 tgCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCCACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCCCCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCAAGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAgtgatg                                                                                                                                                                                                                   >  1:257381/1‑251 (MQ=255)                                   ggaaaCGGCAAAAGCACTGCCGGACATCAGCGCTATCCCTGCTCCAATGGCCGTAAAACATGGCAACTACAGTTAACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAgtgatgtg                                                                                                                                                                                                                 <  2:78807/248‑1 (MQ=17)                                      aaCGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGAtt                                                                                                                                                                                                              <  1:83785/251‑1 (MQ=255)                                                  aCTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGAGCCAGTGCACGCATTATGGGCGTGGGCCTCACCACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGGAACCTCGCGCATACAAcc                                                                                                                                                                                                                              >  1:237195/1‑223 (MQ=17)                                                         gACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGAGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGACGTAGTCGGTAACCTCGCGCATACAACCGGGAAGGGATGTGATTGTCTGCGCTGAAATGGACg                                                                                                                                                                                           >  2:34299/1‑251 (MQ=18)                                                           cATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGa                                                                                                                                                                                          >  1:97394/1‑250 (MQ=255)                                                                           tccctcattGCCGTAAAACAGGGCAACTCCAGTTCCCTAACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACGTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTc                                                                                                                                                                          <  2:97394/243‑1 (MQ=18)                                                                                      cGTAAAACATGGCAACTACAGTTCACTTACAACGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCCGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATcc                                                                                                                                                              >  1:242142/1‑250 (MQ=255)                                                                                                       aCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAACACTCAGGCCGTAGTCGGTAACCTCGCGCATACACCCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGACCAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCgctggctgct                                                                                                                                             >  1:244934/1‑249 (MQ=255)                                                                                                               cTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGt                                                                                                                                    >  1:244583/1‑251 (MQ=255)                                                                                                                       gCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAgg                                                                                                                            >  1:105694/1‑251 (MQ=255)                                                                                                                         ctctCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGAGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTCTGACAGGAt                                                                                                                          >  1:187071/1‑251 (MQ=255)                                                                                                                                 ccGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACggagg                                                                                                                   >  2:81272/1‑250 (MQ=255)                                                                                                                                     tacgCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGTGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGAGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATAAGGAGGACGGt                                                                                                              >  2:104815/4‑251 (MQ=255)                                                                                                                                        gCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTtgt                                                                                                           >  1:225399/1‑251 (MQ=255)                                                                                                                                                     aTTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCTGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTGCGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑GACGCGTATGACAGGAGACGGAGGACGGTCGTGGCGCACGTCTTc                                                                                              <  2:168887/251‑1 (MQ=255)                                                                                                                                                        gatgtGGCCATGAATGGCGTCGGATGTCGCGCCAGGGCGCGTATTATGGGCGTGGGCCTCAACGCGGTTGTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTTACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGGGGGTCTTGGTCGGTGCGGAATCACGTCAGTGTTGGGTATGT‑TACTGGGATGACAGGATACGGAGGACTTTTGGGGCGCagttattcggt                                                                                           <  2:15416/247‑11 (MQ=255)                                                                                                                                                        gATATGGCCATGAATGGCGTCGGATGTCGCGCCAGCGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCggg                                                                                           >  1:91922/1‑250 (MQ=255)                                                                                                                                                              gCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGc                                                                                     >  1:93221/1‑251 (MQ=255)                                                                                                                                                                 aTGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCCTACAACCGGGAAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACTGAGGGCGGTTGTGGCGCACGTCTTCGGTGAACGCAct                                                                                  >  2:65156/1‑251 (MQ=255)                                                                                                                                                                   gAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCActc                                                                                 <  1:207788/250‑1 (MQ=255)                                                                                                                                                                    aaTGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTg                                                                               <  2:225399/251‑1 (MQ=255)                                                                                                                                                                         cgtcggtgggcgcgCCAGTGCACGCATTAGGGGGGGTGGTCTCAACAGGGTTTTTCGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCAGACAACCGTGCAGTGATGTGAGTGACGGCTCTTAAATGGACGGATAGTTGTGGTACGTCGGTGCTAAATTGCGTCAGCGCTTGCTGTTT‑TACGCGTATGACAGGTTTCGTAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCcac                                                                          <  1:155967/241‑1 (MQ=255)                                                                                                                                                                         cGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCcac                                                                          <  1:173304/251‑1 (MQ=255)                                                                                                                                                                                   cgcgCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCg                                                                 <  1:152151/250‑1 (MQ=255)                                                                                                                                                                                        cAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGtcttct                                                           >  2:84049/1‑251 (MQ=255)                                                                                                                                                                                            gCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGATCGTCTTCTGAGc                                                       <  2:128522/251‑1 (MQ=255)                                                                                                                                                                                               cGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCctg                                                    >  1:172117/1‑251 (MQ=255)                                                                                                                                                                                                     aTGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCg                                              <  1:84395/251‑1 (MQ=255)                                                                                                                                                                                                         gCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCt                                          >  1:161750/1‑251 (MQ=255)                                                                                                                                                                                                           gTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGAAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGACGGCCttt                                        >  2:151660/1‑251 (MQ=255)                                                                                                                                                                                                                   cAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGGCAGCGGTGGC‑GTTTGTACGCGTATGACAGGATACGGAGGACGGTTGTTGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                 <  1:244601/250‑1 (MQ=255)                                                                                                                                                                                                                       acGGTTTTACGTCAATTGAAAACCTCCGGCCGTAGTCCCTAACCTCGCGCCTACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGCACAGTGGGGCTACGTCGGTGCTAAATCCCGTCAGCGATGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGCGCGTCTTCTGCGCCTGCTGTCGGCCTTTGAGGTCGTGGt                             >  2:129816/1‑250 (MQ=255)                                                                                                                                                                                                                                    aCTTAAAAAACTCAGGCCGTAGCCGGAAACCTCGCGCATACAACCGGACAGCGCTGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGGTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGGCGGCCTTTGCGGTCGTGGTATAGATGACGGAtg               >  2:116757/1‑251 (MQ=255)                                                                                                                                                                                                                                           aaaCTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGCGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGGATAGATGACGGATGGCTGGcc        >  2:252649/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aGGCCGTAGTCGGTAACCGCGCGCATACAACGGGGCAGTGATGTGATTGTCGGCGCTGAAATGGACGAACAGTTGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGGTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGGCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCTTGGTATAGATGACGGATGGCTGGCCGCTGTa  <  2:244934/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      ATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTT‑TACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTA  >  NC_000913/4518547‑4519037

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