Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,303	T→C	C376R (TGC→CGC)	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,303	0	T	C	88.2%	58.6 / 1.8	36	C376R (TGC→CGC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (11/19);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

GACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAA  >  NC_000913/4639054‑4639538                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             gACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTg                                                                                                                                                                                                                                             <  2:174892/251‑1 (MQ=255)                tACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGGGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTCTTTGCAGGCGGGATTAAAAGGGTGGAGAAACCGCGTgttgtt                                                                                                                                                                                                                              >  2:441974/1‑251 (MQ=255)                  cGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTa                                                                                                                                                                                                                            <  2:232684/251‑1 (MQ=255)                        tCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTc                                                                                                                                                                                                                       <  2:98773/250‑1 (MQ=255)                           gctgATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTcg                                                                                                                                                                                                                      <  2:331507/248‑1 (MQ=255)                            ctgATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCCCGCt                                                                                                                                                                                                                  >  1:208651/1‑251 (MQ=255)                             tgATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTTAAAGGGTGGCGCAACCGCATGTTGttttaatcgggaggg                                                                                                                                                                                                                 >  2:6983/1‑239 (MQ=255)                             tgATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGttt                                                                                                                                                                                                                             <  1:6983/239‑1 (MQ=255)                                   cACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGGTGGAGCAACAGACTGTTGTTTACCCTcgcgatgactgtg                                                                                                                                                                                                           >  1:90736/1‑242 (MQ=255)                                    aCTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGc                                                                                                                                                                                                          >  1:323368/1‑251 (MQ=255)                                             tATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTTCCCTCGCGCAGTTGCTGCTgggatggg                                                                                                                                                                                                 >  1:474480/1‑245 (MQ=255)                                                          gTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACt                                                                                                                                                                                    <  1:61603/251‑1 (MQ=255)                                                                  aaCGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGGGGGGCTTTCATTTGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAAct                                                                                                                                                                            <  2:495646/251‑1 (MQ=255)                                                                                   gTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTAACCTCGCGCTGTTGTTGCTGGATGGTGGGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCtg                                                                                                                                                           >  1:406261/1‑251 (MQ=255)                                                                                     ttACACCCAATGCAATTTTTGCTGGTGGGGCTTTCATTGGTGATGTTTTAGTTGCTCTTGCTGGCGTTTTCTTAACATATCGTTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCtgt                                                                                                                                                          <  2:323368/250‑1 (MQ=255)                                                                                         ccccAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGttgtt‑g                                                                                                                                                    <  2:205389/250‑1 (MQ=255)                                                                                                 cAATATTTGCGGGTGGGTCTGTCATTGGTGATGTTTTAGTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTCCCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGCT‑GGGAACCAg                                                                                                                                            <  2:16395/251‑1 (MQ=255)                                                                                                   atatTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGTTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAgtg                                                                                                                                          <  1:264511/251‑1 (MQ=255)                                                                                                       ttGCTGGGGGGGCTTTCGTTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACTGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCAAGTTGTTTTACCTCGCGGGGGTGTTGATGGATGGAGGGATGTGGGGACTGGTCCAATCTGCCGATAGGGCGGTGTTGTTTGGGCAC‑AGTGTGCt                                                                                                                                      >  2:63999/1‑251 (MQ=255)                                                                                                                  tctatCATGGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGAGTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGAGTGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCtggtg                                                                                                                                  <  2:470304/240‑1 (MQ=255)                                                                                                                  gCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTAGCATGGATAATCGCCAGTCTGATTGGGGCGGTAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGATTTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGGGTGCtgtggggatcgg                                                                                                                           >  1:470304/1‑241 (MQ=255)                                                                                                                       cATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGATCGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGcc                                                                                                                      <  1:58897/251‑1 (MQ=255)                                                                                                                       cATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGc                                                                                                                       <  1:404040/250‑1 (MQ=255)                                                                                                                                    ttATTTGCTGTTTCTGGCGCGTTCTGAACATATCGGTTTTACCGTTGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGttt                                                                                                          <  2:38320/250‑1 (MQ=255)                                                                                                                                                                 atatcGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGttt                                                                            <  1:314011/247‑1 (MQ=255)                                                                                                                                                                     cggtttttcttgtgCATGGCTAATCGTCTGTCTGTTTGGTGCGTTAATGAACGGTATTTATTTGCAGGCGGTAGTGAAAGGTTGGCGTAACAGCTTGTGGTTGTCCCTCGCGCTGTTTTTGCTGGATGTTGTGATGTGGGGACTGCTCATCTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCTGGCATGATGTTTGTGACTCGTAATATCGACTTGTATGCGGTTTCa                                                                         <  2:476684/237‑1 (MQ=255)                                                                                                                                                                                    aTGGCTAATCGCAAGGCTGAGTGGGGATTTAATGATCGGTTTTTATTTGCAGGCGGTATTGAAAGGCTGGCTCAACATCATGTTGTTTTCCCTCTCGCTGGTGTTGCTGGAGGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGTCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAg                                                         <  2:406393/251‑1 (MQ=255)                                                                                                                                                                                        aTAACCGCCAGACTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGAGGGTGTGATGTGGGGACTGCTCAACGCTGcc                                                                                                                                                                       <  2:64429/138‑1 (MQ=255)                                                                                                                                                                                                cAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGGTGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCCGTGTGCTGGTGGTGGCGCTGGCCGGCATGGTGTTTTTGACTCGTAATATCGACTGGGATGCGTTTTCACTGCCGAAAATGAAAACCAGTaaaaaa                                              >  2:185456/1‑247 (MQ=255)                                                                                                                                                                                                          ggggCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTAACCTCGCGCTGTTGTTGCTGGGTGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACaaaggg                                   >  1:371056/1‑247 (MQ=255)                                                                                                                                                                                                                          ggatttatttGCAGTTGGGATTGAAAGGCGGGCGCAACATCATGTGGTTTACCCTCTCGCTTTTGGTGCTGGATGGTGTGAGTTGGGGACTGCTGAACTCTGCCGATAGCGCGCTGTTGTG‑GGGAACCATTTTGCTGGTGGTGGCGCTGGCCGGCATGATGGTTGTGACTGGTAAGATCGACTTGTATGCGTTTGTACGGCCGAAAATGAAAGCCAGTAAAGAAGTTACaatggtcgcgcaattacgtatc                   <  2:22394/249‑21 (MQ=255)                                                                                                                                                                                                                          gTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGGGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCCTGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTAAGTATc                   >  1:390815/1‑251 (MQ=255)                                                                                                                                                                                                                               tatttGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGa               >  1:234528/1‑250 (MQ=255)                                                                                                                                                                                                                                 tttGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAAt            >  1:246786/1‑251 (MQ=255)                                                                                                                                                                                                                                    gcgtgtggTATTGAAAGGCGGGCGCAACAGCATGTTGTGATCCCTCGGGCTGTTGGTGCTGGAGGGTGTTATGTTGGGCCTGCTCAACTCGGCCGATAGCGAGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATAAAGACTGGTATGCGTTTTCACTGCCGAATATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATaa          <  2:163111/244‑1 (MQ=255)                                                                                                                                                                                                                                            tATTGACAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGGTCAATTACGTATCTGGAAATAAGGTTGaaa  >  2:114320/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             GACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAA  >  NC_000913/4639054‑4639538

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: ‑