Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	466,148	G→A	P389S (CCA→TCA)	ybaE ←	uncharacterized protein YbaE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	466,148	0	G	A	100.0%	171.7 / NA	52	P389S (CCA→TCA)	ybaE	uncharacterized protein YbaE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/29);  total (23/29)

GCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTC  >  NC_000913/465913‑466381                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          gCTTCGCCATTTAATCTGTCGCCCATCATGAGGTTTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGAc                                                                                                                                                                                                                             <  2:487734/250‑1 (MQ=255) gCTTCGCCAATTAATCTGTCGCCCAACATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACt                                                                                                                                                                                                                            >  1:4772/1‑251 (MQ=255)   ttCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGAGTCGTGTCGTCCCAGTTTTTTGCGTTATGCAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGACTAGTCCAGCCTGACAGTAATGCATGACTgg                                                                                                                                                                                                                          >  2:72387/1‑251 (MQ=255)        aaagtAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCGGTAGGGTCGTGTCGTCCCATTTTTGTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGGCGCCGTTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATTCGTTACCTGCCAGTGCGGAATAGTCCAGCTTGACCGTATTGCATGAATGGCGGt                                                                                                                                                                                                                      <  2:211353/246‑1 (MQ=255)         aaTTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGa                                                                                                                                                                                                                    <  2:162600/251‑1 (MQ=255)             aaTCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTggtc                                                                                                                                                                                                                  >  1:404819/1‑246 (MQ=255)                  gTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGttt                                                                                                                                                                                                            <  2:437783/250‑1 (MQ=255)                     gCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTtctc                                                                                                                                                                                                        <  2:530565/251‑1 (MQ=255)                      cccATAATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGGTCAGGGTTTCTcc                                                                                                                                                                                                       >  2:554651/1‑251 (MQ=255)                         atcatGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCACGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGCCTGGCGGTGATCAGGTTTTTTCCGAc                                                                                                                                                                                                    >  1:520520/1‑251 (MQ=255)                              gAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGCGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGGTCAGGGTTTTTcccgctttt                                                                                                                                                                                                >  2:289976/1‑243 (MQ=255)                                          tgCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGtaat                                                                                                                                                                                    >  1:5480/1‑250 (MQ=255)                                               gTAGGGTCGTGTCGTTCCAGTTTTTTGCGTTATGAAAAATAATTGGGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGGTCTGCCATGGTATGAATTTCTATCGGTAGGTGATAAACCAGCGTTAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAAc                                                                                                                                                                                <  2:394838/249‑1 (MQ=255)                                                tagtgtagtgtcgtCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGAAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGATTTCTTCGCATTCTAACGGTTGTATTAAAcc                                                                                                                                                                               >  2:212231/8‑249 (MQ=255)                                                            gtcCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAAt                                                                                                                                                                                                                                    >  2:470083/1‑184 (MQ=255)                                                            gtcCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAAt                                                                                                                                                                                                                                    <  1:470083/184‑1 (MQ=255)                                                                      tttGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCAc                                                                                                                                                       >  1:204544/1‑251 (MQ=255)                                                                         gCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCTCTGCCAGTGTCGCCTTTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCTGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACttt                                                                                                                                                    <  1:74031/251‑1 (MQ=255)                                                                                          tgtgAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCCGGGTTTATCCGACTTCTAACGTTTGTACTAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGcccgcggg                                                                                                                                   >  1:283775/1‑245 (MQ=255)                                                                                                  cacaGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCCGCCAGAGGGAGAGTCgg                                                                                                                           >  2:187851/1‑251 (MQ=255)                                                                                                  cacaGCATTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGGATAATGGAGATCACTTTTTGCGCCTGCCAGAAGGAGAGGCgg                                                                                                                           >  2:548146/1‑251 (MQ=255)                                                                                                     aGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGa                                                                                                                        <  2:4772/251‑1 (MQ=255)                                                                                                       ccTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACt                                                                                                                      <  2:61767/251‑1 (MQ=255)                                                                                                           ccGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTg                                                                                                                  <  1:554651/251‑1 (MQ=255)                                                                                                           ccGCGGCCAGTGTCGCCTGTAGGCGTGCTGCCATGTTATGAATGTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTCCCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCGCTTTTCGCGCCTGCCAGAGGTAGAGTCGGGGACtttt                                                                                                                   <  2:542933/250‑1 (MQ=255)                                                                                                                 ccAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAAc                                                                                                            <  1:97771/251‑1 (MQ=255)                                                                                                                   ggtgtCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTGTGCGCAACGt                                                                                                          <  2:219171/250‑1 (MQ=255)                                                                                                                   aGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGt                                                                                                          >  1:95116/1‑251 (MQ=255)                                                                                                                             gTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCaa                                                                                                <  2:404819/251‑1 (MQ=255)                                                                                                                                   gTTCTGCCATGGAATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTACTGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAAAGTTTGTAATAAACCGGATTGTTGAATAATGGAGATCACTTTTCGCGCCCGCCCGAGGGAGAGTCGGGGGATTTTGCGCAACGTCAAATAGGAAAAACCt                                                                                          >  1:576796/1‑251 (MQ=255)                                                                                                                                      cTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCACGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCCACTTCTAACGTTTGTATTAAACCGGAATGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGAATTTTGCGCAACGTCCAATAGCAAAAACCTaa                                                                                        >  1:282933/1‑250 (MQ=255)                                                                                                                                        gCCATGGTATGAAGTTATATCGGTAGGTGAGAAACCGGCGTCAAGGTTTTCGGTAGTTGGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTATCGTTTGTAATAAACCGGATTGATGAAGAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAAGAGCAAAATCCTAAAc                                                                                      <  2:204342/250‑1 (MQ=255)                                                                                                                                             ggTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGc                                                                                <  1:104132/251‑1 (MQ=255)                                                                                                                                                  gAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCt                                                                            <  1:388603/250‑1 (MQ=255)                                                                                                                                                  gAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACGTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTa                                                                           <  1:541239/251‑1 (MQ=255)                                                                                                                                                       tCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGa                                                                       <  1:441672/250‑1 (MQ=255)                                                                                                                                                       tCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCCGAGGGAGAGTCGGGGGCTTTTTCGCAACGTCAAATAGCAAAAACCTAAACCGATGCCGCttcttgc                                                                      >  2:114679/1‑245 (MQ=255)                                                                                                                                                          aTCGGTAGGTTATAAACCAGCGTCAAGGTTTTCGGTATTTTGACTTCATCCGTTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGGTTGTAATAAACCGGATTGATGAATAATGGAGATCTCTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGTCTTTTGCGCAACTTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCt                                                                    <  1:256055/250‑1 (MQ=255)                                                                                                                                                          aTCGGTAGGTGATAAACCAGCTTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCt                                                                    <  1:72387/250‑1 (MQ=255)                                                                                                                                                                                tCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGtttcaagccctcc                                             >  1:363146/1‑240 (MQ=255)                                                                                                                                                                                                 ttGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAAATCACTTTTCGCGCCTGCCCGAGGGAGCGTCGGGGACTTTTTCGCCACGTCAAATAGCAAAAACCTAAAATGATGCCCCTCCTGACCTGGCTGACCCGTTGCAGCCCCTTCGGTTTTCCGATgggggt                            >  1:261750/1‑247 (MQ=255)                                                                                                                                                                                                       tCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGGGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCACAGTCAAAAAGCAAAAAACTAAAATGGTGCCGCTACTGACCTGGCTGACCCGTTGCAGCCCCTCAGGTTTGACGGTGGGGGTTTGCAc                      >  1:134611/1‑251 (MQ=255)                                                                                                                                                                                                        cATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCCGGTTTTCTCCGACTTTTAACGTTTGTAATAAACCGGATTGATGAATAATGGGGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGGATTTTGCACAATGGCAAAAAACAAAAACCTAAACTGATGCCGCTACTGGCCTTGCGGAACCGGTTTAGCTCCCCCAGTTTTCCGATggggtttttcacg                     >  1:226660/1‑240 (MQ=255)                                                                                                                                                                                                        cATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCAc                      <  1:289976/250‑1 (MQ=255)                                                                                                                                                                                                        cATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGGTTGTAATAAACCGGATTGATGAATAATGGAGATCACGTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATGTGCACg                     <  2:261750/251‑1 (MQ=255)                                                                                                                                                                                                        cATCCGGTACCTGCAAATGCGGAATAGGCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCGCCTCCGGTTTGCCGCTGGTGATTTGCAc                      <  2:100506/250‑1 (MQ=255)                                                                                                                                                                                                         ctCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGATTTTTCCGACTTCTAAAACTTTTAATACACCGGATTGATGGATAATTGAGATCCCTTTTCGCCCCTGACAGAGGGAGAGGCGGGGACTTTTTAGCCACGTACAATAGCCAAAAACTAAACTGATGGCGCTAATGACCTCGCTGAACCCTTTCTGCTCCTCCGGTTTGTCGATGGTGATTTGCACgg                    >  2:125133/2‑251 (MQ=255)                                                                                                                                                                                                             ggTACCTGCCAATGCGGAATATTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCTCTTTTCGCGCCTGCCAGAGGGAGTGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGTTGGTGATGTGCACGGGATg                <  2:489475/251‑1 (MQ=255)                                                                                                                                                                                                              gtgccTGCCAATGCGCAATAGTCCAGCCTGACAGTAATGCATGACTTGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGTGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACGGATGCCGCTACTGACCTGGCTGACCCGATGCAGCTCCGCCGGGTTGCCTATTGTGATTTCCACGGGAt                 <  2:262512/246‑1 (MQ=255)                                                                                                                                                                                                                 ccTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGGCTTTTTCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTCATTACCTGGCTGACCCGTTGCAGCCCCTCCGGTTTGCCGATGGTGATTTGCACGGGGTGCCgg            >  1:434556/1‑250 (MQ=255)                                                                                                                                                                                                                       aaTGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTg      <  1:391586/251‑1 (MQ=255)                                                                                                                                                                                                                           cGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCTACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCCCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTc  <  1:212231/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          GCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTC  >  NC_000913/465913‑466381

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