Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,676,106	0	C	A	43.2%	118.8 / 37.9	74	L240F (TTG→TTT)	yphD	putative sugar ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (32/10);  new base A (0/32);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 1.71e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAATCACCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAAC  >  NC_000913/2675959‑2676256                                                                                                                                                     |                                                                                                                                                          cGTAATCACCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCcaacaa                                                                                                                                                         <  3:41050/150‑1 (MQ=255) cGTAATCACCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACaaa                                                                                                                                                        <  1:1133712/151‑1 (MQ=255)  gTAATCACCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACaaa                                                                                                                                                        >  1:266062/1‑150 (MQ=255)    aaTCACCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑g                                                                                                                                                    >  2:80635/1‑151 (MQ=255)     aTCACCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑Gc                                                                                                                                                   <  1:782476/151‑1 (MQ=255)       cACCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑g                                                                                                                                                >  4:277929/1‑151 (MQ=255)       cACCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑g                                                                                                                                                >  4:241754/1‑151 (MQ=255)        aCCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACTGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑gg                                                                                                                                               >  2:692381/1‑151 (MQ=255)          cAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑c                                                                                                                                            >  4:112562/1‑151 (MQ=255)           aGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑c                                                                                                                                            >  4:318815/1‑150 (MQ=255)              acacCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CAccg                                                                                                                                        <  1:318871/151‑1 (MQ=255)                 ccAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CAccgcc                                                                                                                                      <  3:20833/150‑1 (MQ=255)                       aatgccccaaaCAGGAAACCGCGCCCCCCGAAAAGTGCCGTCCCCCCGACCCCGCCCGCGGCGATGACGTCAAACCCCAGCCCTTTTGCGCCCCCGGCTTTCCCAGACCC‑GAGGCGCCCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCtaat                                                                                                                                <  3:244377/141‑1 (MQ=255)                           tACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATcc                                                                                                                           >  2:678281/1‑151 (MQ=255)                                aaCAGGAAACCGCCCCCGCCGAAAAGTCCCGTCCCCCCGCCCCCGCCCGCGGCGATGACGCCAACCTCCAGCCCGTTGGCGGCACCTGCGTTCCCAGAACC‑GAGGCGCGCCGCCAAAAAAAG‑GCG‑GTT‑CACCGCCGCTAATAATCCCGaaa                                                                                                                      <  1:508860/151‑1 (MQ=255)                                  cAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAgg                                                                                                                    >  1:80055/1‑151 (MQ=255)                                    ggcacccGCGCCCGCCGAAAAGCGCCGTCCCGCCGCCCACCCCCGCGGCGATGACGACAAACCCCAGACCTTTGGCGCCACCTGCGTTACCAAACCC‑GAGGGGCCCCCCAAAAAAACC‑GCC‑GTT‑CCCCGCCGCTAATAATCCCGAAAgg                                                                                                                    <  2:324839/144‑1 (MQ=255)                                    ggaccccgcccccgcggGAAAGTGCCGTCCCGCCGCCCGCGCCCGGGGCGATGACGACAAGCCCCACCCCGTTGTCGGCACCTGCGTTCACAGACCC‑GGGGCGCGCCGCCAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAg                                                                                                                     <  2:893730/137‑1 (MQ=255)                                     ggAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGt                                                                                                                  >  2:832203/1‑150 (MQ=255)                                       aaCCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTaa                                                                                                                >  4:308341/1‑150 (MQ=255)                                         ccGCGCCCGCCGAAAAGTGCCGTCCCGCCGCCCACGCCCGCGGCGATGACGTCAACCTCCAGACCGTTTGCGGCCCCTGCGTTCCCAGAACC‑GAGGCGCGCCGCCAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTaa                                                                                                                <  1:149608/148‑1 (MQ=255)                                         ccGCGCCCGCCGAAAAGTGCCGTCCCGCCGACCCCGCCCGCGGCGATGCCGTCAACCCCCAGACCTTTTGCGGCCCCTGCGTTCCCAGAACC‑GAGGCGCGCCGCCAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGa                                                                                                             <  4:285931/151‑1 (MQ=255)                                           gcgACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGAt                                                                                                            >  2:464382/1‑150 (MQ=255)                                            cgACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATa                                                                                                           >  4:154606/1‑150 (MQ=255)                                               ccgcggGAAAGTGCCGTCCCGCCGCCCACGCCCGGGGCAATGACGTCAACCCCCAGCCCTTTTGCGCCCCCTGCTTTCCCAGACCC‑GAGGCGCGCCGCCAAAAAAAG‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATa                                                                                                           <  1:1074537/147‑1 (MQ=255)                                               ccgccgGAAAGCCCCGACCCCCCGCCCACGCCCGGGGCGATGACGACAACCCCCAGACCGTTTGCGGCCCCTGCGTTCCCAAAACC‑GGGGCGCGCCGCCAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGa                                                                                                        <  2:838511/150‑1 (MQ=255)                                                  ccgGAAAGGCCGGTCCCCCCGACCACGCCCGCGGGGGTGACGCCAACCCCCAGCCCGTTTGGGGCCCCGGCTTTCCCAAAACC‑GAGGCGGGCCGCAAACAAAAG‑GCG‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATgc                                                                                                    <  3:156453/151‑1 (MQ=255)                                                      aaaGGGCGGAACCGCCACCCACGCCGGCGGAGAAGACGACAAACCCAAGACCGTTGGCGGCACTTGCGATGAAAGAAGG‑GAGGGGCGCGACAAAAAAAAG‑GCG‑GGC‑CACCGCCGCTAATAATCCCGAAAGGGTaac                                                                                                               <  2:853715/136‑2 (MQ=255)                                                           gCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCcgac                                                                                           >  2:1123668/1‑151 (MQ=255)                                                            ccGTCCCGCCGCCCCCGCCCGCGGCGATGCCGTAAAACCCCAGCCCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCcgacg                                                                                          <  1:843075/151‑1 (MQ=255)                                                               cccccccgcccccgcccGGGGCGAGGCCTTCAACCCCCAGCCCTTTTGCGGCCCCGGCTTCCCCAAAACC‑GGGGCGCCCCGCAAAAAAAAT‑GCC‑GCC‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAAc                                                                                       <  2:293330/146‑1 (MQ=255)                                                                ccccccgACCCCGCCCGCGGCGATGACGTAAACCCCCAGCCCTTTGGCGGCCCCTGCGTTCCCAAACCG‑GGGGCGCGCCGCCAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACg                                                                                      <  3:235023/147‑1 (MQ=255)                                                                 ccgccgCCCCCGACCGCGGCGATGACGTCAAACCCCAGACCGTTTGCGGCACCTGCGTTCCCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAAc                                                                                       <  1:317582/149‑1 (MQ=255)                                                                 ccgccgACCCCGCCCGCGGCGATGACGACAAACCCCAGACCTTTTGCGGCCCCTGCGTTACCAGACCC‑GAGGCGCGCCGCCAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGt                                                                                     <  2:687257/151‑1 (MQ=255)                                                                 ccgccgACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCCCCTGCGTTCCCAGAACC‑GAGGCGCGCCGCAAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAAc                                                                                       <  1:431034/149‑1 (MQ=255)                                                                 ccgccgACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGt                                                                                     >  3:143519/1‑151 (MQ=255)                                                                 cccccgCCCACGACCGCGGCGATGACGTCAAACCCCAGCCCTTTGGCGGCCCCGGCGTTACCAGACCC‑GGGGCGCGCCGCCAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGaa                                                                                        <  1:350113/145‑1 (MQ=255)                                                                  cgccgACCACGCCCGCGGCGATGCCGTAAAACTCCAGCCCGTTTGCGGCACCGGCGTTACCAGACCC‑GGGGCGCGCCGCAAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAAc                                                                                       <  4:207130/148‑1 (MQ=255)                                                                    ccgACCACGACCGCGGCGATGACGTCAACCTCCAGACCGTTTGCGGCACCTGCGTTACCAGACCC‑GAGGCGCGCCGCCAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGa                                                                                  <  4:104773/151‑1 (MQ=255)                                                                       aCCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGc                                                                               >  4:419400/1‑151 (MQ=255)                                                                        ccccgcccGGGGCGATGACGTAAAACCCCAGACCGTTTGCGGCCCCTGCTTTCCCAGAACC‑GAGGGGCGCCCCAAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAGGATAAGAATGCGCACCCGACGAACGtt                                                                                    <  4:298934/139‑1 (MQ=255)                                                                              ccGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGcaca                                                                         >  3:63219/1‑150 (MQ=255)                                                                              ccGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGcaca                                                                         >  1:960337/1‑150 (MQ=255)                                                                                  ggcgATGACGTCAACCCCCAGACCGTTTGCGGCACCTGCGTTACCAGACCC‑GAGGCGCGCCGCCAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTgc                                                                    <  2:72766/151‑1 (MQ=255)                                                                                  ggcgATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTgc                                                                    >  3:378388/1‑151 (MQ=255)                                                                                     gATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGcgcc                                                                 >  3:292469/1‑151 (MQ=255)                                                                                         aCGTCAACCCCCAGCCCGTTGGCGGCCCCGGCGTTACCAGACCC‑GAGGCGCGCC‑CCAAAAAAAA‑GCCGGGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGccgccg                                                             <  1:901046/151‑1 (MQ=255)                                                                                          cGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGt                                                            >  4:255489/1‑151 (MQ=255)                                                                                                 ccccACACGGTTGGCGCCCCCTGCGTTCCAAGACCG‑GGGGCGCGCCCCAAAAAAAAT‑GCG‑GGC‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCAtt                                                      <  3:349686/148‑1 (MQ=255)                                                                                                  ccaagccggTTTGCGCCCCCTGCGTTACCACACCG‑GAGGCGCGCCGCAAAAAAAAG‑CCC‑GGC‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCAtt                                                      <  2:1275466/141‑1 (MQ=255)                                                                                                   ccAGCCCTTTTGCGGCACCTGCGTTACCAAAACC‑GAGGCGCCCCCCAAACAAATC‑GCC‑GGT‑CCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTAc                                                    <  1:765575/150‑1 (MQ=255)                                                                                                       aCCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGcc                                                >  3:245527/1‑150 (MQ=255)                                                                                                        ccGTTGGCGGCCCCTGCTTTACAAGACCG‑GGGGCGCGCCCCAAAAAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCaa                                              <  4:441730/151‑1 (MQ=255)                                                                                                          gTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACg                                            >  2:750758/1‑151 (MQ=255)                                                                                                                     ccggcttACCAAAACCGGGGGCCCGCC‑CCAAAAAAATGGCG‑GGC‑CCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCaa                                        <  2:271494/140‑1 (MQ=255)                                                                                                                      tGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCg                                  >  1:899889/1‑149 (MQ=255)                                                                                                                        cGTTACCAAAGCC‑GAGGCCCCCCGCCAAAAAAAT‑GCC‑GAT‑AACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGa                                 <  1:760237/148‑1 (MQ=255)                                                                                                                         gTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAgcgc                             >  2:214169/1‑151 (MQ=255)                                                                                                                          ttccaaGACCC‑GGGGCGCCCCCCAAAAAAAAG‑CCC‑GTT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAgcg                              <  1:31742/144‑1 (MQ=255)                                                                                                                          ttACCAGACCG‑GGGGCGCGCCCCAAAAAAAAT‑CCG‑GGT‑CCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAgcg                              <  1:773675/149‑1 (MQ=255)                                                                                                                            aCCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCg                          >  4:75745/1‑151 (MQ=255)                                                                                                                              caaaaccgggggCGCGCCCCAAAAAAAAG‑GCC‑GGT‑CCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAgcgc                             <  1:114795/137‑1 (MQ=255)                                                                                                                               aaaccc‑GGGGCGCGCGGCCAAAAAAAT‑GCC‑GGT‑CCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGa                         <  3:277925/145‑1 (MQ=255)                                                                                                                               aGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGa                         >  1:290714/1‑149 (MQ=255)                                                                                                                                 accc‑GGGGCGCCCCGCAAAAAAAAT‑CCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAgg                      <  1:725308/148‑1 (MQ=255)                                                                                                                                   cc‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCgg                   >  4:354328/1‑151 (MQ=255)                                                                                                                                       gggcccccccccAAAAAAAA‑GCC‑CGTCCCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCCCCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAgg                      <  3:320226/136‑1 (MQ=255)                                                                                                                                        ggCGCCCCCCAAAAAAAAT‑GCC‑GTT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGtttt               <  4:426341/151‑1 (MQ=255)                                                                                                                                          cgcgcCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGc             >  2:448615/1‑151 (MQ=255)                                                                                                                                          cgcccCGCAAAAAAAAC‑GCC‑GCT‑CCCCGCCGCTAATAATCCCGAAAGGTTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTg              <  2:368615/150‑1 (MQ=255)                                                                                                                                           gcgcCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAGCGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCg            >  1:195999/1‑151 (MQ=255)                                                                                                                                            ccccgcaAAAAAAAG‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAATGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCgg           <  1:898623/149‑1 (MQ=255)                                                                                                                                             gccgccAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGc          >  2:570621/1‑151 (MQ=255)                                                                                                                                              ccgccAAAAAAAG‑GCC‑GGT‑CCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCt         <  1:559094/151‑1 (MQ=255)                                                                                                                                                 ccAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCtgat      >  4:155921/1‑151 (MQ=255)                                                                                                                                                    acaaACC‑GCC‑GAT‑CCCCGCCGCTAATAACCCCGAAAGGTTAAAGATAAGAAGGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCtgatga    <  2:668063/150‑1 (MQ=255)                                                                                                                                                     caaAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAac  >  2:284792/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                          CGTAATCACCAGCACACCAAGCAATGTACCGAACAGGGAACCGCGACCGCCGGAAAGTGCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACC‑GAGGCGCGCCGCCAACAAAAT‑GCC‑GGT‑CACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAAC  >  NC_000913/2675959‑2676256

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