Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,217	C→T	T347I (ACC→ATC)	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,217	0	C	T	93.6%	127.5 / ‑1.9	47	T347I (ACC→ATC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (22/22);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATT‑GAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTT  >  NC_000913/4638979‑4639434                                                                                                                                                                                                                                               |                                                                                                                                                                                                                            aaTGTTACCGGATGGGAAACCTCCCCCGCGTTTAGCGGCGCAGTGACGTCCCCAGCCGATCAATACCTATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGt                                                                                                                                                                                                                 <  2:382535/251‑1 (MQ=255)       aCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGGGGGGCTTTCATTGGTGGTGTTTTATTTGCTCTTGATGGAGCTTTCTGAACATATCGGTTTTACCGTGGCATg                                                                                                                                                                                                           >  1:16395/1‑251 (MQ=255)         cGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATgg                                                                                                                                                                                                          >  1:416810/1‑250 (MQ=255)                gAAAACTTCCCCGCGTTTAGCGGCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGc                                                                                                                                                                                                  <  1:83155/251‑1 (MQ=255)                  aaaCTTCCCAGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCa                                                                                                                                                                                                <  2:33660/251‑1 (MQ=255)                                ttAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCTCCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGTTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGc                                                                                                                                                                                   <  1:128913/250‑1 (MQ=255)                                     gtggcagggacgacgCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCAGGGATAATCGCCAGTCTGATTGTGGCGTTAAt                                                                                                                                                                             <  2:211676/242‑1 (MQ=255)                                       cGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATg                                                                                                                                                                            <  1:251357/250‑1 (MQ=255)                                                 acgCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGCTAATCGCCAGTCAGATTGGGGCGTTACTGAACGGtattt                                                                                                                                                                  >  1:57454/1‑250 (MQ=255)                                                         cGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCTGTCTGATTGGGGCGTTACTGTACGGAatttattttcagg                                                                                                                                                         >  2:38211/1‑246 (MQ=255)                                                          gATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCaggc                                                                                                                                                        >  1:205389/1‑248 (MQ=255)                                                            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)                                                                                                         gATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCa                                                                                                                                                                                                <  1:128206/164‑1 (MQ=255)                                                                                                         gATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCa                                                                                                                                                                                                >  2:128206/1‑164 (MQ=255)                                                                                                              cACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGGTGGAGCAACAGACTGTTGTTTACCCTcgcgatgactgtg                                                                                                   >  1:90736/1‑242 (MQ=255)                                                                                                               aCTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATT‑GAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGc                                                                                                  >  1:323368/1‑251 (MQ=255)                                                                                                                    cTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTAcc                                                                                                                                                                                                                   >  2:338594/1‑134 (MQ=255)                                                                                                                    cTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTAcc                                                                                                                                                                                                                   <  1:338594/134‑1 (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)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                            AATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATT‑GAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTT  >  NC_000913/4638979‑4639434

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