Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,050,408	T→C	Y128H (TAT→CAT)	insB4 →	IS1 family protein InsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,050,408	0	T	C	100.0%	211.1 / NA	58	Y128H (TAT→CAT)	insB4	IS1 family protein InsB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/25);  total (33/25)

ACGCACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATCCCCCTCTGGAAATACGA  >  NC_000913/1050257‑1050701                                                                                                                                                        |                                                                                                                                                                                                                                                                                                      aCGCACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAAtacgta                                                                                                                                                   >  1:25832/1‑300 (MQ=14)    cACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTCCGc                                                                                                                                                >  1:581710/1‑300 (MQ=11)    cACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGc                                                                                                                                                >  1:13928/1‑300 (MQ=17)    cACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGc                                                                                                                                                >  2:420099/1‑300 (MQ=17)    cACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGAAGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGc                                                                                                                                                >  1:325845/1‑300 (MQ=17)            cGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATc                                                                                                                                        <  1:388407/300‑1 (MQ=21)                 cttgggCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGTGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGTTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAATTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGAt                                                                                                                                   <  2:150556/297‑1 (MQ=255)                 cTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAGTTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGAt                                                                                                                                   <  1:409854/300‑1 (MQ=37)                 cTGGGGCGTCTTATGAGCCGGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATTGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCGGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGAt                                                                                                                                   <  2:322653/300‑1 (MQ=255)                          cTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCGTTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACtt                                                                                                                          >  2:69986/1‑300 (MQ=255)                          cTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAAAATACGTACGCTTGAGATCAGGATTGTCTACtt                                                                                                                          >  2:111589/1‑300 (MQ=255)                           ttATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACttt                                                                                                                         >  1:57398/1‑300 (MQ=255)                             aTGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGa                                                                                                                       >  2:268181/1‑300 (MQ=255)                                     gctgTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCCCGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACa                                                                                                               >  1:469245/1‑300 (MQ=255)                                      ctgTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACAt                                                                                                              >  1:457663/1‑300 (MQ=255)                                            cccTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACANAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTtgt                                                                                                        >  1:532650/1‑300 (MQ=255)                                            cccTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTtgt                                                                                                        >  1:155729/1‑300 (MQ=255)                                                ttGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTTGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGaaa                                                                                                    <  1:229371/300‑1 (MQ=255)                                                ttGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGaaa                                                                                                    <  1:229933/300‑1 (MQ=255)                                                          gATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTc                                                                                          >  1:2081/1‑300 (MQ=255)                                                             atGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTc                                                                                       >  1:215112/1‑300 (MQ=255)                                                               ggATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACNTAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCCCAt                                                                                     >  2:655986/1‑300 (MQ=255)                                                               ggATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCAt                                                                                     >  2:413117/1‑300 (MQ=255)                                                                 aTGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATAc                                                                                   <  1:120757/300‑1 (MQ=255)                                                                 aTGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATAc                                                                                   <  2:391406/300‑1 (MQ=255)                                                                 aTGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAAAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATAc                                                                                   <  1:180545/300‑1 (MQ=255)                                                                  tgaaggATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTTTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATAcc                                                                                  <  2:580978/296‑1 (MQ=255)                                                                    aCGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGGAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGc                                                                                >  2:476413/1‑300 (MQ=255)                                                                         tggctggcCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCa                                                                           >  1:567358/1‑300 (MQ=255)                                                                          ggctggcCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCAc                                                                          <  1:328498/300‑1 (MQ=255)                                                                               gcCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACtatat                                                                     <  1:125453/300‑1 (MQ=255)                                                                               gcCGCTGGATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACtatat                                                                     <  2:28360/300‑1 (MQ=255)                                                                               accaatgcatggaTCCCGCCTGAAGGTAAAGCTGCACTTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTTAGGCAGCACCTGGCACGGCTGGGGCGGAAGTCGCTGTCGTTCTCAGCATCAGTGGAGCTGCACGACAAAGTCATCGTGCATTATCTGAACATAAAACACCATCGATAAGTTTGAGTCATTACCTACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTCTCCTCATAccccctccactcccc                                                                     <  2:28371/288‑14 (MQ=17)                                                                                 cGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTTGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATAc                                                                   <  1:59168/300‑1 (MQ=255)                                                                                 cGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATAc                                                                   <  1:422318/300‑1 (MQ=255)                                                                                 cGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCANGACAAAGTCATCGGGCATTATCTGAANATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATAc                                                                   <  2:516116/300‑1 (MQ=255)                                                                                 cGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCNGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCACGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATAc                                                                   <  2:659468/300‑1 (MQ=255)                                                                                  gCTGTATGAATCCCGCCTGAAGGGAAAGCTGNACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGNTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACt                                                                  >  2:538501/1‑300 (MQ=255)                                                                                  gCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTGTCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTNAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACt                                                                  >  2:557908/1‑300 (MQ=255)                                                                                  gCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACt                                                                  >  1:765000/1‑300 (MQ=255)                                                                                   cTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATACGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTg                                                                 >  1:307671/1‑300 (MQ=255)                                                                                   cTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTg                                                                 <  1:482507/300‑1 (MQ=255)                                                                                    tGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCNGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTNAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGc                                                                <  2:709973/300‑1 (MQ=255)                                                                                    tGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTTTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGc                                                                <  2:361775/300‑1 (MQ=255)                                                                                    tGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGc                                                                <  1:242194/300‑1 (MQ=255)                                                                                     gTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCg                                                               >  1:251260/1‑300 (MQ=255)                                                                                                      aGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCa                                              >  2:431792/1‑300 (MQ=255)                                                                                                             gCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACNAAGTCATCGGGCATTATCTGAACATAGAACACTATCAATAAGTTGGAGTCATTACCGACAANACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTACCTTTCCTGATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGaa                                       >  2:613648/1‑300 (MQ=255)                                                                                                                      aaTCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATaaaat                              <  1:289714/300‑4 (MQ=255)                                                                                                                       aTCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGAtt                             <  1:472589/300‑1 (MQ=255)                                                                                                                               gCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCNCTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGNTAACTa                     >  2:615895/1‑300 (MQ=255)                                                                                                                               gCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGTCGCGCAGAGGGAAATATAAGATTGTTAACTa                     >  1:451295/1‑300 (MQ=255)                                                                                                                               gCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTa                     >  1:643008/1‑300 (MQ=255)                                                                                                                               gCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTAACGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGCGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTCTATACTGCGGGTAGGAAGGGCGCGCAGAGGGGAATCTAAGATTGTTAACTa                     >  1:170257/1‑300 (MQ=255)                                                                                                                                 gATATACGCAGCGAATTGAGCGGCATAACCTNAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTNAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCAGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATc                   >  2:535770/1‑300 (MQ=255)                                                                                                                                             gAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGCGATCAGGGTTGTCTACTTTGAGTATAACATCCTTGTGCAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTCACTATCCACCCCTGGaaa       >  2:195646/1‑300 (MQ=255)                                                                                                                                               aTTGAGCGACATAACCTGAATCTTAGGCAGCACCTGGCACGGCTTGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTTAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATCCCCCTCTGGAAATa     <  2:4273/300‑1 (MQ=255)                                                                                                                                                  gAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTACCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATCCCACTCTGGAAATACGa  >  2:320600/1‑300 (MQ=255)                                                                                                                                                        |                                                                                                                                                                                                                                                                                                      ACGCACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATCCCCCTCTGGAAATACGA  >  NC_000913/1050257‑1050701

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