Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,614,706	G→T	8.7%	A272A (GCC→GCA)	hldE ←	bifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,614,706	0	G	T	8.7%	214.9 / 9.1	69	A272A (GCC→GCA)	hldE	bifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base T (3/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

GCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAGGCGCCCAG  >  NC_002947/5614562‑5614852                                                                                                                                                 |                                                                                                                                                    gccgcGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGTcg                                                                                                                                                 <  4:271624/149‑1 (MQ=255) gccgcGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCcg                                                                                                                                                 <  5:26445/149‑1 (MQ=255)  ccgcGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCcgc                                                                                                                                                <  7:46661/149‑1 (MQ=255)  ccgcGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCcgc                                                                                                                                                <  4:302752/149‑1 (MQ=255)    gcgTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCCATTGCCGCgg                                                                                                                                              <  1:134327/149‑1 (MQ=255)      gTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCCATTGCCGCgggt                                                                                                                                            <  6:319086/149‑3 (MQ=255)           gAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGt                                                                                                                                       >  3:45119/1‑149 (MQ=255)            aGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTg                                                                                                                                      >  6:65730/1‑149 (MQ=255)             gCCGTCCTCACGCTGGCTGGCCCGGCGCAGCTCAGGGGCGCGGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTgg                                                                                                                                     <  4:249057/149‑1 (MQ=255)              ccTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTgg                                                                                                                                     >  5:323415/1‑148 (MQ=255)                ttCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGCCCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGaaa                                                                                                                                  >  7:109524/1‑149 (MQ=255)                ttCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGaaa                                                                                                                                  >  2:32540/1‑149 (MQ=255)                  ccTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAAt                                                                                                                                 <  5:328696/148‑1 (MQ=255)                    tCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAAAGCCCGCGGCAAGGGTGGAAATGa                                                                                                                               >  7:267998/1‑148 (MQ=255)                             aTGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑Gcc                                                                                                                    >  5:238405/1‑149 (MQ=255)                                gCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCgg                                                                                                                  >  1:59618/1‑148 (MQ=255)                                 cccGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑Gc                                                                                                                     <  6:299120/144‑1 (MQ=255)                                 cccGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑Gc                                                                                                                     >  5:299120/1‑144 (MQ=255)                                  gcGGCGCAGCTCAGGGTCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCGGCTGCCAGGTTAGGCAGGGCCTCCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCAc                                                                                                               <  8:109524/148‑1 (MQ=255)                                      cgcAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGt                                                                                                           <  5:260541/149‑1 (MQ=255)                                              aGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCg                                                                                                    <  2:23825/148‑1 (MQ=255)                                                 ggCGCTGATGGCAGCCGTACACAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACACTATC‑GCCGGCACCGGTCACAGCGAaca                                                                                                >  5:333378/1‑149 (MQ=255)                                                      tGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCTAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGaa                                                                                                  <  1:99407/142‑1 (MQ=255)                                                      tGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGaa                                                                                                  >  2:99407/1‑142 (MQ=255)                                                               ccGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCa                                                                                                                <  4:229680/119‑1 (MQ=255)                                                               ccGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCa                                                                                                                >  3:229680/1‑119 (MQ=255)                                                                      cAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCgcagg                                                                             >  5:109058/1‑147 (MQ=255)                                                                      cAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCgcagg                                                                             >  1:269812/1‑147 (MQ=255)                                                                      cAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCgcagg                                                                             >  1:122355/1‑147 (MQ=255)                                                                         cTTGACGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGt                                                                        <  3:208465/149‑1 (MQ=255)                                                                             ccGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCaa                                                                    >  2:15105/1‑149 (MQ=255)                                                                                tccccaATGCCCGCTGCCTGGTACCCCCGGGCAACCGCTTGGGGAAGGGCCCCCCCCCCGCCATTGCCGGGCGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGc                                                                 <  2:141429/145‑1 (MQ=255)                                                                                  cacaATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGc                                                                 >  1:141429/1‑147 (MQ=255)                                                                                   acaATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCgg                                                              <  2:269812/149‑1 (MQ=255)                                                                                     aaTGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGGCCTCCCCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGATCACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCt                                                            <  8:184871/149‑1 (MQ=255)                                                                                      aTGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTg                                                           >  6:257158/1‑149 (MQ=255)                                                                                       tGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGa                                                          <  2:59618/149‑1 (MQ=255)                                                                                             cTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAaca                                                    <  5:103047/149‑1 (MQ=255)                                                                                              tGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGCCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCAc                                                                                                               >  6:163170/1‑89 (MQ=255)                                                                                                ccAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAacac                                                   <  6:238405/147‑1 (MQ=255)                                                                                                  aGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCaa                                               <  1:15105/149‑1 (MQ=255)                                                                                                        gCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGt                                         >  3:50762/1‑149 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCgcgc                                                                                <  4:199127/109‑1 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCgcgc                                                                                >  3:199127/1‑109 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAac                                                     <  3:136000/136‑1 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAac                                                     >  4:136000/1‑136 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAgg                                           <  4:130727/146‑1 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAgg                                           >  3:130727/1‑146 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTc                                        <  1:32540/149‑1 (MQ=255)                                                                                                              gCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTcatgcc                                   <  3:303221/149‑1 (MQ=255)                                                                                                                  ccGGATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTcatgccatgc                               <  7:223683/149‑1 (MQ=255)                                                                                                                  ccGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTcatgccatgc                               >  1:170917/1‑149 (MQ=255)                                                                                                                     cATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGAc                                                         <  4:83516/120‑1 (MQ=255)                                                                                                                     cATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGAc                                                         >  3:83516/1‑120 (MQ=255)                                                                                                                        ggggCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCAcc                          >  1:88276/1‑149 (MQ=255)                                                                                                                           gCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCgg                      <  6:109058/149‑1 (MQ=255)                                                                                                                           gCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGACATGCTCACCGCgg                      >  4:135585/1‑149 (MQ=255)                                                                                                                                  cTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCa               >  4:41968/1‑149 (MQ=255)                                                                                                                                    ttgccGCGGCCATTGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTcatgccat                                 <  3:154151/126‑1 (MQ=255)                                                                                                                                    ttgccGCGGCCATTGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTcatgccat                                 >  4:154151/4‑129 (MQ=255)                                                                                                                                    ttgccGCGGCCATTGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAg              >  6:7535/4‑148 (MQ=255)                                                                                                                                     aCCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGc                                                             <  7:262536/100‑1 (MQ=255)                                                                                                                                     aCCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGc                                                             >  8:262536/1‑100 (MQ=255)                                                                                                                                       ccGCGGCCATTGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAg          >  4:164686/1‑149 (MQ=255)                                                                                                                                       ccGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGAc                                                         >  6:285519/1‑102 (MQ=255)                                                                                                                                       ccGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGAc                                                         <  5:285519/102‑1 (MQ=255)                                                                                                                                       ccGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAg          >  7:190068/1‑149 (MQ=255)                                                                                                                                             cAATGGCCGCGGCAAGGGGGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAGGCGccc    <  8:62216/149‑1 (MQ=255)                                                                                                                                               aTGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAGGCGCccag  >  3:102977/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAGGCGCCCAG  >  NC_002947/5614562‑5614852

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