Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,261	A→T	I362L (ATA→TTA)	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,261	0	A	T	95.2%	124.2 / ‑3.1	43	I362L (ATA→TTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (20/20);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATT‑GAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGA  >  NC_000913/4639015‑4639505                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       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)                                                                          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)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATT‑GAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTT‑GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGA  >  NC_000913/4639015‑4639505

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