Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,123,930	A→T	V176E (GTA→GAA)	cytR ←	DNA‑binding transcriptional repressor CytR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,123,930	0	A	T	100.0%	265.2 / NA	70	V176E (GTA→GAA)	cytR	DNA‑binding transcriptional repressor CytR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (30/40);  total (30/40)

GAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTGATTTG  >  NC_000913/4123683‑4124176                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         gAGCGCCATCACATCGCTATGGCAGAAGACATCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATTCCCTTGCTCATATAAATAATTTTCTg                                                                                                                                                                                                                                                       <  1:347764/251‑1 (MQ=255)     gcCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATc                                                                                                                                                                                                                                                   >  1:359863/1‑251 (MQ=255)               cGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGAGCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATgcggcg                                                                                                                                                                                                                                         <  2:462321/251‑1 (MQ=255)                     ggCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTcca                                                                                                                                                                                                                                    >  1:56383/1‑248 (MQ=255)                     ggCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGAGACACCGCGGGAACTCTTCGGGAACGGCTATACCGCCAACCCGGTTATGCCCTTTCTCATATAAATAATTTTTTGCCTCCAATgcggggggccg                                                                                                                                                                                                                                    >  2:86843/1‑243 (MQ=255)                              cagcagTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGCGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAACGGCATCTCTTCGGGACCGGCTATACAGCCAATCCCTTTATGCCCTTTCTCATATAAATAATTTTCTGCATCAAATTcggcgggccaattttcgat                                                                                                                                                                                                                           >  1:473431/1‑237 (MQ=255)                                   gTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAAc                                                                                                                                                                                                                    <  2:345225/251‑1 (MQ=255)                                            tgtgGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGca                                                                                                                                                                                                            <  1:363669/250‑1 (MQ=255)                                            tgtgGAAAATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTTAAGTCGCCACGGGCGATGTATTGCGAATCAACCATAATGCCGCATCGACTCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACCGCGGCATCCCTTCGGGACCGGCTATACAGCCAATCCGTTTGTGGCCTTTCTTATATAAATAATTTTCTGCAACAAATGCGGCGGTGCGAATGTCGAT‑ATgtaatttttgc                                                                                                                                                                                                             >  1:347978/1‑239 (MQ=255)                                              ttgAAGATCAAGCAGCTGCTGCATCGCTTTGCTCCCTGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCGTGCTCATATAAATAATTTTCTGCATCAAATGCTGCGGCCAGATTGTCGAG‑ATGAACTTTAGGcagc                                                                                                                                                                                                          <  2:279077/248‑1 (MQ=255)                                                 aaGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCTATGTATTGCGGATCAGCCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGcagctcca                                                                                                                                                                                                      <  1:322866/251‑1 (MQ=255)                                                     tCAAGCAGCTGCTGCCTAGCTTTGCCTCCGGCTCCGAAGGCGACCTCGCCACGGGCGATGTATTGCGGATCAACCCTAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGCGACCCAGCGGCCTCTCTTCGGGACCGGCTATACAGCCAATCCGTTTCTGACCTTGCTCATATAAATAATTTTCTTCCTCAAATGCGGCGGGCAGATTTTCGgt‑atttaccagaggacacccccagcc                                                                                                                                                                                                  >  2:31944/1‑229 (MQ=255)                                                           agcTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGc                                                                                                                                                                                            <  1:262320/251‑1 (MQ=255)                                                            gctgctGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGTCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCa                                                                                                                                                                                           <  1:80713/251‑1 (MQ=255)                                                               gcggcATCGCTTTGCTTCAGGCTTCGAAGGGGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCaa                                                                                                                                                                                          <  2:23880/246‑1 (MQ=255)                                                                  gcATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCCGATTTGCGAT‑ATTGACTGTAGGCAACTCCAGCTCCGGGGCAAAttcg                                                                                                                                                                                     >  1:378159/1‑251 (MQ=255)                                                                      cGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAAttcgttcg                                                                                                                                                                                 <  2:6225/251‑1 (MQ=255)                                                                      cGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCGTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAAttcgttc                                                                                                                                                                                  <  1:227889/250‑1 (MQ=255)                                                                          ttGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGGCACAGCGGCATCTCTTCGGGACCGGCTATACCGCCAATCCGTTTATGCCCTTGCTCCTATAAATAATTTTCTGCATCAAATGCGGCGGGCAGGTTTTTGAT‑ATGAATCGTATGCAGCCCCCGCCCCCGTGCACATTCGTTTGcccc                                                                                                                                                                             >  2:24504/1‑249 (MQ=255)                                                                               tCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCAcc                                                                                                                                                                         <  2:132325/250‑1 (MQ=255)                                                                               tCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCAc                                                                                                                                                                          <  1:187173/249‑1 (MQ=255)                                                                               tCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAAGTCGTTCGCCATCAc                                                                                                                                                                          <  1:86843/249‑1 (MQ=255)                                                                                     ttCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGATATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑TTTAACCGTAGGCAGCTCCAGATCCGGGGCAAATTTGTTCGcccttcccctccgcg                                                                                                                                                                  >  1:183877/1‑238 (MQ=255)                                                                                                          ggCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCGTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCTGTCAGATTGTCTAT‑ATGAACTGTAGGCGGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCaa                                                                                                                                             <  1:31944/251‑1 (MQ=255)                                                                                                              aTGTATTGCGGATCAACCCTAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCt                                                                                                                                         <  2:336724/251‑1 (MQ=255)                                                                                                                gTATTGCGGAGCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGGCGAT‑ATTAACTGGAGGCAGCTCCAGCTCCGGTTCAAATTCTTTCTCCATCATCATCGGCTGCAGATTTCGCTGTTCCTCATTGCTg                                                                                                                                        <  2:425141/250‑1 (MQ=255)                                                                                                                         aTCAACAATAAGGCCGCCGCGACGCAGCGACGGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCTAATCCTTTTATGCCATTGCTCATATAATTAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACg                                                                                                                              <  1:403972/251‑1 (MQ=255)                                                                                                                              ccATAAGGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACCCAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑c                                                                                                                        <  2:54954/251‑1 (MQ=255)                                                                                                                                aTAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAAGTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCt                                                                                                                      <  2:161476/251‑1 (MQ=255)                                                                                                                                          aGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACCCAGCGGCATCTCTTCGGGACCGGCTATACAGCCCATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAAATGTAGGCAGCTCCAGCTCCggggc                                                                                                                                                                                            >  1:311011/1‑169 (MQ=255)                                                                                                                                          aGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGc                                                                                                                                                                                            <  2:311011/172‑1 (MQ=255)                                                                                                                                           gCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGca                                                                                                           <  2:359863/251‑1 (MQ=255)                                                                                                                                            cGACGCCGCGCCTTAACAGAGCCTAGCAGGCGGTAGGGACACAGCGGCATCTCTTCGGAACCGTCGATACAGCCAAGCCGTTTATGCCCTTGCTCATATAAATAATTTTCTTCATCAAAAGAGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAACTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATTAAACGGCAG‑CCTTGAACCCAGcaa                                                                                                          <  2:462985/251‑1 (MQ=255)                                                                                                                                               cGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCt                                                                                                                                         <  2:178893/218‑1 (MQ=255)                                                                                                                                                 cAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGAATGTCGGT‑ATTAACTGTAGGCAGCTCCCACTCCGGGGGAAATTCGGTCGCCATCACCATCGGCGGCCGATTAAGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAAACCAGCaacaacc                                                                                                     >  1:142353/1‑250 (MQ=255)                                                                                                                                                  aGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACAt                                                                                                    <  1:250038/251‑1 (MQ=255)                                                                                                                                                  aGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCAGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACAt                                                                                                    <  1:361060/251‑1 (MQ=255)                                                                                                                                                     gcCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGcc                                                                                                 <  1:303712/251‑1 (MQ=255)                                                                                                                                                      ccTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCCTCACTTCGGGACCGGCTATACAGCCAACCCGCTTATGCCCTTCCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGACTGTCGATTAT‑AACTGTAGGCCGCTCCAGCTCCCGGGCAAAATCGTTTCCCCTCACCCTCGGCGGCAAATTACGCAGTTTCTCAACGCTGGCATCAAACAGAAGACTTTAAACCCA‑CAACAACATGccc                                                                                                >  1:428771/1‑250 (MQ=255)                                                                                                                                                          aaCATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCATATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATGCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCTTCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCaa                                                                                            <  1:230078/251‑1 (MQ=255)                                                                                                                                                            cATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAAtt                                                                                          <  1:416183/251‑1 (MQ=255)                                                                                                                                                                          cGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTgatgat                                                                            <  1:273938/251‑1 (MQ=255)                                                                                                                                                                           ggTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATc                                                                           <  2:348777/251‑1 (MQ=255)                                                                                                                                                                             gaGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATGTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACGGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTGCCTCAATTCTGGCATCAAACGGCAG‑CCTTGAGCCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCaa                                                                         <  2:278419/250‑1 (MQ=255)                                                                                                                                                                              aGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGGTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCaa                                                                         <  2:377235/250‑1 (MQ=255)                                                                                                                                                                                  acacAGCGGCATCTCTTCGGGACCGGCTATCCCGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGATCCAGCTCCGGGGCAAATTCGTTCGCCATAACCATCTGCGGACGATTACGCTGGTCCTCAAAGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAAAAACATGCCCTCACTTTGCTTGGTGATGGTCAAATCg                                                                     >  2:231166/1‑250 (MQ=255)                                                                                                                                                                                     caGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGGGATGATCAAATCGATa                                                                  >  1:384156/1‑250 (MQ=255)                                                                                                                                                                                       gCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAg                                                               <  2:171414/251‑1 (MQ=255)                                                                                                                                                                                          gCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGGACTGTAGGCAGCTCCAGCTCCGGGGCCAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCTTCAATTCTGGGATCAAACGGCAT‑CCTTGAACCCAGCAACTACATGCCATCAATTTGCTTGGTGATGAGCAAATCGATAAAGGtt                                                            >  1:329507/1‑251 (MQ=255)                                                                                                                                                                                                ctTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGtttttt                                                        <  1:231166/249‑1 (MQ=255)                                                                                                                                                                                                ctTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGtttt                                                          <  2:378159/247‑1 (MQ=255)                                                                                                                                                                                                     ggACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGGTTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTTCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTTCTTGGTGATGATCAAATCGATAAAGGTTTTTTCctgctg                                                 >  1:377052/1‑251 (MQ=255)                                                                                                                                                                                                           gCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATTATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTg                                           >  1:386925/1‑251 (MQ=255)                                                                                                                                                                                                                   gCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCCTGAACCCAGCAACAACATGCCCTCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTTATGTGcac                                   >  2:138448/1‑251 (MQ=255)                                                                                                                                                                                                                      aaTCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGt                                <  1:186761/251‑1 (MQ=255)                                                                                                                                                                                                                       aTCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGCACTGTAGGCAGCCCCCGCTCCGGAGCAAATTCGTTCGCCATCAACATCGGCGGCAGATTAAGCCGTTCCTCAATGCTCGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATACAGGGTTTTTCCTGCTGATTTTGATGCGCACAGTc                               >  1:450745/1‑251 (MQ=255)                                                                                                                                                                                                                        tCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAACCCCAGCTCCGGTGAAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGGGATGATCAAATAGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCg                              >  1:217875/1‑251 (MQ=255)                                                                                                                                                                                                                         ccGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGc                             >  1:135981/1‑251 (MQ=255)                                                                                                                                                                                                                          cGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCCGATTGTCGAT‑ATGAAATGTAGGAAGCTACAGCTCCGGTGCAAATTCGTTCGCCCTCACCAACGGCGGAAGATTACCCTATTCCTCAATTCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTCGTGGTGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGcc                            >  1:175006/1‑251 (MQ=255)                                                                                                                                                                                                                           gTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCCGATTACGCTGTTCCTCAATGCTGGCATCACACGGCAG‑CCTTGAACCCCGCAACAACATGCCATCAATTTGCTTGGTGATGATCCAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGcc                            >  1:65608/1‑250 (MQ=255)                                                                                                                                                                                                                            tttATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGACCTGTTGGCAGCCCCAGCTCCGGGGGAAATTCGTTCGCCATCACCATCGGCGGAAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑GCTTGAACCCAGCAACAACATGCCATCCATTTGCTTGGGGATGATCAAAACGCTAAAGGTTTTTTCCTGCTGGTTTTGATGCGCACAGTCGCCa                           >  1:478738/1‑250 (MQ=255)                                                                                                                                                                                                                                 gCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTTTCGCT‑ATGAACTGTAGGCAGCTCCCACTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCCGATTACGCTGTTCCTCCATGCTTGCATCAAACGGCAG‑CCTTGTACCCAGCAACAAAATGCCCTCAATTTGCTTTGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGGTTTTTATGCCCCCAGTCCCCAATCAGc                     >  1:92533/1‑251 (MQ=255)                                                                                                                                                                                                                                 gCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGc                     >  1:320325/1‑251 (MQ=255)                                                                                                                                                                                                                                      tGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCa                 <  2:415120/250‑1 (MQ=255)                                                                                                                                                                                                                                          cATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGTCTGTAGAT‑ATGAACTTTAGGCAGCCACAGATCCGGTGCAAATTAGTTCGCCCTCACCCTCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCCTCCATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGatat            >  1:197718/1‑251 (MQ=255)                                                                                                                                                                                                                                          cATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCAAGCTCCGGTGCAAATTCGTTCGCCATCCCCATCGGCGGCCGATTACGCTGGTCCTCAATGCCGGCATCAAAAGGCAG‑CCTTGGACCCAGCAACAACATGCCATCAATTTGCTTGGTGGTGATCAAATCGATAAAGGTTTTTTCCTGCTGGTTTTGATGCGCACCGTCGCCAATctggcacaactat            >  2:88231/1‑239 (MQ=255)                                                                                                                                                                                                                                             ataaataaTTTTCTGCATCAAATGCGGCGGACAGATTGTCGAT‑ATGAACTCTAGGCAGCTCCCGCTCCGGTGCAAATTAGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCCTCAAACGGCAG‑CCTTGAACCCCGAAAAAACATGCCATCAATTTGCTTGGGGATGATCAAATCGATAAAGGGTTTTTTCTGCTGATTTTGATGAGAACAATAGCCAATCAGcccccgaaatccc         >  2:456053/1‑239 (MQ=255)                                                                                                                                                                                                                                                   aatTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAATTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTGAttt   <  2:456425/251‑1 (MQ=255)                                                                                                                                                                                                                                                   aatTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTGAtt    <  2:384156/250‑1 (MQ=255)                                                                                                                                                                                                                                                    atTTTCTGCATCAAATGCGGCGGTCAGATTTTCGAA‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGGGATGATCAAATAGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCACCACCAGATATCCCTGATTTg  >  1:362225/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         GAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGAT‑ATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTGATTTG  >  NC_000913/4123683‑4124176

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