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%	224.1 / NA	72	T50P (ACG→CCG)	aaeB	aromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (33/39);  total (33/39)

AATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGC  >  NC_000913/3387796‑3388286                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          actGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcgg                                                                                                                                                                                                                                                      >  2:147122/3‑251 (MQ=255) aaTGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcgg                                                                                                                                                                                                                                                      <  2:117154/251‑1 (MQ=255)  aTGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCTCAGCAAAGGCCGGACcggc                                                                                                                                                                                                                                                     <  1:182143/251‑1 (MQ=255)     aTCAGAGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGAAGCAAAGGCCGGCCCggcggc                                                                                                                                                                                                                                                  >  2:267261/1‑251 (MQ=255)      tCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTAGTACCAGCAAGGATATCCAGGTACAAACACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGCTCACCCGTCCGGCAATACAGCCACTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAACGGCCGGACCggcggc                                                                                                                                                                                                                                                  >  1:275855/1‑250 (MQ=255)                   aCCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATCTACAGGGACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATcgccgcgg                                                                                                                                                                                                                                    >  1:186064/1‑249 (MQ=255)                    ccGGCCAGCCCCCACGCATTCGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCATAATCATCAATAATGGTGCGCGGATCATCGCAAGGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAAGACGGTTCACCTCCCGCAGCAAATGCCGGACCGGCGGCAACAAGCGCCGCTGt                                                                                                                                                                                                                                   <  2:460414/251‑1 (MQ=255)                     cGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTc                                                                                                                                                                                                                                  <  2:179933/251‑1 (MQ=255)                               ccACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGccc                                                                                                                                                                                                                        >  2:349423/1‑251 (MQ=255)                                 aCGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGAGCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCCGAATACGGTTCAACTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTTAGTAACGcca                                                                                                                                                                                                                        >  1:299703/1‑248 (MQ=255)                                    cATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATCCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCTCGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGGCAGTACCGCCCAGCGt                                                                                                                                                                                                                   <  1:350426/251‑1 (MQ=255)                                       aCGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGc                                                                                                                                                                                                                <  1:284201/251‑1 (MQ=255)                                       aCGAGTTTTCTAGTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCTCTGTCAGTACCGCCCAGCGTgg                                                                                                                                                                                                                 <  1:221471/250‑1 (MQ=255)                                           gTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGttt                                                                                                                                                                                                            <  2:140290/251‑1 (MQ=255)                                           gTTTTCTATTCGTACCAGCGAGGATATCCAGGAACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGttt                                                                                                                                                                                                            >  1:14252/1‑251 (MQ=255)                                             tttCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTc                                                                                                                                                                                                           <  1:285120/250‑1 (MQ=255)                                                     cGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGCGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGa                                                                                                                                                                                                   >  1:84118/1‑250 (MQ=255)                                                      gTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGa                                                                                                                                                                                                   >  2:178662/1‑249 (MQ=255)                                                              gAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGcc                                                                                                                                                                                         <  2:160290/251‑1 (MQ=255)                                                                      ccAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCaacaaaca                                                                                                                                                                                 <  1:4622/251‑1 (MQ=255)                                                                                        ccAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCa                                                                                                                                                                                                         >  1:70979/1‑209 (MQ=255)                                                                                        ccAGATACAGCACACCAGAATCATCAATAATGGGGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCGCCTCCCGCAGCAAAGGCCTGGCCGGCGGCAACAAGCGCCGCTGTCAGTACCGCCCAGCGTGGCTTTTCCa                                                                                                                                                                                                         <  2:70979/209‑1 (MQ=255)                                                                                          aGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGtt                                                                                                                                                                                                             <  2:131160/203‑1 (MQ=255)                                                                                          aGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGtt                                                                                                                                                                                                             >  1:131160/1‑203 (MQ=255)                                                                                          aGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTc                                                                                                                                                             <  1:155313/251‑1 (MQ=255)                                                                                          aGATACAGCACACCAGAATCATCAATAAGGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTc                                                                                                                                                             <  2:369951/251‑1 (MQ=255)                                                                                           aatccAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCg                                                                                                                                                            <  2:144606/247‑1 (MQ=255)                                                                                            aTACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgcg                                                                                                                                                                              <  1:232211/232‑1 (MQ=255)                                                                                            aTACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgcg                                                                                                                                                                              >  2:232211/1‑232 (MQ=255)                                                                                                  cacaCCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGttt                                                                                                                                                     <  1:37501/251‑1 (MQ=255)                                                                                                     acCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAc                                                                                                                                                  >  1:128760/1‑251 (MQ=255)                                                                                                        aGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGc                                                                                                                                               >  1:278550/1‑251 (MQ=255)                                                                                                          aaTCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAATCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCaa                                                                                                                                             <  1:401549/251‑1 (MQ=255)                                                                                                                  ataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGACACGATAGCGAATAGCGCCAGAATCCGGTTCACCTCCCGCAGCAAATGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGCAGTGAAAGCCCACAAACAGCGCCAGTACAATGGCGGGCGCCCGTTTTACCGCCAACCGAat                                                                                                                                      >  1:30503/1‑250 (MQ=255)                                                                                                                     atGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑t                                                                                                                                 <  1:178662/251‑1 (MQ=255)                                                                                                                      tggggcgcGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑t                                                                                                                                 <  2:11785/246‑1 (MQ=255)                                                                                                                      tGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAGTTTTACCGCAAAGCGAATATG‑t                                                                                                                                 >  1:11785/1‑250 (MQ=255)                                                                                                                           cgcgGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGAGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGttt                                                                                                                           >  2:153639/1‑250 (MQ=255)                                                                                                                             cgGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGc                                                                                                                         <  1:216164/251‑1 (MQ=255)                                                                                                                               gATCATCGCAATGATGATACCCAGTCCGGCAATAAAGCCAATAAATGTGCAGAGGAGGCGCAAAAAGCCACGAAAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCa‑a                                                                                                                      <  2:392955/251‑1 (MQ=255)                                                                                                                                     cGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGgag                                                                                                                 <  1:147122/250‑1 (MQ=255)                                                                                                                                          tgatgaTCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGCTGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGACGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGATGGAAGTGAAAGCCAAAAAACAGCGCCAGTACAAAGGCGGTCGCCAGGTTTACCACCACGCGAATATG‑TTGGGTAGAG‑ATGGGGAAAAt                                                                                                            >  1:29148/1‑250 (MQ=255)                                                                                                                                           gatgatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAc                                                                                                          >  1:114383/1‑251 (MQ=255)                                                                                                                                           gatgacCACCAGTCCGGCAATACAGCCAAGAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATa                                                                                                           <  2:434127/250‑1 (MQ=255)                                                                                                                                            atgatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAc‑c                                                                                                        <  1:261342/251‑1 (MQ=255)                                                                                                                                               atCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCAt                                                                                                     <  1:152317/251‑1 (MQ=255)                                                                                                                                               atCACCAGTCCGGCAATACAGCCAATAAAAGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCAAGCATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCCGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGGACAATGGCGGACGACAACTTTTCCGCAAAAGGAAAATG‑TTGGTTAGCA‑Ataggggaaaaaacacc                                                                                                      >  2:193696/1‑236 (MQ=255)                                                                                                                                                   ccAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATCC‑CCATCGTg                                                                                                 <  2:186064/251‑1 (MQ=255)                                                                                                                                                           gCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCAGGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGGTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCa                                                                                         <  2:10091/251‑1 (MQ=255)                                                                                                                                                           gCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCa                                                                                         <  1:192223/251‑1 (MQ=255)                                                                                                                                                               tACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAA‑CCATCGTGATTAACCaaaa                                                                                      >  1:195917/1‑249 (MQ=255)                                                                                                                                                                    ccAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGc                                                                                <  2:343924/251‑1 (MQ=255)                                                                                                                                                                      aaTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTAGCCCGTTTTACCGCAAAGCGAATATG‑GTGGTTAGCA‑ATGGGGGAAATAA‑CCAACGTGATTAACCAAAATCAcac                                                                                >  1:9104/1‑247 (MQ=255)                                                                                                                                                                      aaTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCAGCTGTCAGTACCGACCAGAGTGGCGTTTCCAGATGGAAGTGAAAGCCAACAAACCGCGCACGTACAATGGCGGGCGCCATGTTTAACGCAAAAGGAAAATG‑TTGGTTAACA‑ATGGAGAAAATACACCA‑CGTGTTTAACCaacatcacaca                                                                               >  2:297971/1‑241 (MQ=255)                                                                                                                                                                       aTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAgg                                                                             <  2:257972/251‑1 (MQ=255)                                                                                                                                                                        tAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGGATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTTACCAAACTCACGCAggg                                                                            >  1:376157/1‑250 (MQ=255)                                                                                                                                                                         taaTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAAGACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCTAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACTCAGGc                                                                            <  2:7250/249‑1 (MQ=255)                                                                                                                                                                                                  ccACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAACGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTTCGGTAGAAC‑Ga                                                  >  1:248524/1‑250 (MQ=255)                                                                                                                                                                                                    aCGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCCCAAAGCGCATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAAACAAACTCCCGCAAGCGAGGGGCCATTTTACGGAAGAAC‑GAATc                                               >  1:202626/1‑251 (MQ=255)                                                                                                                                                                                                         aGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGGc                                          <  2:84118/251‑1 (MQ=255)                                                                                                                                                                                                           cGAATAGCTCCATAATTCGGTGCACCTCCAGCAGCAAAGGCCGGACCGGCGGCAACAATCGGCGCTGTCAGTACCGCCCAGCGTGGCGGTTCCAGCTTGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCTGACTCCAGTTTTACCGCAAAGGGAATATGAGTGGT‑AGCA‑ATTTAGAAATTAT‑CCGTCGTGATTAAGCAAACTCACTCAGGCGAGGGTCCAtttgtcgcaagatcggagcctggctt                                         <  2:206365/251‑24 (MQ=255)                                                                                                                                                                                                           cGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGGGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGGCtt                                        >  1:163016/1‑251 (MQ=255)                                                                                                                                                                                                           cGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGCGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTAACGCAAAACGAATATG‑TTTGTTAGGA‑ATTGAGAAAATAC‑CCATCGGGATTAAACAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACAGAA‑CCTGGCtt                                        >  1:361903/1‑251 (MQ=255)                                                                                                                                                                                                                aGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGGCTTtcgtc                                   <  2:52725/251‑1 (MQ=255)                                                                                                                                                                                                                     cAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGGCTTTCGTCGCGAt                              <  1:143092/251‑1 (MQ=255)                                                                                                                                                                                                                     cAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑GTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGGCTTTCGTCGCGAt                              >  1:88377/1‑251 (MQ=255)                                                                                                                                                                                                                      aGAATACGGTTCACCTCCCGCATCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGTCTTtcg                                     <  1:362849/243‑1 (MQ=255)                                                                                                                                                                                                                      aGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGGTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAGGTTTACCGCAAAAAGAATATG‑GTGGGTAACATAGGGGGGAAA‑AC‑CCATCGTGATTAAACACACTAAAGCAGGCGAGGGGGGAGTTTACTGGAGAAC‑Gaagtacggattttgaggttggg                             >  2:362849/1‑231 (MQ=255)                                                                                                                                                                                                                             ggTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGGCTTTCGTCGCGATCTTGTTTg                      >  1:6679/1‑251 (MQ=255)                                                                                                                                                                                                                                      cccGCAGCAAAGGCCGGTCCTGCTGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATTC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGa               <  1:426516/249‑1 (MQ=255)                                                                                                                                                                                                                                      cccGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGGATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGCGATTAACCAAACTCACGCCGGCGATGGGCCATTTTACGGAAGAAC‑GAAACCTGGCTTTCGTCGCGATCTTGTTTGCCAGTgcac             >  1:74237/1‑248 (MQ=255)                                                                                                                                                                                                                                                  gCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCCAAGCGAATATG‑TTGGTTAGCC‑ATGGGGAAAATAC‑CCATCGTGATTAACCAAACTCCCGCCGGCGATGGGGCATTTTACCGGAGAAA‑GAATCCTGGCTTTCGTCGCCATCCTGTTTGCCATTGACCCCCACTGTAGc  >  2:320425/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          AATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATG‑TTGGTTAGCA‑ATGGAGAAAATAC‑CCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAC‑GAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGC  >  NC_000913/3387796‑3388286

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