Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,560,309	T→A	P124P (CCA→CCT)	yjiH ←	uncharacterized protein YjiH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,560,309	0	T	A	95.7%	145.5 / NA	47	P124P (CCA→CCT)	yjiH	uncharacterized protein YjiH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (20/25);  minor base C (1/1);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACG  >  NC_000913/4560071‑4560549                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 cAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCcaac                                                                                                                                                                                                                                      >  1:372027/1‑251 (MQ=255)        aaaCACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCCTCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGtt                                                                                                                                                                                                                               >  1:208793/1‑251 (MQ=255)               aGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAATGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGcac                                                                                                                                                                                                                        <  1:193793/251‑1 (MQ=255)                ggATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGcac                                                                                                                                                                                                                        >  1:336971/1‑250 (MQ=255)                  aTGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGAGACAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAgg                                                                                                                                                                                                                     >  1:265069/1‑251 (MQ=255)                            aaGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAAc                                                                                                                                                                                                           <  1:286398/251‑1 (MQ=255)                                 acgagatctgaagggCCCAGAGATGCAAACCACGCGTCGACGGAAGCGAAGTATGTGGTGATGATGGCGTGATCGCTGGTTTGATATGTGGCGAAAATCGCTTTATCGCGTTCGGTAATTTCACCTGCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGTCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCACAACAACGGTTTTAGCACAGGTATCATTAActgctg                                                                                                                                                                                                      <  2:103721/237‑1 (MQ=255)                                 aCGATTACGGGCGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGACCAACCACGGTTTTAGCACAGGTGTCATTAActgctg                                                                                                                                                                                                      >  2:343811/1‑251 (MQ=255)                                  cGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACtgctgc                                                                                                                                                                                                     <  1:439362/251‑1 (MQ=255)                                    aTTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGTTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCCGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCCAAGGGGTCATTACctgc                                                                                                                                                                                                        >  2:122253/1‑246 (MQ=255)                                    aTTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGGTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAActgc                                                                                                                                                                                                        <  1:122253/246‑1 (MQ=255)                                        cGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcggcg                                                                                                                                                                                               >  1:19919/1‑251 (MQ=255)                                              tGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCAtt                                                                                                                                                                                                              >  1:2516/1‑230 (MQ=255)                                              tGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCAtt                                                                                                                                                                                                              <  2:2516/230‑1 (MQ=255)                                                  aaGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGCTATGCGGCGATAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGTTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGCTCAGCGCCAGCGAGCCGATCCCAGGACTGCCCAACAACGGTTTTAGCACAGGAGTCATTAACTGCTGCCCGGCGCGTAATcca                                                                                                                                                                                      >  2:30639/2‑250 (MQ=255)                                                           aaaCACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTTCCGCTGGTTTGATATGCGGCGAAAATCTCTTTATCGCGTTCGGTAATTTCGCCTTCCTGCGCCAGCGCTTTGGTCATACCGGCGGCGGCACCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCAt                                                                                                                                                                            <  2:179887/251‑1 (MQ=255)                                                            aaCACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATc                                                                                                                                                                           <  1:64528/251‑1 (MQ=255)                                                                   ggACGGCGGAAGAGAAGTAGTTGGTGATGATGGCGCCACCGCTGGGTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGAGCAGCGCCAGCGAGCAGATCCCAGGAAAGCCCAACAA‑GGTTTTAGCACAGGTGTCATTAACTGCGGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTg                                                                                                                                                                   <  2:213041/250‑1 (MQ=255)                                                                            aaGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATc                                                                                                                                                           >  1:230610/1‑251 (MQ=255)                                                                            aaGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATc                                                                                                                                                           <  1:16418/251‑1 (MQ=255)                                                                                           tgctgatgGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCATGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTATCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGCCCATCCGTAATTGCAATAATCACCagaggggcaat                                                                                                                                             >  2:408981/4‑243 (MQ=255)                                                                                              tgatgGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATc                                                                                                                                                           <  1:158888/233‑1 (MQ=255)                                                                                              tgatgGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGGTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGCCCATCGGTAATTGAAATAATc                                                                                                                                                           >  2:158888/1‑233 (MQ=255)                                                                                                  ggCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCCGGAATGCCCCACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACg                                                                                                                                     >  2:435728/1‑251 (MQ=255)                                                                                                        aCCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCTGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCa                                                                                                                               <  1:201621/251‑1 (MQ=255)                                                                                                                      atGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACCACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAACTAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCCCAAGCGCaa                                                                                                                   >  2:307746/1‑249 (MQ=255)                                                                                                                                  tCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACgg                                                                                                                                    >  1:473487/1‑220 (MQ=255)                                                                                                                                  tCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACgg                                                                                                                                    <  2:473487/220‑1 (MQ=255)                                                                                                                                   cgcTTTATCGCGTTCGGTAATTTCACCTTCCGGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACGGCTGCGCGGCGCGTAATccaccat                                                                                                                                                                                  <  2:159507/171‑2 (MQ=255)                                                                                                                                   cACTTTATCGCGTTCGGTAATTTCACCTTCCTGAGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGAGTCATTAACTGCTGCGCGGCGCGTAATCCACCaa                                                                                                                                                                                  >  1:159507/1‑173 (MQ=255)                                                                                                                                    cctgtatctctgtCGGTAATTTCTCCTTCCCGCGCCAGCTCTTTTGTTCTACCGTGGGGTGCAACCGTTGTTTGCAGGGTTGCGGTCAGCGCCAGCGAGCAGATTCCAGGAATGCCCAACAACGGGTTTAGCCCTGGTGGCATTAACTGCTGCGCGGCGCCTAATCCCACAATGCCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCt                                                                                                     <  2:230610/237‑1 (MQ=255)                                                                                                                                    ccTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGTATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCttt                                                                                                   <  1:245339/250‑1 (MQ=255)                                                                                                                                        tATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGc                                                                                                <  1:339623/250‑1 (MQ=255)                                                                                                                                                  ggTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCAc                                                                                     >  1:346049/1‑251 (MQ=255)                                                                                                                                                            ccTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACccgcgccgcg                                                                           >  1:74485/1‑251 (MQ=255)                                                                                                                                                                                ggTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGt                                                       <  2:19558/251‑1 (MQ=255)                                                                                                                                                                                   cctACCGGCTGCGGCACCGGGGTTTTGCAGGTCAGCGAACAGCGCCAGCGAGCAGATCCCAGGAATGTCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCTCGTAATCCTCCAAGACCATCGGTAATTGAAATAACCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCTCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGCGGTGGTGGCTCCGTTc                                                     <  2:392106/248‑1 (MQ=255)                                                                                                                                                                                             gcggcATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTGAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCTAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAg                                          <  1:435728/251‑1 (MQ=255)                                                                                                                                                                                              cggcATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAg                                              >  1:238063/1‑246 (MQ=255)                                                                                                                                                                                              cggcATCGGTGTTTTGCAGGTGAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAg                                              <  2:238063/246‑1 (MQ=255)                                                                                                                                                                                                 cATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGc                                      >  1:103300/1‑251 (MQ=255)                                                                                                                                                                                                             gCAGGTTAGCGAACAGCGCCAGCAAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACATGTGTCATTAACTGCTGCGCTGCGCGTAATCCACCAAGACCATCGTTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGGTGCCCAAACGATcc                           <  2:443133/250‑1 (MQ=255)                                                                                                                                                                                                               aGGTTAGCGATCAGCGCCAGCGAGCAGAGCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCAGCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGtt                        <  2:183040/251‑1 (MQ=255)                                                                                                                                                                                                                       gATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACccg                                                                                  >  1:121200/1‑185 (MQ=255)                                                                                                                                                                                                                       gATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACccg                                                                                  <  2:121200/185‑1 (MQ=255)                                                                                                                                                                                                                                  gCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCAACAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCaaaa     <  2:74485/251‑1 (MQ=255)                                                                                                                                                                                                                                     cgCAGATCCCCGGAAGGCCCAACAACGGTGTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGTCTTACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCCCCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAAcg  <  2:265069/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 CAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACG  >  NC_000913/4560071‑4560549

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