Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,771,979	G→T	intergenic (‑235/‑302)	yibI ← / → mtlA	DUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,771,979	0	G	T	97.5%	118.4 / ‑3.7	40	intergenic (‑235/‑302)	yibI/mtlA	DUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (21/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

ACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTG  >  NC_000913/3771741‑3772212                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           aCATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCCTAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAAcc                                                                                                                                                                                                                                >  2:193982/1‑251 (MQ=255)        tGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATcaaa                                                                                                                                                                                                                          >  1:451318/1‑250 (MQ=255)            gcGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATcaaaacaaa                                                                                                                                                                                                                     <  1:44738/250‑1 (MQ=255)                  tCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGCTCACAAATATTCAACAAACCATCAAAACAAAAAtgtg                                                                                                                                                                                                               >  1:243317/1‑251 (MQ=255)                          aTCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACactact                                                                                                                                                                                                       >  1:273344/1‑251 (MQ=255)                                gttgttTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGATGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAAAAAAAATGTGACACTCCTCACAtt                                                                                                                                                                                                 >  2:355832/1‑251 (MQ=255)                                   gtgtttGGGCTGGCTGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCGCGCCATAACACCAAAGAATAATTTTTATAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGAGCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTaa                                                                                                                                                                                              <  1:31455/248‑1 (MQ=255)                                       tggggCTGGCAGCCATAAGGGAGTTAGGATGATATTTTGACAATAATCCGGGTTCGCTATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAATAAACCATCAAAACAAAAATGTTACACTACTCACATTTAAATg                                                                                                                                                                                           <  2:220788/248‑1 (MQ=255)                                              ggCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAtttt                                                                                                                                                                                    <  1:4970/250‑1 (MQ=255)                                                 aGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGCTTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAACTGTGACAGTCCTCACATTTAAATGCCATTTTTAGc                                                                                                                                                                                >  2:312982/1‑251 (MQ=255)                                                 aGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGc                                                                                                                                                                                >  1:363987/1‑251 (MQ=255)                                                            aGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTACGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATcgccg                                                                                                                                                                     >  1:210505/1‑251 (MQ=255)                                                                     gATATTTTGACAATAATTCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAAGTTGTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCttgtt                                                                                                                                                              <  1:100267/249‑1 (MQ=255)                                                                       aatTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCtt                                                                                                                                                          <  1:139993/250‑1 (MQ=255)                                                                        atTTTGAAAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACTCTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCtt                                                                                                                                                          <  2:464237/250‑1 (MQ=255)                                                                          tttGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCttttt                                                                                                                                                       >  1:396189/1‑251 (MQ=255)                                                                              aCAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTac                                                                                                                                                    >  1:333093/1‑250 (MQ=255)                                                                                 ataatCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAATTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGCTCACAACTATTCAACAAACAATCAAAACAAAAATGTGCCACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTAcaca                                                                                                                                                 >  1:429095/1‑250 (MQ=255)                                                                                             cgcgATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGAGCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCATATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATTCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑Gtgat                                                                                                                                   <  1:310969/251‑1 (MQ=255)                                                                                                 aTTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAac                                                                                                                               <  1:269432/251‑1 (MQ=255)                                                                                                      cGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGGTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGAAACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTAACACAA‑CGTTTTTGTGATGAacgtccc                                                                                                                          >  2:81434/1‑249 (MQ=255)                                                                                                                         aatTTTTAGAGGTGATGAGTTGCTTAGTTACATACCGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATAACAAATATTCAACAAACCATCAAAACACAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTAccc                                                                                                       >  2:209537/1‑251 (MQ=255)                                                                                                                          atTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTAccac                                                                                                      >  1:80590/1‑249 (MQ=255)                                                                                                                             tttAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCt                                                                                                    <  2:451318/250‑1 (MQ=255)                                                                                                                                         gattggctaaGTTACAAAACGAGTGTATGAAGAAGGCAAAACATGCTGTAGAACACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATGCAGATCACAAATATTCAACAAACCTTCAAAGCAAAAATGTGACACTACTCTCATTTAAATGCCATTTTTAGCGAAAATCGCCGACGTGTTGCTTTTTGTCACAAGCGTTTT‑GTGATGAACGTCACGTCAATTAACTCTCTCCCCCCTATATTTAtgtgtt                                                                                       <  2:281987/242‑3 (MQ=255)                                                                                                                                          aattGCTTATTTACATTACTATTGTTTTACGAATGCATAACATGCTGGAGATCACATCTGGTGAACGCCGTAAGAAAATATCTTGTTATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACTCTACTCTCATTTAAATGCCATGTTTAGCGAAAATCGCCGCCTTTTTGCTTTTTTACACAAGCGTTTT‑TTTATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGtgatt                                                                                      <  1:29901/249‑1 (MQ=255)                                                                                                                                                        aTAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACaaa                                                                         >  2:141800/1‑250 (MQ=255)                                                                                                                                                              aTTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGt                                                                  <  2:234139/251‑1 (MQ=255)                                                                                                                                                                tGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGATACTACTCTCATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTCCCCCCtata                                                                                                 <  2:59171/218‑1 (MQ=255)                                                                                                                                                                tGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCtata                                                                                                 >  1:59171/1‑218 (MQ=255)                                                                                                                                                                         aaGGCATAACATGCTGTAGATCACATCAAGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACACACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTCTATTTATGTGATTGATATCACACAAAAGGCCGTCGACTgtac                                                         >  1:400534/1‑246 (MQ=255)                                                                                                                                                                                aaCATGCTGAAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACACATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGACTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCCCGTACATTACCTCTCTACCCCCTATATTTTTGTGATTGTTATCACACAAAAGGCCGTCGACTGAACAGTAAAc                                                  >  2:202568/1‑249 (MQ=255)                                                                                                                                                                                          aGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGcc                                      <  1:114983/251‑1 (MQ=255)                                                                                                                                                                                           gATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGcc                                      >  2:457142/1‑250 (MQ=255)                                                                                                                                                                                                caTCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGAttt                                >  2:259200/1‑251 (MQ=255)                                                                                                                                                                                                   cAGGTGAACGACGTAAGAAAATATCTTGTGATACAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATTCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCCACCCCCTATATTTATGTGAATGATACCACACAAAACGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCg                              <  1:455227/250‑1 (MQ=255)                                                                                                                                                                                                    aGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCa                            <  1:193982/251‑1 (MQ=255)                                                                                                                                                                                                       tGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGCTTTCGCAGt                          >  1:434682/1‑250 (MQ=255)                                                                                                                                                                                                                              gtgATTCATATCACAAATATACAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCGCTACCCCCTATATTTATGTGAGTGATATCACACAAAAGGCCGTCGACTTGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCCACAAGGGCCGGCTCCCTCTg  <  2:29097/251‑1 (MQ=255)                                                                                                                                                                                                                              gtgATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTg  <  2:333093/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           ACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTT‑GTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTG  >  NC_000913/3771741‑3772212

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