Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,388,041	T→G	T50P (ACG→CCG)	aaeB ←	aromatic carboxylic acid efflux pump subunit AaeB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,388,041	0	T	G	100.0%	143.9 / NA	47	T50P (ACG→CCG)	aaeB	aromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/25);  total (22/25)

ATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGC‑CATTTTA‑CGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCA  >  NC_000913/3387800‑3388287                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         aTCAGCGCGGTATAACAGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGCTACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAAGAAATGTGCCGATGATGCGCAAAAAGCCACGCTAGCGAATAGCGCCAGAATCCGGTTCACATCCCGCAGCAAAGGCCGGACCggcggc                                                                                                                                                                                                                                                 >  2:9736/1‑251 (MQ=255)       gcgGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATACAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCCAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGcaccact                                                                                                                                                                                                                                           >  1:184746/1‑246 (MQ=255)        cgGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATc                                                                                                                                                                                                                                          >  1:236266/1‑251 (MQ=255)                     cAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTa                                                                                                                                                                                                                             <  1:307762/251‑1 (MQ=255)                              cGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAAACGGCCTAGATACGGCACACCAGAATCATCAATAATGGAGCGCGAATCCTCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGCTAGCGAATAGCGCCCGCATAAGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTACGGACCGCCCAGc                                                                                                                                                                                                                    >  2:364140/1‑251 (MQ=255)                                         tttCTATTCGTACCAGCGAGGAGATCCAGGTACAAAAACCGGCCCAGATACAGCACTCCAGAAGCATCAATAATGGTTCGCGGATCATCGCAATGAGGATCACCATTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCGACGATAGCGAATAGCTCCAGAATCCTGTTCACCTACCTCAGGAAATGGCGGACTTGCGGCAACAATCGCCGCTGTCAGTACCGCCCATCGTGGCGTTTc                                                                                                                                                                                                          <  2:211515/250‑1 (MQ=255)                                                    aCCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAGTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTg                                                                                                                                                                                              <  1:220501/251‑1 (MQ=255)                                                     ccAGCAAGAATATCCAGGTACAAAAACCGGCCCAGATACAACACAACAGAATCATCAATAATGGTACCCAGCTCATCGCAATGCACATCACCAATCCGGCAATACAGCCAATAAATGTGACGAACATGCGCAAAACACCACGATAGCGAATAGAGACAGCATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCACAAAAAGCCGCAGTCAGTACCGCCCAGCGCGGCCTTTCCAGCTGGAAg                                                                                                                                                                                                >  2:245449/1‑248 (MQ=255)                                                                           aaaaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCGATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTa                                                                                                                                                                       <  1:283331/251‑1 (MQ=255)                                                                             aaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCACTGATGACCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACa                                                                                                                                                                     <  1:85596/251‑1 (MQ=255)                                                                              aCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAACGCCAACAAACAGCGCCAGTACa                                                                                                                                                                     >  2:12759/1‑250 (MQ=255)                                                                                          aCAGCACACCAGAAACATCAATAATGGGGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCa                                                                                                                                                        <  2:341673/251‑1 (MQ=255)                                                                                           cAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCGATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAATGCCGGACCGGCGTCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAg                                                                                                                                                       <  1:18504/251‑1 (MQ=255)                                                                                           cAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGcc                                                                                                                                                         >  1:212880/1‑249 (MQ=255)                                                                                             gCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAAGCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCGGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGtt                                                                                                                                                     <  1:419909/251‑1 (MQ=255)                                                                                                 cccAGAATCATCAATAATGGTGCGCGGATCAGCGCAATGATGATCACCAGTCCGGCAAGACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCTGTTTTAc                                                                                                                                                 <  2:123722/250‑1 (MQ=255)                                                                                                    aGAATCAACAAGAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGAGGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTAGCGTTGCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGc                                                                                                                                              <  2:173002/251‑1 (MQ=255)                                                                                                         catcaATAAGGGTGCGCGGATCAGAGCAATGATGAAAACCAGGCCGGAAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGGTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGc                                                                                                                                         <  2:90600/251‑1 (MQ=255)                                                                                                          atcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCg                                                                                                                                        <  1:251107/251‑1 (MQ=255)                                                                                                             aataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTTGCCCGTTTTACCGCAAAGCGAat                                                                                                                                     >  1:297/1‑251 (MQ=255)                                                                                                                   ggTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGtt                                                                                                                                <  1:12759/250‑1 (MQ=255)                                                                                                                      gcgcGGATCATCGCAATGATGACCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCCGCTGGAAGGGAAAGCCaacaaaccg                                                                                                                                                                               >  1:199637/1‑198 (MQ=255)                                                                                                                       cgcgGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATAAGTTgggt                                                                                                                            >  1:34966/1‑248 (MQ=255)                                                                                                                        gcgGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTAGCCAGTTTTAACCCAAAGCGAATAAGTTGGtt                                                                                                                            >  2:389127/1‑249 (MQ=255)                                                                                                                        gcgGATCATCGCAATGATGATCACCAGTCAGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGACTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCAGCTGACAGTACCGCCCAGCGTGGCGTTTCCAGCTGGACGTGAAAGCCAACAAACAGCGCCAGTACAATGGAGGTCCGCAGTTTTaaccgccagcaaaattgtttggt                                                                                                                            >  2:396400/1‑227 (MQ=255)                                                                                                                            atcatcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCaa                                                                                                                       >  1:231868/1‑250 (MQ=255)                                                                                                                              catcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGGTAGCAATgg                                                                                                                    >  1:275669/1‑251 (MQ=255)                                                                                                                                    aaTGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGaaaa                                                                                                              <  1:59062/251‑1 (MQ=255)                                                                                                                                                        aaaTACAGCCAATAAATGTGCCGATGATGCTCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTGCTCCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGGGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATTGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCaa                                                                                          <  2:385044/250‑1 (MQ=255)                                                                                                                                                          aTACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCCAAGCGAATATGTTTGTTAGCAATGGAGAAAATACCCATCGTGATTAACCaaa                                                                                         >  1:338990/1‑250 (MQ=255)                                                                                                                                                                 cAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCa                                                                                 <  1:252705/251‑1 (MQ=255)                                                                                                                                                                  aaTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAg                                                                                <  1:389127/251‑1 (MQ=255)                                                                                                                                                                     aaaTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGc                                                                              <  1:22947/250‑1 (MQ=255)                                                                                                                                                                      aaTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGGGAAAATACCCCACGGGGTTAACCAAAATCACACAGGca                                                                             >  1:192117/1‑249 (MQ=255)                                                                                                                                                                           gCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATggg                                                                        <  1:93767/250‑1 (MQ=255)                                                                                                                                                                             cGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGGTAGCAAAGGGGAAAAAACCCAACGCGAGTAACCAAACTCACGCAGGGGATggggcc                                                                     >  2:217188/1‑248 (MQ=255)                                                                                                                                                                                           aagcccagATACCGAATAGCGCCAGAATGCGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTGTACCGCAAAGCAAATATGTTGGTTAGCAATGGAGAAAATACCCAACGTGATTAACCAAACTCACGCAGGCGATGGGC‑CATTTTA‑CGGAAGAAc                                                     <  2:340842/244‑1 (MQ=255)                                                                                                                                                                                                  gATAGCGAATAGCGCCAGAATACAGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCACGGTTTAACACAAAACGAATATGTTGGTTTACAATGGAGAGAATACACATCGTGGATAACCCAAATAACGCAGGCGGTGGGCGCATTT‑A‑CGGAAGAACGAATCCt                                              >  2:113451/1‑251 (MQ=255)                                                                                                                                                                                                         aaTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGC‑CATTTTA‑CGGAAGAACGAATCCTGGCTTtc                                       <  1:175203/251‑1 (MQ=255)                                                                                                                                                                                                         aaTAGCGCCAGAATACGGTTCACCCCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGC‑CATTTTA‑CGGAAGAACGAATCCTGGCTTtc                                       <  2:245024/251‑1 (MQ=255)                                                                                                                                                                                                               gcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGC‑CATTTTA‑CGGAAGAAAGAATCCTGGCTTTCGTcgcg                                 >  1:291436/1‑251 (MQ=255)                                                                                                                                                                                                                        cGGTTCCCCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGGACCGCCCAGCGTGGCGTTTCCAGCGGGAAGTGAAAGCCAACAAACAGCGCCAGGACAATGGCGGTCGCCAGTTGTACCGCAAAGCGAATATGGTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAGCTCACGCAGGCGATGGGC‑CATTTTA‑CGGAAGAACGAATCCGGGCTTTCGTCGCGATCTTGttt                        <  2:164559/251‑1 (MQ=255)                                                                                                                                                                                                                        cGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGgaga                                                                                                                 <  1:288846/164‑1 (MQ=255)                                                                                                                                                                                                                            tcgccgcccGCAGCACAGGCCGGACCGTCGGCAACAATCGCCGCTGTCAGTACCGCCCCGCGTGGCGTGTCCAGCTGGAAGTGAAAGCCATCAAACAGCGCCAGGACAATGGCGGTCGCCAGTTTTACCGCAAAGCTAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTTATTAACCAAACTCGCGCAGGCGATGGTC‑CATTTTA‑CGGAAGAACTAATCCTGGCTTTCGTCGCGATCTTGTTTGcc                     <  1:17807/244‑1 (MQ=255)                                                                                                                                                                                                                                 tCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGCTGGGC‑CATTTTA‑AGGAAGAACGAATCCAGGCTTTCGTCGCGATCTTGTTTGCCAGTGa                >  1:45822/1‑250 (MQ=255)                                                                                                                                                                                                                                     gcagcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGCACCGCCCAGCGGGGCGTTTCCAGCTGGCAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCCACAAGAATCTGGTGGTTAGGACTGGTGAAAATACCCATCGTAGTTAACCAAAATCACGCAGGCGATGGGC‑CCTTTTA‑AGGCAGAAAGAATCCCAGCGTTCGTCGCCATCTTGATTGGCAGGGaccacc           >  2:32764/1‑251 (MQ=255)                                                                                                                                                                                                                                              gCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGACTCCAGTTTTACCCCAAAGCCAAATTGTTTGGTAGCAATGGAGAAAACACCCCTTCCGATTAACCAAACTCACGCAAGCGATGGGC‑CATTTTAACGGA‑GGACGCATCCTGGCCTTCGTCCCGCTCCTGttttcccacgaccccccccgtagcc  >  2:352531/1‑229 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         ATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGC‑CATTTTA‑CGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCA  >  NC_000913/3387800‑3388287

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