Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,888,227	T→C	E29G (GAA→GGA)	cysH ←	phosphoadenosine phosphosulfate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,888,227	0	T	C	100.0%	179.9 / NA	52	E29G (GAA→GGA)	cysH	phosphoadenosine phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/29);  total (23/29)

TCGGTAGCACGGTACACTTTCAGGTTGAGCTTGAGTTTGTCCGTTAACTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTC  >  NC_000913/2887990‑2888464                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               ccGGTAGCACGGTACACTGTCAGGTTGAGCTTGAGTTTGTCCGTTAACTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGTAAGAAAGTACATATTCACCGTGCAGATGAGCCAGCGCCCAGGCTACGCGGCCTTCACCGTCCAGTGTTCCCAGTTCGGCGtt                                                                                                                                                                                                                                   <  2:204183/250‑1 (MQ=255)        cACGGTACACTTTCAGGTTGAGCTTGAGTTTGTCCGTTAACTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTcc                                                                                                                                                                                                                            <  1:211542/251‑1 (MQ=255)                    tCAGGTTGAGCTTGAGTTTGTCCGTTAACTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCcgccagcg                                                                                                                                                                                                                     <  1:16290/246‑1 (MQ=255)                    tCAGGTTGAGCTTGAGTTTGTCCGTTAACTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCcgccagcg                                                                                                                                                                                                                     >  2:16290/1‑246 (MQ=255)                      aGGTTGAGCTTGAGTTTGTCCGTTAACGCGGCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTAGCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAAt                                                                                                                                                                                                              <  2:337254/251‑1 (MQ=255)                            aGCTTGATGTTGTCCGTTAACGCGTTAATAAAGAGGTAAGTTTCCGGGAACAAGTAACTCGTATCGGTGAGGATCACCGGAAGATCCGGGCTAAGTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCTTGCAGATTATCCAGCGCCCGGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATc                                                                                                                                                                                                        <  1:413752/251‑1 (MQ=255)                             gCTTGAGTTTGTCCGTTAACTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt                                                                                                                                                                                                       <  2:110409/251‑1 (MQ=255)                                       gtacgtgaaCTCGTCAATAAAGCAGTTGGTTCCCGGGAACAAGAAACACGTAGCGGTGAGGATTACCGGAATATCCGGGCGTATTTGATTCACCTGATGCAGGCTCACCGCCGCTTGAATGCCAAAGCTGGAAGAACGTACATATTCACCGGGCACATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGGTAGATGCCGTCAGCGCCAGAATGCGATCTTCCTTCGGcg                                                                                                                                                                                             <  2:353310/244‑2 (MQ=255)                                          cGTTAACTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAATTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAgtt                                                                                                                                                                                          <  1:221546/251‑1 (MQ=255)                                           gTTAACTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAgttc                                                                                                                                                                                         <  2:305599/251‑1 (MQ=255)                                              aaCTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAgttcgtt                                                                                                                                                                                      <  1:96428/251‑1 (MQ=255)                                               aCTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAgttcgttc                                                                                                                                                                                     <  1:248208/251‑1 (MQ=255)                                               aCTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAgttcgttc                                                                                                                                                                                     >  1:378256/1‑251 (MQ=255)                                                       aTAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGtt                                                                                                                                                                             >  1:146798/1‑251 (MQ=255)                                                         aaaGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGttt                                                                                                                                                                            <  1:169678/250‑1 (MQ=255)                                                                 aGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCATTTCGTTCAGGGCGTTTAGATCGAGt                                                                                                                                                                   <  2:36017/251‑1 (MQ=255)                                                                   gTTTCCGGGAACAAGTAACCCGTATTGGGGCGGCCCCCCGGCCTATCCGGGCGCATTTGATTCACCAGCTGCACGCTAACAGCCGCCTGAATGCAACAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCAAGCCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTCGTTTCCGCCAGCCCCAGAATGCGATCTACCTACGGCCGTTCGTTCAGAGCGGTTCGCTCGAGttt                                                                                                                                                                 >  2:412098/1‑251 (MQ=255)                                                                   gTTTCCGGGAACAAGTAACCCGTAGCGGTGAGGATCACCTGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCTCCTTAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTGTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGttt                                                                                                                                                                 <  1:266795/251‑1 (MQ=255)                                                                     ttCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTg                                                                                                                                                                <  1:135811/250‑1 (MQ=255)                                                                       ccGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGAc                                                                                                                                                             <  1:11604/251‑1 (MQ=255)                                                                             aCAAGTAACCCGTATCGGGGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGACTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGCATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAgct                                                                                                                                                       >  2:78024/1‑249 (MQ=255)                                                                                  tAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTc                                                                                                                                                  >  1:41283/1‑251 (MQ=255)                                                                                          aTCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCCGAATGCGATCTCCCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGCTTCCTCACTGTttg                                                                                                                                          >  2:15059/1‑251 (MQ=255)                                                                                            cGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgcc                                                                                                                                        >  2:42913/1‑251 (MQ=255)                                                                                                    aTCACCGGACTCTCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCCGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGATTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTg                                                                                                                                 >  2:311310/1‑250 (MQ=255)                                                                                                      gtgcggagtagccGGGCGAATTGGATTCACCAGAAGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAATTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTTCCTTCGGCAGATCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATg                                                                                                                              <  1:213543/240‑1 (MQ=255)                                                                                                      cACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGAt                                                                                                                               <  1:3116/250‑1 (MQ=255)                                                                                                                       gAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCt                                                                                                             >  1:210283/1‑251 (MQ=255)                                                                                                                            tGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATc                                                                                                         <  2:350895/250‑1 (MQ=255)                                                                                                                            tGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGGTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATc                                                                                                         <  2:47290/250‑1 (MQ=255)                                                                                                                                                 aCCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGAt                                                                                  <  1:146850/251‑1 (MQ=255)                                                                                                                                                        ccTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAAt                                                                           >  2:90010/1‑251 (MQ=255)                                                                                                                                                             aTGCCAAAGCTGGAAGAAAGTGCATATGCGCCGGGCAGAGTATCCAGCGGCCAGGCTACGCGGCCTGCAGCGTCCAGTTTTCCCAGTGCGGCGTTAGTTTCCGCCAGCGCCAGAAGGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGGCATAGATTCCTGACGGTTTTCCTTGCCTGACGCGACGCGCACGCGTGCTTATCAGGCACT‑CAAGACCGGGCTGATGGTTAATCCcac                                                                      <  2:429285/251‑1 (MQ=255)                                                                                                                                                                    aGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGCTGGTTAATCCCACAAATccc                                                               >  2:7582/1‑249 (MQ=255)                                                                                                                                                                      cTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCAcgc                                                             >  1:122624/1‑251 (MQ=255)                                                                                                                                                                         gaagaaAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGcc                                                          >  1:209759/1‑251 (MQ=255)                                                                                                                                                                            gaaAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCgg                                                        >  1:354031/1‑250 (MQ=255)                                                                                                                                                                             aaaGTACATATTCACCGGGCAGATTATCCAGCGCACAGGCTACGCGGCCTTCAGCGTCCAGCTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTAGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCCCACTGTTTGCCTTGCCTGATGCGCCGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGAt                                                      >  1:237173/1‑251 (MQ=255)                                                                                                                                                                              aaGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAAGCACGCGCCGGATc                                                     <  2:311342/251‑1 (MQ=255)                                                                                                                                                                               aGTACATATTAACCGGGCAGATTATCCAGCGCCCAGGCTACGAGGCATTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGCTTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCg                                                    >  1:66525/1‑251 (MQ=255)                                                                                                                                                                                gTACATATTCACCGGGCAGATTATACAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCg                                                    >  1:347476/1‑250 (MQ=255)                                                                                                                                                                                        tCACCGGGCAGATTATCCAGCGCCCAGGCGACGCGGCCTTCAGCGTCCAGTTTGCCCGGTTCGGCTTTAGTTTCCGCCAGCGCCAGAATGCGATCGACCGTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCTCTGGTTGCCTTGCCGGATGCGACGCTCACGCGTCCTTATCAGGC‑CTTCAAGACCGCGCTAATTGTTAATCCCACAAATCACGCGCCGGAGCGAGCACCgg                                            <  1:412098/250‑1 (MQ=255)                                                                                                                                                                                            cGGGCAGATTATCCAGCGACCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCCGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGac                                       >  1:346609/1‑250 (MQ=255)                                                                                                                                                                                              ggCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAatg                                     >  1:252017/1‑251 (MQ=255)                                                                                                                                                                                                cAGATTATCCAGCGCCCAGGCGACGCGGCCAGCACCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTGCCTTCGGCAGTTCGTTCAGGGCGTTTATATCGAGTTTGGACATAGATTCCTCACTGGTGGCCTTGCCTGATGCGACGCTCACGCGTTCGTATCAGGC‑CTACAAGACCTGGCTGAGGGTTAATCCCACAAATCACGCGCCGGATTGAGCTCCGTGCGAatgat                                   <  2:76424/251‑1 (MQ=255)                                                                                                                                                                                                   aTTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGcc                                >  1:350200/1‑251 (MQ=255)                                                                                                                                                                                                    ttATCCAGCGCACAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGccc                               >  1:332051/1‑251 (MQ=255)                                                                                                                                                                                                        ccAGCGCCCAGGCTGCGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCcgcac                           <  2:367354/251‑1 (MQ=255)                                                                                                                                                                                                               ccAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGt                    <  1:273325/251‑1 (MQ=255)                                                                                                                                                                                                               ccAGGCTACGCGGCCTTCAGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGACCACCGGGCGAATGATGCCCGCACGCACCgg                    >  2:377299/1‑250 (MQ=255)                                                                                                                                                                                                                                 aGCGTCCAGTTTTCCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTc  <  1:42913/251‑1 (MQ=255)                                                                                                                                                                                                                                 aGCGTCCAGTTTTCCCACTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCCCTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGGCCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCGCCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTc  <  2:167266/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               TCGGTAGCACGGTACACTTTCAGGTTGAGCTTGAGTTTGTCCGTTAACTCGTCAATAAAGCGGTAGGTTTCCGGGAACAAGTAACCCGTATCGGTGAGGATCACCGGAATATCCGGGCGAATTTGATTCACCAGATGCAGGCTCACCGCCGCCTGAATGCCAAAGCTGGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGC‑CTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTC  >  NC_000913/2887990‑2888464

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