Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,518,501	C→T	intergenic (‑15/+42)	aroB ← / ← aroK	3‑dehydroquinate synthase/shikimate kinase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,518,501	0	C	T	100.0%	152.3 / NA	46	intergenic (‑15/+42)	aroB/aroK	3‑dehydroquinate synthase/shikimate kinase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/21);  total (25/21)

ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTC  >  NC_000913/3518262‑3518743                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   atactccgggCCGTAAGGGAGGATAACGCGATCATGGTGAACCCCCGCCTGTTCAATTACGCCGCGGACCTTATCGAGATACAGAGGAGCAAGGGTTTCGGTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTGACGACAATCCGCTCCATAACGCTACATCCACTTTAATTACTGt                                                                                                                                                                                                                                         <  2:141857/242‑1 (MQ=255)       cTCGCCGTCAGGGAGGATAACGCTATCGACGTTACCACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGCTTTCAGCGGTAAGAATGAAGCTGGTTCATTCAACAAACCAGCTGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATACCGCGACATCCACTTTAATTACTGTACCCGc                                                                                                                                                                                                                                   >  2:134568/1‑251 (MQ=255)             gTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACg                                                                                                                                                                                                                              >  2:385299/1‑250 (MQ=255)                         aagccTATCGACGTTAACACCCGCCTGTTCAAGTCCGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTGTAATTACTGTACCCGCAGACGAGTGTATATAAAg                                                                                                                                                                                                                 <  2:474724/247‑1 (MQ=255)                                      ttAACCCCCGCCTGTTCAAGTACGCCTCGGACCTGATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTTGGTCATTAAACAAACCAGATGCGATGGTAATAGTGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTGAATTACTGTGCCCGCAGACGAGTGTATATAAAGCCAGAAGTAGtt                                                                                                                                                                                                     <  2:236807/250‑1 (MQ=255)                                           acCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTc                                                                                                                                                                                               >  2:507029/1‑251 (MQ=255)                                                   gTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCAt                                                                                                                                                                                         <  1:253877/249‑1 (MQ=255)                                                      cAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGa                                                                                                                                                                                     <  1:476520/250‑1 (MQ=255)                                                            cgccgcGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCt                                                                                                                                                                              >  2:314744/1‑251 (MQ=255)                                                            cgccgcGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATCACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCt                                                                                                                                                                              <  1:19510/251‑1 (MQ=255)                                                             gccgcGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTg                                                                                                                                                                             >  1:80130/1‑251 (MQ=255)                                                              ccgcGGACCTTATCGAGCTACAGAGGAGCAAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTgg                                                                                                                                                                            >  1:122612/1‑251 (MQ=255)                                                                 cgGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGttc                                                                                                                                                                         >  1:546396/1‑250 (MQ=255)                                                                  ggACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTg                                                                                                                                                                        >  2:367553/1‑251 (MQ=255)                                                                                  cAGAGGAGCTAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTTGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGTCAATCCTCTCCATAACGCGACATCCACTTTAATTCCTGTCCCCGCAGACGAGTGTATATAAATCCAGAATTAGCTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCt                                                                                                                                                        <  1:54276/251‑1 (MQ=255)                                                                                            aGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTCACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAATGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCCCTTTAGCGCTTTGAtcatca                                                                                                                                              <  2:16434/251‑1 (MQ=255)                                                                                                tttCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTAc                                                                                                                                          >  2:251490/1‑251 (MQ=255)                                                                                                           aCCAACATGACCTGCTCTCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCTCTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACg                                                                                                                               <  1:231903/251‑1 (MQ=255)                                                                                                               aCATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCgg                                                                                                                           >  2:277393/1‑251 (MQ=255)                                                                                                                   gACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCaa                                                                                                                        <  1:404846/250‑1 (MQ=255)                                                                                                                           cGCCCGATTTCATCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAATCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCATTACGAATGGTCACGTCGGCAATCTCTTCa                                                                                                                <  1:134568/250‑1 (MQ=255)                                                                                                                                gATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGACATTGGGTAACTACGTTCCCCGAGAGTCACGACAATCCTCTCAATAACGCGCCATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTACTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATAGTCACGTCGGCAATCTCTTCATACAGc                                                                                                          >  1:28261/1‑251 (MQ=255)                                                                                                                                        cGGTAAGAATGACGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATACACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAg                                                                                                                                                                                                                 >  1:442922/1‑140 (MQ=255)                                                                                                                                        cGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAg                                                                                                                                                                                                                 <  2:442922/140‑1 (MQ=255)                                                                                                                                        cGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCg                                                                                                  >  2:245783/1‑251 (MQ=255)                                                                                                                                              gAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCAtt                                                                                            <  1:368666/251‑1 (MQ=255)                                                                                                                                                     gCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGAAATCCACTTTAATTACTGTACCCGCAGACGCGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAAc                                                                                     >  1:405175/1‑251 (MQ=255)                                                                                                                                                       tGGTTCATTAAACCAACCAGATGCGATGGTAATTGGGTAACTACGTCCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGAATATAAAGCCAGACTTAGTTGATTTCCAGAATGTGACTAATCTGGTTTGCAACCACTTTAGAGCTTTACTCATCAGTACGAATGGTCACGTAGGCAATCTCTTCATACAGCGGCTTACGTTCATTGGCCACCGc                                                                                   >  1:375215/1‑251 (MQ=255)                                                                                                                                                        ggtgcatgaaacaaatCAGTTGCGACGGTAATGGGGTAACTACGTTCCCCTAGAGTATCGACAATCCTCTCCATAACGCGACATCCACTTTAAGTACTGTACCCGCAGACGAGTGTATAGAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTGTGATCATCAGTCCGAATGGGCACGTCGGCAACCTCTTCATACAGCGGGTTGCTTTCATTGGCCAACg                                                                                    <  2:28261/241‑1 (MQ=255)                                                                                                                                                          ttCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTc                                                                                >  1:570068/1‑251 (MQ=255)                                                                                                                                                             aTTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGc                                                                                                          <  2:252043/222‑1 (MQ=255)                                                                                                                                                             aTTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGc                                                                                                          >  1:252043/1‑222 (MQ=255)                                                                                                                                                                   caaacCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTc                                                                       >  2:161806/1‑251 (MQ=255)                                                                                                                                                                    aaacCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGAGACATCCCCTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCACTCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTc                                                                       >  2:360954/1‑250 (MQ=255)                                                                                                                                                                     aacCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCAc                                                                     <  1:425613/251‑1 (MQ=255)                                                                                                                                                                              cGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGgtg                                                            >  1:197949/1‑251 (MQ=255)                                                                                                                                                                                     aaTTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCaa                                                      <  1:339303/250‑1 (MQ=255)                                                                                                                                                                                         gggTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGGACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGAGTGCGTTCATTGGCCAACGCTTCCAGAACGTCACGCGGCGGTGTTTTAACGTg                                                  <  2:62894/250‑1 (MQ=255)                                                                                                                                                                                             aaCTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTgcagc                                              >  1:273097/1‑250 (MQ=255)                                                                                                                                                                                               cTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTgcagcacc                                           >  2:463420/1‑249 (MQ=255)                                                                                                                                                                                                   gTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGCGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCACAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGc                                       >  2:447281/1‑251 (MQ=255)                                                                                                                                                                                                                   gACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACGTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCtg                       <  1:268969/251‑1 (MQ=255)                                                                                                                                                                                                                     cAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCtgcg                     >  1:381914/1‑251 (MQ=255)                                                                                                                                                                                                                         ccTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCtgcgtgc                  <  2:58491/250‑1 (MQ=255)                                                                                                                                                                                                                                    cgcgACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCt      >  2:486645/1‑251 (MQ=255)                                                                                                                                                                                                                                        aCATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGGGCAAGTTGCTTTTc  <  1:463420/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTC  >  NC_000913/3518262‑3518743

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