Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,360	C→T	8.8%	intergenic (‑193/+54)	fbp ← / ← glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,743,360	0	C	T	8.8%	153.7 / 4.7	57	intergenic (‑193/+54)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (23/28);  new base T (2/3);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

TACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTT  >  NC_002947/5743212‑5743506                                                                                                                                                       |                                                                                                                                                   tACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGcc                                                                                                                                                    >  6:255937/1‑149 (MQ=255)   cAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                                                  <  5:202621/149‑1 (MQ=255)   cAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                                                  <  8:280071/149‑1 (MQ=255)         gagTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCt                                                                                                                                            >  2:84479/1‑149 (MQ=255)            tGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                         <  8:176746/149‑1 (MQ=255)              aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                       <  4:56372/149‑1 (MQ=255)              aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                       <  8:222923/149‑1 (MQ=255)              aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                       <  3:192747/149‑1 (MQ=255)              aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCg                                                                                                                                        <  7:250648/148‑1 (MQ=255)                gggCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCgaga                                                                                                                                     <  1:294796/149‑1 (MQ=255)                    tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtacc                                                                                                                                 >  7:112013/1‑142 (MQ=25)                    tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                                                 >  4:247290/1‑133 (MQ=255)                    tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                                                 <  3:247290/133‑1 (MQ=255)                    tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTGCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGCCCATTCGCTCGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                                                  <  1:84479/149‑1 (MQ=255)                        gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                                             <  5:255937/149‑1 (MQ=255)                        gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGT‑GTCTGGTTGTGTGTCCATTCGC‑GGGCACGCCC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                            <  2:127323/149‑1 (MQ=18)                             atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGc                                                                                                                        >  3:236987/1‑149 (MQ=255)                              tgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGATCTGGTTGTGTGTCCATTCGC‑GGGCACGCCC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                            <  1:189265/144‑1 (MQ=18)                              tgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGATCTGGTTGTGTGTCCATTCGC‑GGGCACGCCC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                            >  2:189265/1‑144 (MQ=18)                                agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCCGTGGGAGCga                                                                                                                      <  8:112013/148‑2 (MQ=21)                                  aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCg                                                                                                                   >  7:203493/1‑149 (MQ=255)                                       ccACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGccc                                                                                                              >  3:257986/1‑149 (MQ=255)                                         aCGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGGCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcgc                                                                                                            <  7:301037/149‑1 (MQ=255)                                            gCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGAGAGACACCTGTGGGAGCGGGCGTACCcgaac                                                                                                          >  6:92021/1‑145 (MQ=255)                                            gCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCAATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTACCcgaac                                                                                                          >  6:96796/1‑145 (MQ=255)                                              tctcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTAGAACCTTGCGAGACACCTGTGGGAGCGCGCGTGCCCGCGAACg                                                                                                       >  6:264262/1‑149 (MQ=255)                                                  aGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcct                                                                                                   >  7:299883/1‑149 (MQ=255)                                                    aCGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctca                                                                                                 >  3:147991/1‑149 (MQ=255)                                                                     tAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCtt                                                                                <  4:129546/149‑1 (MQ=255)                                                                        ttCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTCGGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTc                                                                               >  1:93614/1‑148 (MQ=255)                                                                         tCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCa                                                                             <  8:299883/148‑1 (MQ=255)                                                                               ttCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCg                                                                      >  7:161828/1‑149 (MQ=255)                                                                                  gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                                   <  2:93614/149‑1 (MQ=255)                                                                                       gCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTAGTTGGCata                                                              >  6:314382/1‑149 (MQ=255)                                                                                            tggtCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCAGAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGACATATTcgc                                                         <  4:236987/149‑1 (MQ=255)                                                                                                 tGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCACCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATg                                                    >  1:35484/1‑149 (MQ=255)                                                                                                   gttgttTGTCCTTTGGG‑TGGCGCGGTC‑GCTCCCACAGGTACAGCTCTAGCATTGAAGTTTGTGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATgg                                                   <  1:134678/148‑1 (MQ=255)                                                                                                    ttGTGTCGCCCTTCGC‑GGGCGCGCTC‑GCTCCCACCGGTACAGCCCAAGCCTTGAAGCTTGCGAGACACCTGTTGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                                  <  4:257986/148‑1 (MQ=38)                                                                                                                 cgc‑gGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACgag                                     <  2:35484/148‑1 (MQ=255)                                                                                                                  gc‑gGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACgaga                                    <  8:203493/148‑1 (MQ=255)                                                                                                                        cgcggTC‑GCTCCCAGAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCGGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGc                              <  5:264262/149‑1 (MQ=255)                                                                                                                           gctc‑gctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCa                           >  3:27610/1‑149 (MQ=255)                                                                                                                             tc‑gctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATc                         >  6:52934/1‑149 (MQ=255)                                                                                                                              c‑gctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc                        <  5:1895/149‑1 (MQ=255)                                                                                                                                gctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc                        <  5:52934/148‑1 (MQ=255)                                                                                                                                  tcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtg                     >  5:71571/1‑149 (MQ=255)                                                                                                                                    ccACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGc                              >  4:242290/1‑138 (MQ=255)                                                                                                                                    ccACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGc                              <  3:242290/138‑1 (MQ=255)                                                                                                                                    ccACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc                   >  6:93694/1‑149 (MQ=255)                                                                                                                                    ccACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc                   >  6:88914/1‑149 (MQ=255)                                                                                                                                     cacaGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                 <  2:66632/118‑1 (MQ=255)                                                                                                                                        aGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCgtgt               <  4:147991/149‑1 (MQ=255)                                                                                                                                         ggTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt              <  6:71571/149‑1 (MQ=255)                                                                                                                                         ggTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt              <  4:27610/149‑1 (MQ=255)                                                                                                                                          gTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt              >  4:238041/1‑148 (MQ=255)                                                                                                                                                   aaGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGcc    >  3:56665/1‑149 (MQ=255)                                                                                                                                                     gCCTTGAATCTTGCTAGACACCTGTGGGAGCGGGCTTGCCCGCGAACTGCCTCACCTCAATCCAATGCCTTCCAGATCTCTGTGGCATATTCTCGAATGGTGCTGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCtt  <  5:314382/149‑1 (MQ=255)                                                                                                                                                     gCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCtt  >  8:207325/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                   TACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTT  >  NC_002947/5743212‑5743506

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