Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	474,192	Δ1 bp	5.8%	coding (793/1983 nt)	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	474,192	0	C	.	5.8%	219.5 / 6.4	69	coding (793/1983 nt)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base . (2/2);  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 = 1.00e+00

TGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGAT  >  NC_002947/474044‑474336                                                                                                                                                     |                                                                                                                                                  tGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTcgacgac                                                                                                                                                   >  4:301889/1‑149 (MQ=255)       ccGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACgatgat                                                                                                                                             >  2:46372/1‑149 (MQ=255)                 agaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGA‑GATCATCGTCGTCGAgg                                                                                                                                  <  2:305527/149‑1 (MQ=255)                   aacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGACGtt                                                                                                                                >  5:293120/1‑149 (MQ=255)                      agCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCg                                                                                                                             >  8:106243/1‑149 (MQ=255)                        cACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACgatg                                                                                                                                               <  8:1992/130‑1 (MQ=255)                        cACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACgatg                                                                                                                                               >  7:1992/1‑130 (MQ=255)                         aCGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGtt                                                                                                                          <  5:253837/149‑1 (MQ=255)                           gCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTtgt                                                                                                                        <  1:46372/149‑1 (MQ=255)                                   ccacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgc                                                                                                                <  3:268745/149‑1 (MQ=255)                                    cacGCGGAACAGGCGCATGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgcg                                                                                                               >  6:58352/1‑149 (MQ=255)                                     acGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgcgc                                                                                                              >  4:36106/1‑149 (MQ=255)                                     acGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgcgc                                                                                                              >  3:179323/1‑149 (MQ=255)                                        cgGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTc                                                                                                           >  6:228282/1‑149 (MQ=255)                                        cgGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTc                                                                                                           >  6:233368/1‑149 (MQ=255)                                              aGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAg                                                                                                     >  7:202003/1‑149 (MQ=255)                                                gcgcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGcc                                                                                                   <  3:36106/149‑1 (MQ=255)                                                gcgcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGcc                                                                                                   >  4:219197/1‑149 (MQ=255)                                                 cgcTTGAGGGGCTCTGCGCGGGTGGCGGGGCCAAGGGTGGCCACGGCATGGCGCAGGCCTTGCTGGGCCAGGGCATTGACGTACATGTAGCCCGCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGccc                                                                                                  <  8:243916/149‑1 (MQ=255)                                                         gTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACGGGCTGTAcgag                                                                                           <  5:233368/148‑4 (MQ=255)                                                         gTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACGGGCTGTAcgag                                                                                           <  5:228282/148‑4 (MQ=255)                                                         gTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCt                                                                                          >  2:67241/1‑149 (MQ=255)                                                           gCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTcc                                                                                        >  5:144853/1‑149 (MQ=255)                                                              ctTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGc                                                                                     >  4:147048/1‑149 (MQ=255)                                                                 cGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGccc                                                                                   <  5:58352/148‑1 (MQ=255)                                                                  gCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCt                                                                                 >  8:97370/1‑149 (MQ=255)                                                                      gtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgtg                                                                             <  3:301889/149‑1 (MQ=255)                                                                      gtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                              <  1:53831/148‑1 (MQ=255)                                                                       tggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTgg                                                                            <  1:67241/149‑1 (MQ=255)                                                                       tggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTgg                                                                            <  3:148938/149‑1 (MQ=255)                                                                        ggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgc                                                                                                                >  6:334398/1‑112 (MQ=255)                                                                        ggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgc                                                                                                                <  5:334398/112‑1 (MQ=255)                                                                         gCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGa                                                                           >  2:225542/1‑148 (MQ=255)                                                                           ggTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCt                                                                                       >  4:197315/1‑134 (MQ=255)                                                                           ggTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCt                                                                                       <  3:197315/134‑1 (MQ=255)                                                                               cgaaGGGGGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                    <  5:126097/147‑1 (MQ=255)                                                                                 aaGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGTAAACCggg                                                                   <  8:240735/148‑1 (MQ=255)                                                                                 aaGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGGAAACCggg                                                                   >  7:240735/1‑148 (MQ=255)                                                                                 aaGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgggg                                                                  >  5:182232/1‑149 (MQ=255)                                                                                      tGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGt                                                                                                <  8:19501/114‑1 (MQ=255)                                                                                      tGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGt                                                                                                >  7:19501/1‑114 (MQ=255)                                                                                         ccACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCggg                                                          >  5:286257/1‑149 (MQ=255)                                                                                             ggCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCggg                                                          >  3:110889/1‑145 (MQ=255)                                                                                             ggCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCggg                                                          <  4:110889/145‑1 (MQ=255)                                                                                              gCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGtt                                                     <  6:286257/149‑1 (MQ=255)                                                                                               cATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAAGTTGCGGGTGTGCTTGCGCGCCGCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTc                                                    <  6:293120/149‑1 (MQ=255)                                                                                                 ttGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgc                                                                                                                <  4:88639/87‑1 (MQ=255)                                                                                                 ttGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgc                                                                                                                >  3:88639/1‑87 (MQ=255)                                                                                                 ttGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                    >  3:205466/1‑131 (MQ=255)                                                                                                 ttGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                    <  4:205466/131‑1 (MQ=255)                                                                                                    cgcAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                   <  4:251585/129‑1 (MQ=255)                                                                                                    cgcAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                   >  3:251585/1‑129 (MQ=255)                                                                                                          ccTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTggg                                         <  6:182232/149‑1 (MQ=255)                                                                                                           cTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGc                                        >  3:117223/1‑149 (MQ=255)                                                                                                              gCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCtt                                      <  3:43066/148‑1 (MQ=255)                                                                                                                      aGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGcc                             <  4:179323/149‑1 (MQ=255)                                                                                                                       gggCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTAGAGGTTGCGGTTGTGCTTGCGCGCCTAGTACAGACCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGcc                             >  3:279685/1‑148 (MQ=255)                                                                                                                        ggCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGcc                             >  4:278293/1‑147 (MQ=255)                                                                                                                                cGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCg                    <  7:97370/148‑1 (MQ=255)                                                                                                                                 gTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGc                  <  8:202003/149‑1 (MQ=255)                                                                                                                                 gTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGc                  >  5:99992/1‑149 (MQ=255)                                                                                                                                 gTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGc                  >  7:200250/1‑149 (MQ=255)                                                                                                                                   ccATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCtt                                      <  7:278458/127‑1 (MQ=255)                                                                                                                                   ccATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCtt                                      >  8:278458/1‑127 (MQ=255)                                                                                                                                           cccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGTTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGc        <  6:144853/149‑1 (MQ=255)                                                                                                                                            ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCggg                                                          <  5:123322/98‑1 (MQ=255)                                                                                                                                            ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCggg                                                          >  6:123322/1‑98 (MQ=255)                                                                                                                                            ccTCGACGA‑GATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGc        >  1:104657/1‑148 (MQ=255)                                                                                                                                            ccTCGACGA‑GATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCg       >  5:50049/1‑149 (MQ=255)                                                                                                                                             cTCGACGA‑GATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCga      <  2:104657/149‑1 (MQ=255)                                                                                                                                                gacgacgaTGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgat     <  6:99992/147‑1 (MQ=255)                                                                                                                                                 acgacgaTGATC‑TCGTCGAGGTTGCGGTTGTGCTTGAGCGCATCGTACAGCCCTTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatgat  <  6:316221/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGAT  >  NC_002947/474044‑474336

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