Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	162,339	A→G	T79A (ACC→GCC)	hrpB →	RNA‑dependent NTPase HrpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	162,339	0	A	G	100.0%	201.9 / NA	66	T79A (ACC→GCC)	hrpB	RNA‑dependent NTPase HrpB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/31);  total (35/31)

TCTGGAGCGTTTTGTGTCGTCGTTGCCCGTTGCTGCCGTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGG  >  NC_000913/162092‑162582                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    tctGGAGCGTTTTGTGTCGTCGTTGCCCGTTGCTGCCGTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCg                                                                                                                                                                                                                                                  >  2:21787/1‑251 (MQ=255)       gCGTTTTGTGTCGTCGTTGCCCGTTGCTGCCGTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCt                                                                                                                                                                                                                                            >  2:82286/1‑251 (MQ=255)                              ttGCTGCCGTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATTCGTGCGCAAAACTGc                                                                                                                                                                                                                     >  2:366789/1‑251 (MQ=255)                                   gcCGTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATGAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCGACCGGATGCGTGCGCAAAACTGCGTCg                                                                                                                                                                                                                 <  2:191151/250‑1 (MQ=255)                                      gTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCTCAGCTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCCACCTGACTGCCGCTGCAACTGCTGGCGCCTCCCGGCATTAACGGGAACCTTATCCTGCTGCAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCACCGGCTGGCGGCGCTGCTTAACGAAAAGCCACGCGATGCCGTTGTCTACCGGATCCGCGCGCAAACCTGCGTCGGGcc                                                                                                                                                                                                             >  1:60604/1‑251 (MQ=255)                                           aCCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAAt                                                                                                                                                                                                         <  2:20399/250‑1 (MQ=255)                                               gAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACGGCGGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCTGGCCGAATACCCg                                                                                                                                                                                                    <  2:339540/251‑1 (MQ=255)                                                 aTTACTTACCGCTCTCGATTGTGCGCCGCAGGGATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGcc                                                                                                                                                                                                  <  2:119811/251‑1 (MQ=255)                                                      ttacCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCGTTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGGGCGCAAAACTGCGTCGGGCCGAATACCCGCCTggc                                                                                                                                                                                              >  2:238924/1‑249 (MQ=255)                                                                            gcAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTg                                                                                                                                                                        >  1:107594/1‑250 (MQ=255)                                                                             cAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCTTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCTTTGGCTACGGGATGCGTGCGCAAAACTGCGTCGGTCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGAc                                                                                                                                                                      <  1:165504/251‑1 (MQ=255)                                                                                     attaAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATCCCCGCCTGGAAGTGGTTACCGAAGGCGGGCTGACGCGCATg                                                                                                                                                               >  2:364466/1‑250 (MQ=255)                                                                                            gcgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGc                                                                                                                                                        <  1:82286/250‑1 (MQ=255)                                                                                             cgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGt                                                                                                                                                      <  1:158232/251‑1 (MQ=255)                                                                                             cgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAAAGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGGATACCCGCCTGGCAGTGGTGACCGAAGGCGTGCTGACGCGCATGATCCAGCgg                                                                                                                                                      >  2:104446/1‑250 (MQ=255)                                                                                                cGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGAc                                                                                                                                                   >  1:130591/1‑251 (MQ=255)                                                                                                    cGGGGCCGGGAAATCAAACTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGGGCACAAATGCGTCGGGCCGAATACCCGCCTGGGAGTGGGTAACGAAGGCGTGCTGACGCGCATGAtaaaccgagatccgg                                                                                                                                               >  2:241213/1‑237 (MQ=255)                                                                                                        gCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACt                                                                                                                                           >  1:109246/1‑251 (MQ=255)                                                                                                            ggAAAGCAACCGGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCGGGAGCCGCTTCGTCTGGCGGCGCTTAACGTCGCGCAACGTCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTGGGCTACCGTATGCGTTCGCAAAACGGCGTCGGGCCAAATTGCCGCCTGTAAGTGGTTTCCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGc                                                                                                                                       <  1:211063/251‑1 (MQ=255)                                                                                                                aTCACCCTGGCTGCCGCGGCACCTGCTGGCGCAGCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGTGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGAGGCCGTGGGCTACCGGATGGGTGCGCAAAACGGCGTCGGGCCTAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGgtg                                                                                                                                   <  2:74327/251‑1 (MQ=255)                                                                                                                aTCAACCGGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTACCCTGCTGTGGCGGCGTCGTCGGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAACCAGGCGAGGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGgtg                                                                                                                                   <  2:261946/251‑1 (MQ=255)                                                                                                                 ccaccctggtggCCGCTGCACCTGCTGGCGCATCCCGACATTAACGGGAAAATTATCCGGCTGGGGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATGCCCGCCTAGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGgtgt                                                                                                                                  <  2:204269/240‑1 (MQ=255)                                                                                                                   aaCCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTg                                                                                                                                <  2:109246/251‑1 (MQ=255)                                                                                                                   aaCCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTg                                                                                                                                >  1:210893/1‑251 (MQ=255)                                                                                                                       gggCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACt                                                                                                                            <  1:33402/250‑1 (MQ=255)                                                                                                                            gCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGGTGGACTGGTg                                                                                                                        <  1:105389/250‑1 (MQ=255)                                                                                                                             ccGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGAt                                                                                                                      <  1:95825/251‑1 (MQ=255)                                                                                                                                 tGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCtt                                                                                                                  >  1:185241/1‑251 (MQ=255)                                                                                                                                          tgGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAAttt                                                                                                         >  1:233511/1‑251 (MQ=255)                                                                                                                                               cATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATg                                                                                                     >  1:208693/1‑250 (MQ=255)                                                                                                                                                       cATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGCAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGcc                                                                                            >  2:206244/1‑250 (MQ=255)                                                                                                                                                       cATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGCGCGGTGTTGGACTGGTGATCCTTGATTAATTTCCTGAGCGCAGcc                                                                                            >  2:172978/1‑250 (MQ=255)                                                                                                                                                       cATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCt                                                                                            >  1:38665/1‑251 (MQ=255)                                                                                                                                                           aaCGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCa                                                                                        <  1:78669/251‑1 (MQ=255)                                                                                                                                                                 aaaaTTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCCTGAGCGCAGCTTGCAGGCggg                                                                                  >  1:121602/1‑250 (MQ=255)                                                                                                                                                                   aaTTATCCGGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGAtt                                                                                <  2:233511/251‑1 (MQ=255)                                                                                                                                                                    aTTATCCTGCTGGAGCCGAGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCCTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGGGCGCAAAACTGCGTCGGGCCCAATACCCGCCTGGAAGGGGTTACCGAAGGCGTGCTGACGCGCCTGATCCAGCGGGACCCGGAACTGAGCGGTGTTGGACTGGTGCTCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGAttt                                                                               >  2:63728/1‑251 (MQ=255)                                                                                                                                                                       aTCCTGCGGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGTATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttgg                                                                             <  1:199578/250‑1 (MQ=255)                                                                                                                                                                           tgctgGAACCGCGTCGGCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGGTGCGGGCGCAAAACTGCGGCGGGCCGGATACCCGCCTGGAGGTGGTTACCGAAAGCGTGCTGACGCGCATGATCCAGCGTGACCCGGGACTGGGCGGGGTTGGGCCGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGGAGGCGGATttggcgttg                                                                        >  2:51229/1‑251 (MQ=255)                                                                                                                                                                               ggAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggcgttggcg                                                                     <  1:132723/250‑1 (MQ=255)                                                                                                                                                                               ggAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggcgttggc                                                                      <  1:364466/249‑1 (MQ=255)                                                                                                                                                                               ggAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggcgttggc                                                                      <  1:35098/249‑1 (MQ=255)                                                                                                                                                                                 ggcggggtcgtcTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGAGGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTg                                                                  <  2:347773/245‑1 (MQ=255)                                                                                                                                                                                       gtcgtcTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCtt                                                                                                                  <  1:289663/197‑1 (MQ=255)                                                                                                                                                                                       gtcgtcTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCtt                                                                                                                  >  2:289663/1‑197 (MQ=255)                                                                                                                                                                                            cgggcgtcGCGTAACTTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCGGGCGATGCCGTTGGCGACCGGAGGCGTGCGCAAACCTGCGGCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGAtgt                                                       <  2:357233/249‑1 (MQ=255)                                                                                                                                                                                                        aaCGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCt                                           <  2:253771/251‑1 (MQ=255)                                                                                                                                                                                                           gTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCAAGGCGCTGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCACGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCg                                        >  1:338773/1‑251 (MQ=255)                                                                                                                                                                                                             cgcgCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGTGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAgggc                                            >  2:319016/1‑243 (MQ=255)                                                                                                                                                                                                             cgcgCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGTGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTc                                            <  1:319016/245‑1 (MQ=255)                                                                                                                                                                                                                 cAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGtgatga                                  <  1:327748/251‑1 (MQ=255)                                                                                                                                                                                                                  aaCGGCTGGCGGAGCTGCTTAACGAAACGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGAcg                                                                                                                                                                     >  1:295872/1‑119 (MQ=255)                                                                                                                                                                                                                  aaCGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGAcg                                                                                                                                                                     <  2:295872/119‑1 (MQ=255)                                                                                                                                                                                                                           cGGAGCTGCTTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACTCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAActg                        <  1:206244/251‑1 (MQ=255)                                                                                                                                                                                                                            ggAGCTGCTTAACGAAAAGCCAGGCGAGGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAActgc                       <  2:20894/251‑1 (MQ=255)                                                                                                                                                                                                                               gctgctTAACGAAAAGCCAGGAGATGCCGTTGGCTACCGGCTGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGACTGCGTGATGACCTTAAActgctg                     >  1:369463/1‑250 (MQ=255)                                                                                                                                                                                                                                   ctTAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTAt                >  2:344039/1‑251 (MQ=255)                                                                                                                                                                                                                                     tAACGAAAAGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTAt                <  2:121602/249‑1 (MQ=255)                                                                                                                                                                                                                                          caaaGCCAGGCGCTGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCAGGAACTGAGCGGTGTTGGACTGGGGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCGACAAGGTCTGCGTGCTGACCTTAAACTGCTGATTATGTCggat         >  2:85349/2‑249 (MQ=255)                                                                                                                                                                                                                                          aaaaGCCAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCt         >  2:269382/1‑251 (MQ=255)                                                                                                                                                                                                                                               ccAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCt    >  2:169008/1‑251 (MQ=255)                                                                                                                                                                                                                                                cAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTg   >  2:332340/1‑251 (MQ=255)                                                                                                                                                                                                                                                cAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTg   >  1:226417/1‑251 (MQ=255)                                                                                                                                                                                                                                                cAGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTg   >  2:213978/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaGGCGATGCCGTTGGCTACCGGATGCGAGCGCAAAACTGCGTAGGGCAGAATAAACGCCAGGAAGTGGTTACCGAAGGCGTGCTGACGAGCATGATCCAGAGTGACACGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTAACCTg   >  1:226416/2‑251 (MQ=255)                                                                                                                                                                                                                                                 aGGCGATGCCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGACTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTgg  >  1:60690/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    TCTGGAGCGTTTTGTGTCGTCGTTGCCCGTTGCTGCCGTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGG  >  NC_000913/162092‑162582

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