Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,301,917	C→T	Q246Y (CAG→TAT)	oppA →	oligopeptide ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,301,917	0	C	T	100.0%	281.3 / NA	74	Q246Y (CAG→TAT)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/38);  total (36/38)

GATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATA  >  NC_000913/1301677‑1302160                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     gatgatCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCt                                                                                                                                                                                                                                             <  1:536599/250‑1 (MQ=255)        acacaTTAGAAGTCAACTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCt                                                                                                                                                                                                                                     <  2:43727/251‑1 (MQ=255)          ataTTAGAAGTTATCTTAAGTGAACCCGTTCCGTACTTCTATCAATGCCTTGTTCACCCATCAACTTCACCGGTGCCAAAAGACGCTATCGAGAACTTCGGCGAAAAATGGACCCAGCCTGGTATTATCGTCACCAACGGTGCCTATACCTTCAACGATCGGGTCGTAAACGAACGAATCGTTATTGCCCGCAGCCCGACCTGCTGGAACAACGCGAAAACCGTTATTAACTATCTAACCTATTTGCCTa                                                                                                                                                                                                                                    >  2:484968/3‑250 (MQ=255)           caTTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTAtt                                                                                                                                                                                                                                  >  1:562202/1‑251 (MQ=255)                     cACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGa                                                                                                                                                                                                                          >  1:205998/1‑249 (MQ=255)                          tAAGTGAACCCGTTCCGTTCTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCCCCAACGGTGTCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCTCGGCCTACTGGAACAACGCGAAATCCGTTATTAAATATGTAACCTATTTGCCTATTGGTTCTGAAGTTAc                                                                                                                                                                                                                    <  2:92102/250‑1 (MQ=255)                          tAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAc                                                                                                                                                                                                                    <  2:201667/250‑1 (MQ=255)                                          gTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTAcc                                                                                                                                                                                                   >  1:398452/1‑251 (MQ=255)                                            aCTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGt                                                                                                                                                                                                 >  1:184529/1‑251 (MQ=255)                                                        tACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATAGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATc                                                                                                                                                                                     >  2:476501/1‑251 (MQ=255)                                                         aCTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCg                                                                                                                                                                                    >  1:566619/1‑251 (MQ=255)                                                               tCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGTTAATATCGGCTCCCACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGGTCTTGAACGCAGCCCCACCTACTAGAACAACGCGAAAACCGTTATTAACTATTTAACCTATTTGCCTATTGCTTCTGAAGTTACCGAGGTCAGCCGCTACCGTAGTGGGGAAATTGACATg                                                                                                                                                                               <  2:505686/250‑1 (MQ=255)                                                                 ccccATCAACCTCACCGGTGCCAAAAGCCGCTATCGAGAACTTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGTCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTTAACGCAGCCCGACCTACTCGAACAACTCGAAAACCGTTATTAACTATGTAACCGATTTGCCTATTGCTTCTTAAGTTACCGATGTCAACCGCTACCGGAGTGGTGAAATCGACATGACg                                                                                                                                                                            <  2:525082/250‑1 (MQ=255)                                                                 aCCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACg                                                                                                                                                                            >  1:50831/1‑251 (MQ=255)                                                                 aCCCAACAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACg                                                                                                                                                                            <  1:565282/251‑1 (MQ=255)                                                                    cATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAAt                                                                                                                                                                         <  2:132894/251‑1 (MQ=255)                                                                       cAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATTTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAAc                                                                                                                                                                      <  1:65585/251‑1 (MQ=255)                                                                            tCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACTGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGc                                                                                                                                                                   <  1:76785/249‑1 (MQ=255)                                                                                       aaaaGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAGCGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGTTTAAATCGACGTGTCGAATAACGGCATGCCGATCGAAt                                                                                                                                                      <  2:127810/251‑1 (MQ=255)                                                                                       aaaaGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAAt                                                                                                                                                      >  1:28550/1‑251 (MQ=255)                                                                                            ccGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTc                                                                                                                                                 <  2:50467/251‑1 (MQ=255)                                                                                                    gagaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATGGTTCCAGAAGCt                                                                                                                                         <  2:241496/251‑1 (MQ=255)                                                                                                     agaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTg                                                                                                                                        >  1:280452/1‑251 (MQ=255)                                                                                                     agaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGAGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTg                                                                                                                                        >  2:621411/1‑251 (MQ=255)                                                                                                        aaTTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGaa                                                                                                                                      >  2:77801/1‑250 (MQ=255)                                                                                                           tCGGCGAAAAATGGACCCAGCCTGGTAATAACGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGaaaaa                                                                                                                                   <  2:332054/250‑1 (MQ=255)                                                                                                                   aaaTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAATCTGAAAAAAGAGATccc                                                                                                                          <  1:128640/251‑1 (MQ=255)                                                                                                                     aTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCgg                                                                                                                        <  1:438049/251‑1 (MQ=255)                                                                                                                       ggACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGAc                                                                                                                      >  1:336597/1‑251 (MQ=255)                                                                                                                            cAGCCTGGTAATATCGTCACCAACGGGGACTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCGGCCCGACCTACTGGAACGAGGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGCAGTTACCAATGTCAACCGCTACCGTAGTGTGGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGCCGacgt                                                                                                                 >  1:202931/1‑248 (MQ=255)                                                                                                                                     aaTATCGGCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTg                                                                                                         <  1:353534/250‑1 (MQ=255)                                                                                                                                           ttCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCAt                                                                                                  >  2:395750/2‑250 (MQ=255)                                                                                                                                              tccAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCt                                                                                              >  2:202235/2‑251 (MQ=255)                                                                                                                                                aaaCGTTGCATATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCTACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGACATCGACATGACGAATAACAGCCTGCCGATCTAATTGTTCCAGAAGCTGAAACACGAGATCACGGCCGAAGTTCACGTTGAT‑CCATACCtgt                                                                                            >  2:590208/2‑251 (MQ=255)                                                                                                                                                   cGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCa                                                                                         <  1:262430/251‑1 (MQ=255)                                                                                                                                                        ccTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTAc                                                                                    <  2:170037/251‑1 (MQ=255)                                                                                                                                                            taCCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTAtt                                                                                >  2:571929/1‑251 (MQ=255)                                                                                                                                                             aCCTTAAAAGATTGGGACGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTa                                                                               <  2:223478/251‑1 (MQ=255)                                                                                                                                                               cTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACg                                                                             >  1:197886/1‑251 (MQ=255)                                                                                                                                                                  aaaaGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGaaa                                                                          <  1:433301/251‑1 (MQ=255)                                                                                                                                                                     aGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAAtt                                                                        >  1:565891/1‑250 (MQ=255)                                                                                                                                                                      gATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTaa                                                                      <  1:2000/251‑1 (MQ=255)                                                                                                                                                                      gATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTa                                                                       >  2:116459/1‑250 (MQ=255)                                                                                                                                                                        ttGGGTCGTAAACGAACGAATCGTTCTTGAACGCATCCCGACCTACTGGAACAACGCCAAAACCGTTATTAACTATGTCACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTTTTCCCGTAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGCT‑CCATACCTGTGCACTTACTATTACGAAAttcac                                                                     >  2:167698/1‑247 (MQ=255)                                                                                                                                                                            gTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCTCGTTGAT‑CCATGCCTGTGCACTTACTATTACGAAATTAACAACCa                                                                <  1:476433/251‑1 (MQ=255)                                                                                                                                                                            gTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAAcc                                                                 >  2:397200/1‑250 (MQ=255)                                                                                                                                                                                tAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGcc                                                                                                                                                              <  1:266401/154‑1 (MQ=255)                                                                                                                                                                                tAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGcc                                                                                                                                                              >  2:266401/1‑154 (MQ=255)                                                                                                                                                                                   acgaacgaATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCg                                                         <  2:552010/251‑1 (MQ=255)                                                                                                                                                                                    cgaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGc                                                        <  2:572816/251‑1 (MQ=255)                                                                                                                                                                                     gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGCACAACGCGAAAACAGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCACGGACGAAGTTCACGTTGAT‑CAATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGcc                                                       >  2:228594/1‑251 (MQ=255)                                                                                                                                                                                     gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTAt                                                                                 <  1:626483/225‑1 (MQ=255)                                                                                                                                                                                     gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTAt                                                                                 >  2:626483/1‑225 (MQ=255)                                                                                                                                                                                       acgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAACCCGTTATTAACTATGTAACCTTTTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTCCCGTAGTGGTGAAATCAACATGACGAATAACAGCATGCCGATCGAATTGTCCCAGAAGCTGAAAAAAGAGATCCCGGCCGAAGTTCACGTTGAT‑CCATCCCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCAt                                                     <  1:621411/251‑1 (MQ=255)                                                                                                                                                                                                    gAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcg                                        >  1:436287/1‑251 (MQ=255)                                                                                                                                                                                                    gAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcg                                        >  1:470330/1‑251 (MQ=255)                                                                                                                                                                                                    gAACGCAGCCCGACCTACGGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTAAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgc                                         <  1:592935/250‑1 (MQ=255)                                                                                                                                                                                                    gAACGCAGCCCGACCGACTGTAAACACGCGAAAACCGTTATTAACTAGGTAGCCTATTCGCCTATTGCTTCTGAAGTTCACGATGACAACCGCGACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGCTCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCTGACGAAGTCCACTTTGAT‑CCATACCAGTTCACTTCCAATAACGAAATTAACAACCATAACCCGCCATTCAACGAtgtgc                                         <  2:9330/250‑1 (MQ=255)                                                                                                                                                                                                       cGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcgtgt                                     >  1:193998/1‑251 (MQ=255)                                                                                                                                                                                                          aGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcgtgtgcg                                  >  1:50502/1‑251 (MQ=255)                                                                                                                                                                                                           gCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATgtgcgtgtgcgt                                 <  2:184529/251‑1 (MQ=255)                                                                                                                                                                                                             ccGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTAc                               <  1:543283/251‑1 (MQ=255)                                                                                                                                                                                                              cGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTAcc                              >  1:536839/1‑251 (MQ=255)                                                                                                                                                                                                                 ccTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGc                            >  1:450468/1‑250 (MQ=255)                                                                                                                                                                                                                    aCTGGAACAACGCAAAAACAGTTATTAACTATGTACCCTTGTTGCCTATTGCTTCTGAAGTGACGGATATCAACCGGTACCGTAGTGTTGAAAACGACATGACAAATAACAGCACGCCGGTCGAATTTTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTGGAT‑CCATAACTGTGCACTTTGTATTACGCAGTTAACATCCAAAAGCCGCCATTCAACGATGTGCGTGTGCGTACCGCACTg                        <  2:406765/251‑1 (MQ=255)                                                                                                                                                                                                                         aacaacGCGAAAACCGTTATTAACTAGGAAACCTATTGGCCTATTGCTTCTGAAGTTACCGAGTTCAACCGCTACCGTAGGGGTGAAAGCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCATTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCTCTGAAACt                   <  2:193998/251‑1 (MQ=255)                                                                                                                                                                                                                          acaacGCGAAAACCGTTATTAACTAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCACTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCGGAAAAAAGAGACCCCGGACGAATTCCACGTTGAC‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTa                  <  1:590208/251‑1 (MQ=255)                                                                                                                                                                                                                            cacGCGAAAACAGTGATTAACTATGTAACCTATTTGCCTATTGCTTCAGTAAATACCGATGTCAACCGCTACCGTAGTTGCGAAATCGACATGCCGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCTCGTTGATACAATAC‑TGTGCACTTCCTGTCAAGAAATTCACAACCAGAAGCGGCCATTCAACTATGTGCTTGTGCGTACCGCACTAAAACTAgg                >  2:432718/2‑251 (MQ=255)                                                                                                                                                                                                                            aacGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAgg                >  1:294479/1‑251 (MQ=255)                                                                                                                                                                                                                             acGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGt               <  2:623171/251‑1 (MQ=255)                                                                                                                                                                                                                               gcgAAAACCGTTATTAACTTTGTAACCTTTTTGCCTTTAGCTTCTTAAGTTCCCGATGTCAACCTCTCCCGGAGGGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCGGAAAAAAGAGATCCCAGTCGAATTTCGCGTTGAT‑CCATGCCTGTGCCCTCTCTATTACGAAATTAACAACCAGAAACCTCCATTCAACGATGTGCGTTTTCGTACCGCACTGAAAcgagtt               <  2:489684/249‑6 (MQ=255)                                                                                                                                                                                                                                    aaCCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTCCTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGAcc        >  1:223046/1‑251 (MQ=255)                                                                                                                                                                                                                                        gTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGAGGAGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGATGAT‑CCATACGTGTGGACTTACTATTACGAAATTAACAACTAGAAAGCGCCATTCAACGATACTCGTGTGGGAACCGCAGTAAAACTAGGGATGGACcgcg     >  2:95755/1‑250 (MQ=255)                                                                                                                                                                                                                                           attaACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTGCCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGata  <  1:263289/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     GATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAT‑CCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATA  >  NC_000913/1301677‑1302160

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: ‑