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%	238.5 / NA	62	V176E (GTA→GAA)	cytR	DNA‑binding transcriptional repressor CytR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (29/33);  total (29/33)

AGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGAT  >  NC_000913/4123684‑4124164                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             aGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTTGTGACACAGCGGCCTCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGGCCTTGGTCATATAAATAATTTTCTg                                                                                                                                                                                                                                           >  2:187938/1‑250 (MQ=255) aGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTg                                                                                                                                                                                                                                           <  1:187938/250‑1 (MQ=255)         cacaTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAAt                                                                                                                                                                                                                                   <  1:163712/250‑1 (MQ=255)              cGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATgcggcg                                                                                                                                                                                                                             <  1:407827/251‑1 (MQ=255)               gCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCCGCTATACAGCCAACCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCaaaagacgcgg                                                                                                                                                                                                                            >  2:233795/1‑243 (MQ=255)                 tATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACAAAGCCTTGCAGGCGGTAGTGACACAACGGCATCTCTTCGGGGCCGGCTATACAGCCAATCCGCTTATGCCCTTGATCATATAAATAATTTTCTGCATCAAATGCGGCGGTc                                                                                                                                                                                                                          >  1:237221/1‑251 (MQ=255)                  aTGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCCAATGCGGCGGTc                                                                                                                                                                                                                          >  1:373359/1‑250 (MQ=255)                       agaagaCAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGGCAGAttt                                                                                                                                                                                                                    >  1:37398/1‑250 (MQ=255)                           gaCAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAAACCCTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGcggcgggcacaattgtcc                                                                                                                                                                                                                >  2:365173/1‑239 (MQ=255)                             cagcagTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATTCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGttcgatttttggt                                                                                                                                                                                                               >  1:490326/1‑238 (MQ=255)                                    aGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACt                                                                                                                                                                                                        <  1:281792/250‑1 (MQ=255)                                     ggCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCCTCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGcggcgggccgggtttcggtcttaaatgg                                                                                                                                                                                                      >  1:428572/1‑229 (MQ=255)                                            gtgGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAACGGCATCTCTTCGGGACCGGCTATACAGCCCATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGGCAgcatgtcggtgttaaatttaaggccc                                                                                                                                                                                               >  2:192280/1‑226 (MQ=255)                                                     cAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCCGATTGCTGAAATGAACTGTAGGCagctcccgctcc                                                                                                                                                                                      >  2:134278/1‑251 (MQ=255)                                                         cagcTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCCGCTCCGGTg                                                                                                                                                                                  >  1:518218/1‑251 (MQ=255)                                                         cagcTGCGGCATCGCTGTGCTTCCGGCTTCGAAGGTGAAGGCGCCACGGGCGCTGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGGGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTg                                                                                                                                                                                  <  1:402594/251‑1 (MQ=255)                                                                 gcATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAAtt                                                                                                                                                                            <  1:192280/249‑1 (MQ=255)                                                                     cGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAAttcgttc                                                                                                                                                                       <  1:321141/250‑1 (MQ=255)                                                                              tgcGGCTTCGAAGGTGACGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCAc                                                                                                                                                               <  2:155777/247‑1 (MQ=255)                                                                              tCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCAcc                                                                                                                                                              <  1:260936/250‑1 (MQ=255)                                                                                             gAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTAc                                                                                                                                              <  2:37398/251‑1 (MQ=255)                                                                                                    ccACGGGCGATGGAATGGGGATCAACCATAATGCCGCAGGCG‑CGCAGCGCCTGAAAATAGCCTTGGAGGCGGGAGTGACACAGCGGCAACTCTTCGGGACCGGCTAAACAGGCAATCCGTTTATGCCCTTGGTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGGCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAAATCGTTCGCCATCACCATCGGGGGCAGATTACGCTGTTc                                                                                                                                       <  2:522589/251‑1 (MQ=255)                                                                                                      aCGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCt                                                                                                                                     <  1:79208/251‑1 (MQ=255)                                                                                                             aTGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCt                                                                                                                              <  1:220145/251‑1 (MQ=255)                                                                                                                   tGCGGATCAACCATAATGACGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAAATTAACTGTAGGCAGCTTCAGCTCCGGGGGAAATTTGTTCGCCATCCCCATTGGCGGCAGATTACGCCGTTCTCCAATGCTGGCATc                                                                                                                        >  1:171375/1‑251 (MQ=255)                                                                                                                     cGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGGGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATcc                                                                                                                       >  1:178080/1‑249 (MQ=255)                                                                                                                     cGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGGGACACCGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCCTATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACAGTAGGGAGCTCCAGCTCCGGTGCAAATTCGTTTGCCATCACCATCGGCGGCAGATTACGCTTTTCCTCAATGCTGGCATCa                                                                                                                       >  1:503368/1‑250 (MQ=255)                                                                                                                       gATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAAc                                                                                                                    <  1:14193/251‑1 (MQ=255)                                                                                                                         tCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCCTATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCCCCAGCTCCGGTGCAAATTCGGTCGCCATCACCATCGGCGGGAGATTACGCTGTTCCCCAATGCTTGCCTCAAACgg                                                                                                                  >  2:230057/1‑251 (MQ=255)                                                                                                                             ccATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAgn                                                                                                              <  2:275398/251‑2 (MQ=255)                                                                                                                                        cAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGGTATAAACTGGAGGCAGCTCCAACCCCGGGGCAAATTCGTTTGCCAACACCATCCGCGGCAGATTTCGCTGTTTCTCAATGCTGGCATCAAACGGCAA‑CCTTGAAccc                                                                                                    >  1:124041/1‑249 (MQ=255)                                                                                                                                         aG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCCTCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAAGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGGGGCAGATTACGCTGTTCCTCCATGCTGGCATCAAACGGCAG‑CCTTGAAccccgc                                                                                                 >  1:470005/1‑248 (MQ=255)                                                                                                                                                gCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGGCAGATTGTCGATATGAACTGGTGGCAGCCCCAGCCCCCGTGCAAATTCCTTCGCCATCACCACCGGCGGGCGATTACGCTGTTCCCCAATTCTGGCATCCAACGGCAGGCCTT‑AACCCaacaacca                                                                                            >  2:28756/1‑248 (MQ=255)                                                                                                                                                           aCATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTTTCGATATGAACTGTAGGCCGCTCCAGCTCCGGTGCCAATTCGTTCGCCATCACCATCGGCGGCAAATTACGCTGTTCCTTCATTCTTGCATCCAACGGCAGGCTTTGAAccagccaacaccagccatccat                                                                                 >  1:17406/1‑231 (MQ=255)                                                                                                                                                                     gCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGt                                                                       <  2:178080/250‑1 (MQ=255)                                                                                                                                                                          gggTAGTGAAACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTgatgat                                                                 <  2:17406/250‑1 (MQ=255)                                                                                                                                                                           ggTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATc                                                                <  1:99162/251‑1 (MQ=255)                                                                                                                                                                            gTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCa                                                               <  1:215838/251‑1 (MQ=255)                                                                                                                                                                              aGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCaaa                                                             <  2:373359/251‑1 (MQ=255)                                                                                                                                                                              aGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCaaa                                                             <  2:203021/251‑1 (MQ=255)                                                                                                                                                                               gTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAAt                                                            >  1:244082/1‑251 (MQ=255)                                                                                                                                                                                  acacAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTTGCATCAAACGGCAG‑CCTTGAAACCAGCAACAACATGCCATCAATTTGCTTGGGGATGATCAAATCg                                                          >  2:480850/1‑250 (MQ=255)                                                                                                                                                                                        cGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAgg                                                  <  1:166885/251‑1 (MQ=255)                                                                                                                                                                                        cGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAgg                                                  <  1:166718/251‑1 (MQ=255)                                                                                                                                                                                         ggCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGGTCCTCAATGCTTGCATCAAACGGCAG‑CCTTTAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAggg                                                 >  2:282808/1‑250 (MQ=255)                                                                                                                                                                                           cATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGttt                                               <  2:518218/251‑1 (MQ=255)                                                                                                                                                                                                    gggACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCctgct                                      <  1:282808/251‑1 (MQ=255)                                                                                                                                                                                                    gggACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCctgct                                      <  1:365173/251‑1 (MQ=255)                                                                                                                                                                                                     ggACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCctgctg                                     <  1:220474/251‑1 (MQ=255)                                                                                                                                                                                                      gACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCctgctg                                     >  1:379005/1‑250 (MQ=255)                                                                                                                                                                                                       aCCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGAt                                   <  2:428572/251‑1 (MQ=255)                                                                                                                                                                                                         cGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCGCCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAATGTTTTTTCCTGCTGAttt                                 <  2:592056/251‑1 (MQ=255)                                                                                                                                                                                                          ggCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAAGCAGCTCCAGCTCCGGGGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCCGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAAAATGCCATCAATTTGCTTGGGGCTGATCAAATCGATCAAGGGTTTTTTCC‑GCTGAgttt                                >  2:146685/1‑251 (MQ=255)                                                                                                                                                                                                              ataCAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGATTTTGAt                             <  2:244082/250‑1 (MQ=255)                                                                                                                                                                                                                       aTCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATAAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTc                   <  2:124041/251‑1 (MQ=255)                                                                                                                                                                                                                         ccGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCg                  <  2:8729/250‑1 (MQ=255)                                                                                                                                                                                                                               aTGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGATGTTGATGCGCACAGTCGCCAATc            <  2:508693/250‑1 (MQ=255)                                                                                                                                                                                                                                  cccTTGCTCATATAACTAATTTTCTGCATCAAATGCGGCGGTCAGGTTGTCGATATGAAATGTAGGCAGCTCCAGCCACGGGGCAAATTCGTTCGCCATCACCATCCGCGGCAGCTTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTTATTGGTGATGATCCAATCGAT‑AAAGGTTTTTTCCTGCTTAtttttattgcagcagacgcccattcgc         >  2:169337/1‑227 (MQ=255)                                                                                                                                                                                                                                     ttGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCAcc      >  1:516514/1‑250 (MQ=255)                                                                                                                                                                                                                                     ttGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCCTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAAACGAT‑AAAGGGTTTTTCCTGCTGCTTTTTATGCGCACAGTCGCCAATCAGCAccc     >  2:410430/1‑250 (MQ=255)                                                                                                                                                                                                                                       gCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCAcc      >  1:121953/1‑248 (MQ=255)                                                                                                                                                                                                                                         tCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCCTCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCCGat  >  2:492358/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             AGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAG‑CGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAG‑CCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAT‑AAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGAT  >  NC_000913/4123684‑4124164

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