Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,735	T→C	P186P (CCT→CCC)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,735	0	T	C	97.4%	122.9 / ‑5.1	38	P186P (CCT→CCC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (16/21);  minor base A (1/0);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGC  >  NC_000913/4638490‑4638978                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    aTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGGACCGGGTTAGAGCAAGTCGGGCAGGGCGTGCATATCCCGTTACCCGAAgg                                                                                                                                                                                                                                                <  2:82772/251‑1 (MQ=255)    gggATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGa                                                                                                                                                                                                                                             <  1:302400/251‑1 (MQ=255)          tataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGGAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCGGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGc                                                                                                                                                                                                                                        <  1:227906/250‑1 (MQ=255)                ggTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTAGATGTTTCACGTCTTAGCGAGCTAAACGCGCCAAATATCATCTTAGGCCAGCCCTTTATTGTGATTAGCGACGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCAGAAAGGGACTGGAGGaacc                                                                                                                                                                                                                                  >  1:290543/1‑248 (MQ=255)                             acagtgagttcacgctatacgCCGATTTCGATGTTTCGCGTCTTACCGAGCTCAACGCCCCAAATATCATATTAGTCACGCCCTTTATTGTGATTAGCGTCGGGGCTGCGCGTGGTAT‑TGGCGTGGTGATAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAACAGACGCTGAAGCTg                                                                                                                                                                                                                    >  2:95029/21‑251 (MQ=255)                                       aaCGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGccc                                                                                                                                                                                                          >  1:431291/1‑251 (MQ=255)                                               aaGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGCACGGCGCTGACCCTTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGAGTGCATATACCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTAAAGCTTAATATGGCCCTGAAtttt                                                                                                                                                                                                  >  2:557990/1‑250 (MQ=255)                                                aGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGCACGGCGCTGACCATTGAACCCGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAGGGGGACTGGCGGAAGCAGAACCTGGAGccggataatgccctgtatata                                                                                                                                                                                                  >  2:110076/1‑230 (MQ=255)                                                         cGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTAcc                                                                                                                                                                                         <  2:345456/250‑1 (MQ=255)                                                          gATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCACGCCCTTTATTGTGATTAGCGTCAGGGATGCGCGTGGCAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCGGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCCGGGCGTGCATATCCCTTTACCCGAAGGGGCCTGGCGGAAGCAGAAACTGAAGCTGGATATGGCCCTGAATTTAAGCGGTGCCgg                                                                                                                                                                                       >  1:234160/1‑251 (MQ=255)                                                                    ggctgaGCGAGCTCAACGCGGCAAATATCATCTTAGGCAAGCCCTTTAGTGTGATTAGCGTCGGGGATGCGCGTGTTATGTGG‑GTGGTGAAAGCGCCAGAAGTTAACGGAACGGCGCTGACCATTGAGCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGGATATCCCGTGATCCGAAGGGGTCTGTCGGAGGCAGAGCCTGAAGCTCAATATGGGCCTGAATGTGAGCGGTCCGGGCGCTCTTTCt                                                                                                                                                                             <  2:538537/246‑1 (MQ=255)                                                                      cTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCCTTGAACCCGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGGAGCAGAACCTGAAACTGAATAATGCCCTGAATTTAAAGGGTACACGCGATCCTTCtgg                                                                                                                                                                           >  2:354458/1‑250 (MQ=255)                                                                                    cgcgCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCCTATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGGATATGGCCCTGAATTTTAGCGGTACCGGCGATCTTTCTTTGGGGCCTgccgggg                                                                                                                                                             >  2:155259/1‑250 (MQ=255)                                                                                            atatCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGc                                                                                                                                                     >  1:236351/1‑251 (MQ=255)                                                                                              atcatcTTAGGCAAGCCCTTTATTGTGAGTAGCGTCGGAGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGTTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGa                                                                                                                                                   <  2:170958/251‑1 (MQ=255)                                                                                                   cTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATg                                                                                                                                               >  1:160885/1‑250 (MQ=255)                                                                                                          aaGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGCCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATTTCCCTTTACCCGAAGGGGCCTGGCGGCAGCCGAAACTGACGCTGAATGTGGCCCTGAATTTAAGCGGTACCGGGGGTCTTTCTTTGGTGCCTGCCGGGCGGAATAGCGAAATGACCTTAAc                                                                                                                                       >  2:611817/1‑251 (MQ=255)                                                                                                                 ttATTGTGATTAGCGGCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCaa                                                                                                                                 <  1:354458/250‑1 (MQ=255)                                                                                                                   aTTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTg                                                                                                                              <  1:507680/251‑1 (MQ=255)                                                                                                                     tgtgATTAGCGTCGTGGAGGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCGGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCATCTgg                                                                                                                             <  2:409751/250‑1 (MQ=255)                                                                                                                            aGCGTCGTGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGGTACCGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGTGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCGGAATTTAAGCGGTACCGGCGATCTTTCGGTGGTGCCTGCCGGTCGTAATAGCGAAATGACCTTAACCAGCAGCTTGCCGCATcc                                                                                                                     <  2:582033/251‑1 (MQ=255)                                                                                                                               gTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAACGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGCATATGGCCCTGAATTTAAGCGGTACCGGCGCTCTTTCTGTGGTGCATGCCGGGCGTAATCGCGAAATGACCTTAACCAGCAACTGTCCGcctcac                                                                                                                    >  1:329079/1‑244 (MQ=255)                                                                                                                                 cGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTTCCTGCCGGGCGGAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGtt                                                                                                                >  2:28416/1‑251 (MQ=255)                                                                                                                                    ggATGCGCGTGGTAT‑TTGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTGAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGTACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGTGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGttttt                                                                                                             <  2:587707/251‑1 (MQ=255)                                                                                                                                                                                     ggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTc                                                             <  2:196988/251‑1 (MQ=255)                                                                                                                                                                                            aCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGgcgc                                                       >  1:470130/1‑250 (MQ=255)                                                                                                                                                                                                  gAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGAGCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTGTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCa                                                <  1:572509/251‑1 (MQ=255)                                                                                                                                                                                                    aCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCaaa                                              <  2:337832/251‑1 (MQ=255)                                                                                                                                                                                                     cccGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCTCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAag                                             <  1:611817/251‑1 (MQ=255)                                                                                                                                                                                                       cGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTGTACCCGAATGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAACGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCAGCATCCCAGGTGGTTAGGTGATTTTTTGCCAGCCAAACTGAAAGGTAGCGAGTCAGGTTGTTAGGCGCAATGGCAAagca                                           <  2:229981/251‑1 (MQ=255)                                                                                                                                                                                                        ggCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTGTCAGGCGCAATGGCAAagcag                                          <  1:213149/251‑1 (MQ=255)                                                                                                                                                                                                               ggTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGttt                                   <  2:535000/251‑1 (MQ=255)                                                                                                                                                                                                               ggTTAGAGCAAGGCGGGCAGGGCGTGCATAGCCCTTTACCCGAAGGGGACGGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTTGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTGTTTAGGTGATGTTTTGCCAGCCAAACGGGAAGTTATCGAGTCAGGTGGGCAGGCGCAATGGCAAAGCAACTGGttt                                   <  2:362060/251‑1 (MQ=255)                                                                                                                                                                                                                gTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTg                                  <  1:370931/251‑1 (MQ=255)                                                                                                                                                                                                                    gagCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGGTATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCta                               >  2:98825/1‑250 (MQ=255)                                                                                                                                                                                                                        aaGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTaataat                          >  1:274633/1‑251 (MQ=255)                                                                                                                                                                                                                             gggCAGGGCGTGCATATCCCTTTACCCGCAGGGGACTGGCGGAAGCAGAACCTGAAGCTGCATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCCGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCgg                     >  2:74660/1‑251 (MQ=255)                                                                                                                                                                                                                                                  ttACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGc  <  1:174587/249‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGC  >  NC_000913/4638490‑4638978

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