Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,515,876	C→T	G13D (GGC→GAC)	dam ←	DNA adenine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,515,876	0	C	T	100.0%	143.7 / NA	45	G13D (GGC→GAC)	dam	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/22);  total (23/22)

CAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCCTC  >  NC_000913/3515654‑3516115                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                cagtttGGTTCGGGAACAAACAGCTCTCGTGCGGGCTGTACGTACTCATCATTACGCATCTTCACAATGTGATAGAGGCTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCGCCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGAtt                                                                                                                                                                                                                     <  2:202223/248‑1 (MQ=255)   aTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATAGCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAGATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGAtttt                                                                                                                                                                                                                   <  1:48396/251‑1 (MQ=255)       gTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTttctt                                                                                                                                                                                                               <  1:549429/251‑1 (MQ=255)             ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCt                                                                                                                                                                                                         <  1:542326/251‑1 (MQ=255)               aacaaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGa                                                                                                                                                                                                       <  1:74785/251‑1 (MQ=255)                     caGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCAGGCGGactaact                                                                                                                                                                                                 <  1:158871/251‑1 (MQ=255)                                 ggCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTtc                                                                                                                                                                                     <  1:32551/251‑1 (MQ=255)                                        acgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACATTCTCCGGCt                                                                                                                                                                              >  2:24305/1‑251 (MQ=255)                                                    gTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCCGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCt                                                                                                                                                                  >  2:329639/1‑251 (MQ=255)                                                     tACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTc                                                                                                                                                                 >  1:274121/1‑251 (MQ=255)                                                             ttCACAATGTTATAGAGACGGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCCACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCCGTGGATACTGGTCCCCTGCCCACTGCAAAAAAGCGCGATTTTTCCACATGCTGCCTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGATCCGACAGc                                                                                                                                                          <  2:317768/250‑1 (MQ=255)                                                             ttCACAATCTTATAGAGACTGATCAGGTCGCTATTGATATCGGAAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCTAGCAGGCAATCACAAAAGGCTCAACCAGAGCTTCGACCTTGTGCAAAAGCCGATTAATATCATCAAGCAGGGAATACTTGTTCGCTGCCAAATTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAAGTAATTCCACATACTCCGGTGGTGGAGaacatcagggcagc                                                                                                                                                          >  2:448240/1‑238 (MQ=255)                                                              tCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCAt                                                                                                                                                        >  1:665829/1‑251 (MQ=255)                                                                                 gATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACACTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACCCCTTCTCCTGCTGTGGAGAAAGCTCCGACAGCATCCATCGCTTCAGGCAAc                                                                                                                                       <  2:505150/249‑1 (MQ=255)                                                                                         cGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATGCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTGCAAAAAAGCGCGATTTTTCGTCATGCTGACTGACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACGCCCGCAAc                                                                                                                              <  1:77851/250‑1 (MQ=255)                                                                                                 atatCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAAcc                                                                                                                      >  1:602510/1‑250 (MQ=255)                                                                                                   atCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAGTCGCCCTTGGGCAAATGCCGTTTAAGATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTGCGCCGGCTGTGGAGAAAGCTCCGACAGCAGCCAGCGCTTCAGGCATCCCCCGCATCGATAAACCAt                                                                                                                    <  2:450820/250‑1 (MQ=255)                                                                                                      ggCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCg                                                                                                                >  1:25210/1‑251 (MQ=255)                                                                                                            gATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCggg                                                                                                              >  2:537528/1‑247 (MQ=255)                                                                                                                 aaCGAAAAAAGTCGGTGTTGAGAAACACCCACCCGGCACCTCCACCAAGCTCAACCAGAAATTCACCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGAAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCACTGCGGGAGATAACCt                                                                                                     >  2:50375/1‑251 (MQ=255)                                                                                                                              ggTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCTCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGa                                                                                        <  1:315965/251‑1 (MQ=255)                                                                                                                                 gtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTc                                                                                           >  1:571734/1‑245 (MQ=255)                                                                                                                                        aaCACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGt                                                                               >  1:20260/1‑250 (MQ=255)                                                                                                                                                cccGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGc                                                                      >  2:610403/1‑251 (MQ=255)                                                                                                                                                    gCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCg                                                                  <  2:700396/251‑1 (MQ=255)                                                                                                                                                      aCCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGAGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGCTTACTTCAGATAGGCCTGTACCTGACGCAGCGGt                                                                >  2:250202/1‑251 (MQ=255)                                                                                                                                                        cTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAAc                                                                                                                                       >  1:61519/1‑178 (MQ=255)                                                                                                                                                        cTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTGCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCGCCGGCTGTGGAGAAAGCTCCGACAGCAGCCAGCGCGTCAGGCAAc                                                                                                                                       <  2:61519/178‑1 (MQ=255)                                                                                                                                                              aaGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCa                                                        <  1:49638/251‑1 (MQ=255)                                                                                                                                                                ggcgcAACCACACATTCGCCCACGGGCAAATGACGTTTCATATCATCAAGCAGGGGATACTTTTCACCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACATTCTCCGGCTGTGGCGCAAGCTCAGACAGCATCCAGCGCTTCAGGCAACTCCAGAAACGATAAACCATAGCGGGAGCTAACCTTGCTTCCTTCAGATCGGCATGTACCTGACGCAGCGGTTTCGCCCACg                                                      >  2:498538/5‑251 (MQ=255)                                                                                                                                                                    cAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGTGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTGTCTTCATGCTGACTAACAAATTACACCTTCTCCGTCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACCCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTCCCTGACGCAGCGGTTTCGCCCACggg                                                    <  1:480241/249‑1 (MQ=255)                                                                                                                                                                    cAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGt                                                   <  2:631798/250‑1 (MQ=255)                                                                                                                                                                     aaCCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATAAAACAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGATCCGACAGCATCCAGCGCTTAAGGCAACTCCCGCAACGATAAACCATAGCGGGAGATAAACTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGttt                                                 >  1:58260/1‑251 (MQ=255)                                                                                                                                                                           aCATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTc                                                                                                                                     <  1:371304/161‑1 (MQ=255)                                                                                                                                                                           aCATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGAGCGGCCTGTCCCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGcc                                           >  2:83348/1‑251 (MQ=255)                                                                                                                                                                           aCATTCGCCCTTGGGCAAATGACGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTc                                                                                                                                     >  2:371304/1‑161 (MQ=255)                                                                                                                                                                                 gCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCt                                                                           >  1:444645/1‑213 (MQ=255)                                                                                                                                                                                 gCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGAGCGGCCTGTACCt                                                                           <  2:444645/213‑1 (MQ=255)                                                                                                                                                                                        ggCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAGCCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTGCGCCCACGGGTTTTTGGCCTGGACATctgctg                              <  1:612525/251‑1 (MQ=255)                                                                                                                                                                                               gccttttAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCGTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACTATAATCCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATg                        <  1:370993/246‑1 (MQ=255)                                                                                                                                                                                                 cGTTTAATATCATCAAGCAGGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGt                       >  1:78109/1‑249 (MQ=255)                                                                                                                                                                                                    ttAATATCATCAAGCAAGGGATACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTCACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGATTTTGGCCTGGACATCTGATGGCAATGTAGAt                   >  1:613717/1‑250 (MQ=255)                                                                                                                                                                                                           catcaAGCATGGGATACTTTTCCCCTGCCCACTTCAAAAAAGCGCGATGTTTCTTCATGCTGACTAACTAATTACGCCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCttt           <  2:520775/251‑1 (MQ=255)                                                                                                                                                                                                             tcaAGCAGGGGATACTTGTCCCCAGCCCACTTCAAAAAAGCGCGATTGTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCttttt         >  1:209273/1‑251 (MQ=255)                                                                                                                                                                                                                    ggggGTACTTGTCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACATTCTCCGGCTGTGGAGAAAGCTCCGACAGCCTCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCctc  >  1:674392/1‑251 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                CAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCCTC  >  NC_000913/3515654‑3516115

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