Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,614,703	A→C	8.8%	I273M (ATT→ATG)	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,703	0	A	C	8.8%	202.5 / 10.3	68	I273M (ATT→ATG)	hldE	bifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (31/31);  new base C (3/3);  total (34/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 = 9.55e-01

GCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAGGCGCCC  >  NC_002947/5614562‑5614850                                                                                                                                              |                                                                                                                                                     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                                                                                                              >  3:229680/1‑119 (MQ=255)                                                               ccGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCa                                                                                                              <  4:229680/119‑1 (MQ=255)                                                                      cAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCgcagg                                                                           >  1:269812/1‑147 (MQ=255)                                                                      cAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCgcagg                                                                           >  5:109058/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                                                                              >  3:199127/1‑109 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCgcgc                                                                              <  4:199127/109‑1 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAac                                                   <  3:136000/136‑1 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAac                                                   >  4:136000/1‑136 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAgg                                         >  3:130727/1‑146 (MQ=255)                                                                                                         ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAgg                                         <  4:130727/146‑1 (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                               >  4:154151/4‑129 (MQ=255)                                                                                                                                    ttgccGCGGCCATTGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTcatgccat                               <  3:154151/126‑1 (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)                                                                                                                                              |                                                                                                                                                     GCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATC‑GCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAGGCGCCC  >  NC_002947/5614562‑5614850

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: ‑