Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	358,962	G→A	G55D (GGC→GAC)	cynT →	carbonic anhydrase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	358,962	0	G	A	95.6%	208.0 / ‑5.1	69	G55D (GGC→GAC)	cynT	carbonic anhydrase 1
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (46/19);  minor base C (2/1);  total (48/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AACTCATCATAAATGAGACTTTTACCTTATGACAATCGGCGAGTAGTCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCG  >  NC_000913/358717‑359205                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    aaCTCATCATAAATGAGACTTTTACCTTATGACAATCGGCGAGTAGTCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCCACGTGAGCCTGACGATc                                                                                                                                                                                                                                                >  1:154852/1‑251 (MQ=255)  aCTCATCATAAATGAGACTTTTACCTTATGACAATCGGCGAGTAGTCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCt                                                                                                                                                                                                                                               <  1:320595/251‑1 (MQ=255)        catAAATGAGACTTTTACCTTATGACAATCGGCGAGTAGTCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGCGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGTGCTGGTGACGCAACGTGAGCCTGACGCTCTGTTCGt                                                                                                                                                                                                                                         >  2:174885/1‑251 (MQ=255)            aaTGAGACTTTTACCTTATGACAATCGGCGAGTAGTCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTAtt                                                                                                                                                                                                                                     >  1:65295/1‑251 (MQ=255)                  aCTTTTACCTTATGACAATCGGAGAGTAGTCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGAGGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCAACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAAc                                                                                                                                                                                                                               >  2:260470/1‑251 (MQ=255)                        aCCTTATGACAATCGGCGAGTAGTCTGCCTCTCATTCCCGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAACTTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGc                                                                                                                                                                                                                         >  2:292688/1‑251 (MQ=255)                              tggcAATCGGCGAGTAGTCTGCCTCTCATTCCAGAGACAGACAGAGGTTCACGGTGAAAGAGATTATTGCTGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGCAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACCCTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCGGAGCTGGTGCCGCAACATGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATc                                                                                                                                                                                                                   >  2:317886/4‑251 (MQ=255)                              tGACAATCGGCGAGTAGTCTGCCTCTCATTCCAGCGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGAGAGCCTGACGATCTGTTCGTTATTCGCCACGCGGGCAATATc                                                                                                                                                                                                                   >  2:258161/1‑251 (MQ=255)                              tGACAATCGGAGAGTAGTCTGCCTCTCATTCCAGAGACAGACCGAGGTTAACGGTGAAAGCGCTTATTGATGGATTCCTTAAATTCCCGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAACCAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGCCGCAACGTGAGCCTGACGATCTGTTCGGTATTCGCAACGCGGGCAatatt                                                                                                                                                                                                                   >  2:57137/1‑250 (MQ=255)                                     cGGCGAGTCGTCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGCAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGCTCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCtt                                                                                                                                                                                                            >  2:556481/1‑251 (MQ=255)                                     cGGCGAGTAGTCTGCCTCTCATTCCAGAAACAGACAGAGGATAACGGTGAAAGAGATTATTGATGGATTCCTTAACTTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGCCACAGCAACGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGCGCCTTACGATCTGTTCGTTATTAGCCACGCGGGCAATATCGTCACtt                                                                                                                                                                                                            >  2:357637/1‑251 (MQ=255)                                           gtagtCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACg                                                                                                                                                                                                      <  1:104959/251‑1 (MQ=255)                                             agtCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACgg                                                                                                                                                                                                     <  2:222700/250‑1 (MQ=255)                                              gtCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAGAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGc                                                                                                                                                                                                   >  2:518704/1‑251 (MQ=255)                                                      ctcATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAAc                                                                                                                                                                                             <  1:372909/249‑1 (MQ=255)                                                           tCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTg                                                                                                                                                                                       <  1:160723/250‑1 (MQ=255)                                                             cAGAGACAGACAGAGGTTCACGGTGTAACAGATTATTGCTGGATTACTCCAATTGCAGCGCGAGGCATTTCCGAAGCGAGAAGCCTTGTTTAATCAGCCTGCGACACAGCAAAGCCCGCGCACACCTTATATCTCCTGCTCCGACAGCCGTCTGGTCCCTCATCTGGTGACGCAACGCGAGACTCCCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCCACGGGCCGGATCCCGgcggcg                                                                                                                                                                                    >  2:558031/1‑251 (MQ=255)                                                                              tAACGGTGAAAGAGATTCTTGCTGGCTTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGATCGCGTTTCTGCTTCGGTgg                                                                                                                                                                     >  2:295715/1‑249 (MQ=255)                                                                                 cGGTGAAAGAGCTTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTAt                                                                                                                                                                >  1:146322/1‑251 (MQ=255)                                                                                       aaGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTc                                                                                                                                                          <  1:206656/251‑1 (MQ=255)                                                                                          agaTTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGTGTATGCCGTCGCt                                                                                                                                                       >  2:451732/1‑251 (MQ=255)                                                                                          agaTTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCt                                                                                                                                                       <  1:518704/251‑1 (MQ=255)                                                                                              tattGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTgcg                                                                                                                                                    <  2:340475/250‑1 (MQ=255)                                                                                                 tGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGGGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTc                                                                                                                                                >  1:320022/1‑251 (MQ=255)                                                                                                  gATGGATTCCTTCAATTCCAGCGAGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGCGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTCTGCCGTAGCTGCGCCTCg                                                                                                                                               >  2:59360/1‑251 (MQ=255)                                                                                                  gATGGATTCCTAAAATTCCAGCGCGAGGCATTTCCCAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACGCAGCAAAGCCCGCGCGCACTTTTTATCTCCTGCCCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTTACGAGCTGTTCGTTATTCGCAACTCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGTTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTc                                                                                                                                                <  2:523255/250‑1 (MQ=255)                                                                                                    tGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACCGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCggg                                                                                                                                             >  1:98135/1‑251 (MQ=255)                                                                                                     ggATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTT‑ATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGt                                                                                                                                            <  2:292497/250‑1 (MQ=255)                                                                                                      gATTCCTTAAAGTCCAGCGCGAGGCATTGCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTa                                                                                                                                           <  2:356248/251‑1 (MQ=255)                                                                                                        ttCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGAAACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATc                                                                                                                                         >  1:49354/1‑251 (MQ=255)                                                                                                                   ccAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCAGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATtgtg                                                                                                                               >  1:355051/1‑250 (MQ=255)                                                                                                                       cgcgAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCTTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTtg                                                                                                                          <  1:306771/251‑1 (MQ=255)                                                                                                                               aTTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGCTTTGTGGTCAtt                                                                                                                  >  2:173313/1‑251 (MQ=255)                                                                                                                                tttCCGAAGCGGGACGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATTTGCCATTGTGATTTGTGGTCATTc                                                                                                                 >  2:381426/1‑251 (MQ=255)                                                                                                                                       aGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGAGCTGTTCGTTATTCGCAACGCGTGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTGCGGTGGAGTATGCCGGCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACtgt                                                                                                          <  1:575660/251‑1 (MQ=255)                                                                                                                                        gCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCGCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCACCtgtg                                                                                                         <  2:427238/251‑1 (MQ=255)                                                                                                                                                     tttAAACAGCTGGCGACACCGCAAAGCCCGCGCACACTTTTTATCTACTGCTCCGACAGCCGTCTGCTCCCTGAGCTGGAGACACAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGCCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCACCTGTGGCGCGATGACCGc                                                                                            >  1:507451/1‑251 (MQ=255)                                                                                                                                                        aaaCAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGCCGATCTGTTCGTTATTCGCACCGCGGGCAATATCGTCCATTCCTACGGGCCGGAAACCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCAt                                                                                         >  1:528035/1‑251 (MQ=255)                                                                                                                                                          aCAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTg                                                                                       >  1:417886/1‑251 (MQ=255)                                                                                                                                                           cAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGc                                                                                      >  1:189453/1‑251 (MQ=255)                                                                                                                                                            aGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGcc                                                                                     >  2:297842/1‑251 (MQ=255)                                                                                                                                                                 gCGACACAGCAAAGCCCGCGCACTCTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTg                                                                                <  1:38337/251‑1 (MQ=255)                                                                                                                                                                    acacAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTc                                                                              >  2:166358/1‑250 (MQ=255)                                                                                                                                                                     cacaGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGAGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATTACCGCCATTGCCAGCTGTc                                                                              >  1:208420/1‑249 (MQ=255)                                                                                                                                                                           aaaGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATg                                                                      >  1:36777/1‑251 (MQ=255)                                                                                                                                                                                ccGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGCTCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGAcc                                                                  >  1:380201/1‑250 (MQ=255)                                                                                                                                                                                   cgcACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCatat                                                              >  1:194495/1‑251 (MQ=255)                                                                                                                                                                                        acTTTTTATCTCCTGCTCCGCCAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATAtgcctg                                                         >  1:489906/1‑251 (MQ=255)                                                                                                                                                                                              tATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCCGTGCATGGACCATATGCCTGCCGtct                                                   >  1:19985/1‑251 (MQ=255)                                                                                                                                                                                               aTCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGtctc                                                  >  2:550884/1‑251 (MQ=255)                                                                                                                                                                                                 ctcCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTccc                                                >  2:225334/1‑251 (MQ=255)                                                                                                                                                                                                  tcCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGACGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTccc                                                >  1:132591/1‑250 (MQ=255)                                                                                                                                                                                                   ccTGCTCCGACAGCCGTCGGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTGCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTccc                                                <  1:292688/249‑1 (MQ=255)                                                                                                                                                                                                     tGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCCCTg                                            >  2:385692/1‑251 (MQ=255)                                                                                                                                                                                                      gCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACt                                             <  1:206505/249‑1 (MQ=255)                                                                                                                                                                                                        tCCGACAGCCGTCTGGTTCCTGAGCTGGTGACGCACCGTGAGCCTGAAGATTTGTTCGTTATTCCCAACGCGCGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCCTTTCTACTTCGGTGGAGTCTGCCGTCGATCCGCTTCGGGCATCTGCCATTGCGATTTGTGGTCATCCCCACTGCGGAGCGATGACCGCCATTGCCAGCTGTACGGGCATGGACCATATGCCTGCCGTCTCCCCCTgccc                                         >  2:343575/1‑248 (MQ=255)                                                                                                                                                                                                         ccGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTg                                        >  2:78048/1‑251 (MQ=255)                                                                                                                                                                                                             cAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGtt                                    >  1:416500/1‑251 (MQ=255)                                                                                                                                                                                                             cAGCCGTCGGGTCCCTGATCTGGTGACGCAACGTGAGCCTGACGATCTGGTCGTTAGTCGCAACGCGGGCAATCTCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCGGCTTCGGTGGGGTATTCCGTCGCTGCGCTTCGTGTATCTGACATGGTGATTTGTGTTCATGCCAACGGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGtt                                    <  1:557939/251‑1 (MQ=255)                                                                                                                                                                                                                 cGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGc                                >  1:274243/1‑251 (MQ=255)                                                                                                                                                                                                                  gTCTGGTCCCTGGGGTGGGGACGCAACGTGAGCCTGCCGTTCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTGCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGcc                               <  2:65295/251‑1 (MQ=255)                                                                                                                                                                                                                     tGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGCGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCgct                            >  2:493827/1‑249 (MQ=255)                                                                                                                                                                                                                               gCTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATAGTCCCTTCCTACGGGCCGGAACCCGGAGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGACGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCcgcg                  >  2:108556/1‑251 (MQ=255)                                                                                                                                                                                                                                cTGGTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGcgt                 >  2:239108/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggTGACGCAACGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGcgtc                <  1:78048/250‑1 (MQ=255)                                                                                                                                                                                                                                           aCGTGAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAgg      >  2:419972/1‑251 (MQ=255)                                                                                                                                                                                                                                            cgggAGCCTGACGATCTGTTCGGTATGCGCAACGCGGGCAATATCGTCCCTTCCTACTGGCCGGAACCTGGTGGCGTGTCTGCTTCGGTGGAGTATGCCGGCGCGGCGCTTCGGGTATCTGACATTGTGGTTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCt                                                                                 <  1:551690/172‑1 (MQ=255)                                                                                                                                                                                                                                            cGTGCGCCTGACGCTCTGTTCGTTATTCGACACGCGGGAAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTCTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGCATCTGACCTTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCt                                                                                 >  2:551690/1‑175 (MQ=255)                                                                                                                                                                                                                                               gAGCCTGACGATCTGTTCGTTATTCGCAACGCGGGCAATAGCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGGCGCTGCGCTTCGGGTATCTGACCTTGTGATTTGTGGTCATTCCAACTGTGGCGCGCTGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACAGGCTGCGTTATCCCGATTCGGCCCGCGTCGTTAATGAGgcgcg  >  2:449815/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    AACTCATCATAAATGAGACTTTTACCTTATGACAATCGGCGAGTAGTCTGCCTCTCATTCCAGAGACAGACAGAGGTTAACGGTGAAAGAGATTATTGATGGATTCCTTAAATTCCAGCGCGAGGCATTTCCGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCG  >  NC_000913/358717‑359205

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