Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,518,627	C→A	P34H (CCT→CAT)	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,627	0	C	A	100.0%	150.0 / NA	51	P34H (CCT→CAT)	insA7	IS1 family repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (36/15);  total (36/15)

AGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTA  >  NC_000913/4518381‑4518866                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                aGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGc                                                                                                                                                                                                                                             >  1:285727/1‑251 (MQ=255)  gCCTCTGCTTCTGGCATTTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGcc                                                                                                                                                                                                                                            <  2:171207/251‑1 (MQ=255)      ctGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGGAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAACAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTa                                                                                                                                                                                                                                         >  1:129714/1‑250 (MQ=255)         cTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAAc                                                                                                                                                                                                                                     >  1:38325/1‑251 (MQ=255)              ggCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGc                                                                                                                                                                                                                                >  1:501433/1‑251 (MQ=255)                      tCGGTCAGAATGAATGATGATATGGACTGGTATTCTTGACATCCTAACTGCACTAAAGGATGAGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCAATTTCAATCAGATGTCCTTCATGCTCCGCTACTGACGGCGTGGTGCGTAACGGCAAAAGCACTGACGGACATCAGCGCTATCACTGCTCTCATTGCCGTAAAACGTGGCCACTAc                                                                                                                                                                                                                          >  2:328761/1‑249 (MQ=18)                         gTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTc                                                                                                                                                                                                                     >  1:12127/1‑251 (MQ=255)                            aGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCACTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACt                                                                                                                                                                                                                  >  2:409456/1‑251 (MQ=255)                                  aCTCATGATGTGGTCTGCTATTATTGACATCCTCTCTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACAcc                                                                                                                                                                                                            <  1:17427/251‑1 (MQ=255)                                    tCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACCGTTCACTTACACCGc                                                                                                                                                                                                          >  2:445878/1‑251 (MQ=255)                                           gtgGTCTGCTATTTTTGACCTACTCACTGCCCTGACGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTCAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGCTGGCGTGGTGCGTAACGGAAACAGCACTGCCGGACATCACCGCTCTCTATGCTCTCATTGCCGTAAAACCTGGCAACTACAGTTCACTTCCACCGCCcctctg                                                                                                                                                                                                   >  2:256008/1‑249 (MQ=21)                                           gtgGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGCCATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCtctc                                                                                                                                                                                                     >  2:393799/1‑249 (MQ=255)                                                     attattGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTacgc                                                                                                                                                                                         >  2:174137/1‑248 (MQ=255)                                                      ttattGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGACAAATTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCACAAGCACTGACGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGACTCTCAGACCGGTacga                                                                                                                                                                                         >  2:226975/1‑247 (MQ=255)                                                               tCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAAt                                                                                                                                                                               >  1:337143/1‑251 (MQ=255)                                                                ccTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAAt                                                                                                                                                                               <  2:155763/250‑1 (MQ=255)                                                                        cccTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGatat                                                                                                                                                                      >  2:378346/1‑251 (MQ=255)                                                                        cccTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGatat                                                                                                                                                                      >  1:488526/1‑251 (MQ=255)                                                                              aGGATGGGGATTTCGGTAATGCTGCCAACTTACTGCTTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCAt                                                                                                                                                                >  1:113946/1‑251 (MQ=255)                                                                                                   gctcttccgatcTGATTTAGTGGATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACGCCGCCTCTCAGCCCGGTACGCACCAGAAAATAAATGATATGGCCATGAATGGCGTCGGATGTcaggc                                                                                                                                           >  2:130900/12‑247 (MQ=18)                                                                                                                 aTTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGCAAATCATTGATATGGCCATGAATGGCGTCGGATGGCGCGCCAGAGCACGCAtt                                                                                                                              >  2:277702/1‑250 (MQ=255)                                                                                                                               ggTGATTTTAAGGTGCTTGCGTGTCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATTAATGTCGTTGGATGCCGGGCAACCGCCCGCATTATGGGCGTTGGCCt                                                                                                                <  1:142990/250‑1 (MQ=11)                                                                                                                                       tAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGt                                                                                                       >  2:418207/1‑251 (MQ=255)                                                                                                                                             gCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCATCAACACGGTTTTACg                                                                                                 >  1:461965/1‑251 (MQ=255)                                                                                                                                                tGCGTGGCTTCCATTTCCATCAGATTTCCTTCCTGCCCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACAGGGCAACTCCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGt                                                                                                <  1:25033/249‑1 (MQ=255)                                                                                                                                                  cGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACt                                                                                            >  2:149074/1‑251 (MQ=255)                                                                                                                                                    tGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACtt                                                                                           >  2:411366/1‑250 (MQ=255)                                                                                                                                                    tGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGGTTTACGTCACtt                                                                                           >  2:207713/1‑250 (MQ=255)                                                                                                                                                               tccatCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAgg                                                                               >  1:179293/1‑251 (MQ=255)                                                                                                                                                                 catCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAgg                                                                               >  2:124520/1‑249 (MQ=255)                                                                                                                                                                 catCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAgg                                                                               <  1:124520/249‑1 (MQ=255)                                                                                                                                                                      gctGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGACCATGAATGGCGTCGGATTTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAACAAACTCAGGCCGTAGt                                                                        >  2:487692/3‑251 (MQ=18)                                                                                                                                                                        tGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGAGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGCTATGGCCATGAATGGCGTCGGATGTCGCGCCCGTGCACGCATTATGGGCGTTGGCATCAACCCGGTTTTAAGTCACTTAACAAACTCAGGCCGTAGTCg                                                                      >  2:400151/1‑251 (MQ=18)                                                                                                                                                                            cttcctGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGCCATCAGCGCTATCTCTGCTCTCATAGCCGTAAACCATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATAATTGATATGGCCATGAATGGCGTCGTATGTCGCGCCAGTGCACGCCTTATGGGCGTTGGCCTCAACACGGTTTTACGTCCCTTAAAAAACTCAGGCCGTAGTCGGTa                                                                   >  1:474378/1‑250 (MQ=255)                                                                                                                                                                            cttcctGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCCGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCCTTGTTATGGCCATGAATGGCGTCGGATGTCGCGCCAGAGCCCGCATTCTGGGCGTTGGCCTCGACACGGTATTACGTCACTTAAAACACTCAGGCCGTAGACGGTa                                                                   >  2:272014/1‑250 (MQ=17)                                                                                                                                                                             ttcctGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACACCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTa                                                                   <  1:393799/249‑1 (MQ=255)                                                                                                                                                                                 tGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTcg                                                             >  1:304860/1‑251 (MQ=255)                                                                                                                                                                                             gagggcgtgtgtcgaaaCGACAAAAGCCCTCCCGAACACCCGCGCTATCTCTGCTCTCATTGCCGAAAAACAGGGCAACTCAAGTTCCCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAAc                                                   <  2:285727/234‑1 (MQ=17)                                                                                                                                                                                             gAAGGCGTGGTTCGTAACGGCAAAAGCACTGCCGCAAACCAGCGCTATCTCTCCTCTCACCGCCGTAAAACTTGGCAACTACAGTTCACTTATCCCGCCTCTCAGCCCGCTACGCACCAGAAAATCATTGATATGTCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCg                                                 <  1:143141/251‑1 (MQ=255)                                                                                                                                                                                                ggCGTGGTGCTAAACGAAAAAATCACTCCCGGACATCCGCGCTATCTCTGCTCCCATTGCCGTAAAACGGGGCAAATACAGTTCACTTACCCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCTCTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGc                                              <  2:129714/251‑1 (MQ=17)                                                                                                                                                                                                      gtgCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTAGCGCCAGTGCACGCATTCTGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAAATCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAgtgat                                        >  2:461387/1‑251 (MQ=255)                                                                                                                                                                                                         cGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTCCGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAgtgatg                                       >  1:207950/1‑249 (MQ=255)                                                                                                                                                                                                               ggCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTc                               >  2:459227/1‑251 (MQ=255)                                                                                                                                                                                                                gCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACGACAGTTCACTTACACCGCCTCTCAGCCCGGGACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGACAGTGCACGCATTATGTGCGTTGGCCTCAACACGGTTTTACGACACTTAAAACACTCAGGCCGTAGTCGGTAACCTCGCACATACCACCGAGCAGGGATATGATTGTCt                              >  1:149026/1‑251 (MQ=17)                                                                                                                                                                                                                  aaaaGCACTGCCGGAAATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTgc                            <  1:445878/251‑1 (MQ=255)                                                                                                                                                                                                                   aaaGCACTGCCGAAGATCAGCCCTTGCCCTGCTCTCATGGCCGTAAAACATGGCACCTACAGTTCCCTTCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTgcg                           <  2:180956/251‑1 (MQ=18)                                                                                                                                                                                                                              gcaCATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACAACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCAACAGAAAATAATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACAAGGTTTTACGTCACTTCAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAGTGGAc                >  2:357526/3‑251 (MQ=255)                                                                                                                                                                                                                                     aGCGCTATCTCTGCTCCCATTGCCGAAAAACATGGCAACTACATTTCCCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGt         <  1:300614/251‑1 (MQ=18)                                                                                                                                                                                                                                          tATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGc    >  1:397448/1‑251 (MQ=255)                                                                                                                                                                                                                                            cctctGCTCTCATTGCCGAAAAACATGGCAACTCCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTa  <  1:409456/250‑1 (MQ=255)                                                                                                                                                                                                                                               ctGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTa  <  1:174137/248‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                AGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTA  >  NC_000913/4518381‑4518866

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