Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,284	A→G	Q369Q (CAA→CAG)	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,639,284	0	A	G	90.5%	103.7 / 1.3	42	Q369Q (CAA→CAG)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (19/19);  total (21/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 = 8.82e-01

TCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATT‑GAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAAT  >  NC_000913/4639038‑4639528                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       tCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCgg                                                                                                                                                                                                                                                    >  1:480687/1‑251 (MQ=255)  cAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCgg                                                                                                                                                                                                                                                    <  1:38211/250‑1 (MQ=255)      aCCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGAGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCTTTGGAGATGTTTTATTTGCTCTTGCTGGAGCTTTCTGAACATATCGGATTTACCGTGGCATGGTTAATCGCCAGTCTGATTGTGGCGTTAATGAACGGTATTTATTTGCAGGCGGGATT‑g                                                                                                                                                                                                                                              >  2:190141/1‑251 (MQ=255)                ggACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGCTTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTAGAAAgatggg                                                                                                                                                                                                                                     >  2:404040/2‑246 (MQ=255)                 gACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATT‑GAAAGGTTGGTg                                                                                                                                                                                                                                   <  2:174892/251‑1 (MQ=255)                  aCCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTggcgc                                                                                                                                                                                                                                  >  1:232684/1‑248 (MQ=255)                   gcGGGCAACTAAGTACGCCACGCTGCTGATTGCTCTAACTTTTATGGCGTTTTTTGTATCTGAAACGCTCACTGCGCAACGTTTATCCCCACTGCAATATTCGCTGTCGCGGCTTCCATTGGCGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTTCCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAGTGAGAGATTTTTATTTGCAGGCGGTGTT‑GACAggaaggcgac                                                                                                                                                                                                                                 >  2:286148/2‑243 (MQ=255)                         gACTAAGGACGCCATTCTGCTGATTGCCCTGACTTTTATGACGTTCTTTGATTTTGAAACGCTCACAGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGGGATGTTATATTTGATCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCTTGGATAATCGCACGTCTTATTGGGGCGTTAATGACAGGTATTTTTTTGAAGGCGGTAtgctgtggattgggtaaaagct                                                                                                                                                                                                                            >  2:371925/1‑230 (MQ=255)                           cTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGGATT‑GAAAggcgggggcaaaacgc                                                                                                                                                                                                                           >  2:264511/1‑235 (MQ=255)                                tACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGGGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTCTTTGCAGGCGGGATT‑AAAAGGGTGGAGAAACCGCGTgttgtt                                                                                                                                                                                                                    >  2:441974/1‑251 (MQ=255)                                  cGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTa                                                                                                                                                                                                                  <  2:232684/251‑1 (MQ=255)                                        tCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTc                                                                                                                                                                                                             <  2:98773/250‑1 (MQ=255)                                           gctgATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTcg                                                                                                                                                                                                            <  2:331507/248‑1 (MQ=255)                                            ctgATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCCCGCt                                                                                                                                                                                                        >  1:208651/1‑251 (MQ=255)                                             tgATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑TAAAGGGTGGCGCAACCGCATGTTGttttaatcgggaggg                                                                                                                                                                                                       >  2:6983/1‑239 (MQ=255)                                             tgATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGttt                                                                                                                                                                                                                   <  1:6983/239‑1 (MQ=255)                                                   cACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGGTGGAGCAACAGACTGTTGTTTACCCTcgcgatgactgtg                                                                                                                                                                                                 >  1:90736/1‑242 (MQ=255)                                                    aCTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGc                                                                                                                                                                                                >  1:323368/1‑251 (MQ=255)                                                             tATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTTCCCTCGCGCAGTTGCTGCTgggatggg                                                                                                                                                                                       >  1:474480/1‑245 (MQ=255)                                                                          gTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACt                                                                                                                                                                          <  1:61603/251‑1 (MQ=255)                                                                                  aaCGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGGGGGGCTTTCATTTGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAAct                                                                                                                                                                  <  2:495646/251‑1 (MQ=255)                                                                                                   gTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTAACCTCGCGCTGTTGTTGCTGGATGGTGGGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCtg                                                                                                                                                 >  1:406261/1‑251 (MQ=255)                                                                                                     ttACACCCAATGCAATTTTTGCTGGTGGGGCTTTCATTGGTGATGTTTTAGTTGCTCTTGCTGGCGTTTTCTTAACATATCGTTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCtgt                                                                                                                                                <  2:323368/250‑1 (MQ=255)                                                                                                         ccccAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGttgtt‑g                                                                                                                                          <  2:205389/250‑1 (MQ=255)                                                                                                                 cAATATTTGCGGGTGGGTCTGTCATTGGTGATGTTTTAGTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTCCCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGCT‑GGGAACCAg                                                                                                                                  <  2:16395/251‑1 (MQ=255)                                                                                                                   atatTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGTTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAgtg                                                                                                                                <  1:264511/251‑1 (MQ=255)                                                                                                                       ttGCTGGGGGGGCTTTCGTTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACTGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCAAGTTGTTTTACCTCGCGGGGGTGTTGATGGATGGAGGGATGTGGGGACTGGTCCAATCTGCCGATAGGGCGGTGTTGTTTGGGCAC‑AGTGTGCt                                                                                                                            >  2:63999/1‑251 (MQ=255)                                                                                                                                  tctatCATGGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGAGTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATT‑GAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGAGTGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCtggtg                                                                                                                        <  2:470304/240‑1 (MQ=255)                                                                                                                                  gCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTAGCATGGATAATCGCCAGTCTGATTGGGGCGGTAATGAACGGTATTTATTTGCAAGCGGTATT‑GAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGATTTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGGGTGCtgtggggatcgg                                                                                                                 >  1:470304/1‑241 (MQ=255)                                                                                                                                       cATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGATCGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGcc                                                                                                            <  1:58897/251‑1 (MQ=255)                                                                                                                                       cATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGc                                                                                                             <  1:404040/250‑1 (MQ=255)                                                                                                                                                    ttATTTGCTGTTTCTGGCGCGTTCTGAACATATCGGTTTTACCGTTGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGttt                                                                                                <  2:38320/250‑1 (MQ=255)                                                                                                                                                                                 atatcGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGttt                                                                  <  1:314011/247‑1 (MQ=255)                                                                                                                                                                                     cggtttttcttgtgCATGGCTAATCGTCTGTCTGTTTGGTGCGTTAATGAACGGTATTTATTTGCAGGCGGTAGT‑GAAAGGTTGGCGTAACAGCTTGTGGTTGTCCCTCGCGCTGTTTTTGCTGGATGTTGTGATGTGGGGACTGCTCATCTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCTGGCATGATGTTTGTGACTCGTAATATCGACTTGTATGCGGTTTCa                                                               <  2:476684/237‑1 (MQ=255)                                                                                                                                                                                                    aTGGCTAATCGCAAGGCTGAGTGGGGATTTAATGATCGGTTTTTATTTGCAGGCGGTATT‑GAAAGGCTGGCTCAACATCATGTTGTTTTCCCTCTCGCTGGTGTTGCTGGAGGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGTCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAg                                               <  2:406393/251‑1 (MQ=255)                                                                                                                                                                                                        aTAACCGCCAGACTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGAGGGTGTGATGTGGGGACTGCTCAACGCTGcc                                                                                                                                                             <  2:64429/138‑1 (MQ=255)                                                                                                                                                                                                                cAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGGTGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCCGTGTGCTGGTGGTGGCGCTGGCCGGCATGGTGTTTTTGACTCGTAATATCGACTGGGATGCGTTTTCACTGCCGAAAATGAAAACCAGTaaaaaa                                    >  2:185456/1‑247 (MQ=255)                                                                                                                                                                                                                          ggggCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTAACCTCGCGCTGTTGTTGCTGGGTGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACaaaggg                         >  1:371056/1‑247 (MQ=255)                                                                                                                                                                                                                                          ggatttatttGCAGTTGGGATT‑GAAAGGCGGGCGCAACATCATGTGGTTTACCCTCTCGCTTTTGGTGCTGGATGGTGTGAGTTGGGGACTGCTGAACTCTGCCGATAGCGCGCTGTTGTG‑GGGAACCATTTTGCTGGTGGTGGCGCTGGCCGGCATGATGGTTGTGACTGGTAAGATCGACTTGTATGCGTTTGTACGGCCGAAAATGAAAGCCAGTAAAGAAGTTACaatggtcgcgcaattacgtatc         <  2:22394/249‑21 (MQ=255)                                                                                                                                                                                                                                          gTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGGGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCCTGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTAAGTATc         >  1:390815/1‑251 (MQ=255)                                                                                                                                                                                                                                               tatttGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGa     >  1:234528/1‑250 (MQ=255)                                                                                                                                                                                                                                                 tttGCAAGCGGTATT‑GAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAAt  >  1:246786/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       TCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATT‑GAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAAT  >  NC_000913/4639038‑4639528

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