Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,237,544	A→G	10.7%	S2G (AGC→GGC)	PP_4614 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,237,544	0	A	G	10.7%	170.9 / 8.9	65	S2G (AGC→GGC)	PP_4614	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base G (2/5);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAA  >  NC_002947/5237400‑5237692                                                                                                                                                 |                                                                                                                                                      ttagCTGGACCAGGCAGTCGACGGAGCCCCCACCAAGCTCAATCAAGCGCTGGAGCACTTCTACGCCAAGCTGCGCAGCAAGCGCGGCGCCTACGAACTGAGAGACCCTGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCg                                                                                                                                                   >  4:32056/4‑149 (MQ=255)   cGCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCgg                                                                                                                                                  <  5:285324/148‑1 (MQ=255)      tGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGt                                                                                                                                              >  1:85536/1‑149 (MQ=255)        gACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCTCAAGGCCGGGCGtt                                                                                                                                             <  3:347488/148‑1 (MQ=255)        gACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCTCAAGGCCGGGCGtt                                                                                                                                             <  8:27904/148‑1 (MQ=255)        gACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGttt                                                                                                                                            <  8:118818/149‑1 (MQ=255)        gACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGttt                                                                                                                                            >  5:210738/1‑149 (MQ=255)        gACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGt                                                                                                                                              <  2:89220/147‑1 (MQ=255)         accctgCAGTCGACGGGGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCAACCGCAAGCTGCGCAGCAAGCGCTGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCTCAAGGCCGGGCGtttt                                                                                                                                           <  8:68052/144‑1 (MQ=255)         aCCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGtttt                                                                                                                                           >  7:37767/1‑149 (MQ=255)            aGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGTAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCTCCTGGGAACTGCGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑Gt                                                                                                                                       <  7:122452/149‑1 (MQ=255)                aGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTtgt                                                                                                                                   >  5:47045/1‑149 (MQ=255)                  tCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGc                                                                                                                                 <  2:246113/149‑1 (MQ=255)                  tCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGc                                                                                                                                 <  8:37767/149‑1 (MQ=255)                           ccgccACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATg                                                                                                                        >  1:40335/1‑149 (MQ=255)                            cgccACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATgg                                                                                                                       <  7:21453/149‑1 (MQ=255)                              ccaccaAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGgcg                                                                                                                     <  2:85536/149‑1 (MQ=255)                                accaAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAAt                                                                                                                         >  6:41244/1‑143 (MQ=255)                                accaAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAAt                                                                                                                         <  5:41244/143‑1 (MQ=255)                                             aCGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGc                                                                                                       >  2:358630/1‑148 (MQ=255)                                                      aGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCgg                                                                                             <  4:20764/149‑1 (MQ=255)                                                      aGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCgg                                                                                             <  3:206092/149‑1 (MQ=255)                                                          cTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCaa                                                                                         >  2:76723/1‑149 (MQ=255)                                                           ttCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAAc                                                                                        <  1:76723/149‑1 (MQ=255)                                                           ttCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAAc                                                                                        >  6:362207/1‑149 (MQ=255)                                                             cTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAAc                                                                                        <  5:362207/147‑1 (MQ=255)                                                             atACGGCAAGGTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGTCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAAcc                                                                                       <  3:32056/147‑1 (MQ=255)                                                                 ggCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTg                                                                                   >  3:14533/1‑148 (MQ=255)                                                                   cAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCTGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCAt                                                                                <  7:12197/149‑1 (MQ=255)                                                                    aaGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCAtc                                                                               <  2:40335/149‑1 (MQ=255)                                                                     aGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCAtca                                                                              <  4:269046/149‑1 (MQ=255)                                                                       cTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCAtcact                                                                            <  3:254319/149‑1 (MQ=255)                                                                        tGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCAtcact                                                                            >  2:133440/1‑148 (MQ=255)                                                                            cagcaAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCg                                                                       >  2:56244/1‑149 (MQ=255)                                                                               caAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGGGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCg                                                                       >  7:66913/1‑146 (MQ=255)                                                                                          ccTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCa                                                         <  3:282761/149‑1 (MQ=255)                                                                                                         aCCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTc                                           >  4:204261/1‑148 (MQ=255)                                                                                                         aCCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCATGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTc                                           >  1:176439/1‑148 (MQ=255)                                                                                                                    tGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGc                                                                                                                                 <  4:93486/51‑1 (MQ=255)                                                                                                                    tGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGc                                                                                                                                 >  3:93486/1‑51 (MQ=255)                                                                                                                    tGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAGCCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCAGTCGAACGCATCAGt                               >  8:237800/1‑149 (MQ=255)                                                                                                                    tGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATAAGt                               >  5:232565/1‑149 (MQ=255)                                                                                                                       tCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGtgct                            >  3:112518/1‑149 (MQ=255)                                                                                                                          acaaacgaaaaCGCCCGGCTCACGGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTCAACCAACATGCGAACGTCCGGGCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGtgct                            <  6:337133/139‑1 (MQ=255)                                                                                                                          acaaacgaaaaCGCCCGGCTCAAGGCCGGGCGGTTTC‑GGGTGAGCCGATGAATGGCGCTCAACCAACATGCGAACGTCCGGTGCACCTCTGCCTACCTCACCGCAGACGCAAGACCATGAACAACCCATTCGAAGGAATCAGtgct                            >  5:337133/8‑146 (MQ=255)                                                                                                                           cAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAt                        <  8:66913/149‑1 (MQ=255)                                                                                                                            aaTGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAt                        >  1:53498/1‑149 (MQ=255)                                                                                                                            aaTGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATGAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAtt                       >  6:6318/1‑149 (MQ=255)                                                                                                                            aaTGCAAAACGCCCGGCCTTGAGCCGGGCGTTATC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAtt                       >  1:225184/1‑149 (MQ=255)                                                                                                                             aTGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTc                      <  1:133440/149‑1 (MQ=255)                                                                                                                             aTGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTc                      <  2:53498/149‑1 (MQ=255)                                                                                                                                cAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGcc                   <  2:225184/149‑1 (MQ=255)                                                                                                                                cAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGcc                   >  4:138909/1‑149 (MQ=255)                                                                                                                                  aaaCGCCCGGCTCAAGGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTc                                           >  6:357148/1‑123 (MQ=255)                                                                                                                                  aaaCGCCCGGCTCAAGGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTc                                           <  5:357148/123‑1 (MQ=255)                                                                                                                                      gCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGtgctg                           <  2:280757/135‑1 (MQ=255)                                                                                                                                      gCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGaa             >  8:166521/1‑149 (MQ=255)                                                                                                                                      gCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACAGGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGtgctg                           >  1:280757/1‑135 (MQ=255)                                                                                                                                         cGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATAc          >  6:129695/1‑149 (MQ=255)                                                                                                                                         cGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATAc          >  5:203643/1‑149 (MQ=255)                                                                                                                                          ggCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATAcc         <  1:358630/149‑1 (MQ=255)                                                                                                                                            ccTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCgg       <  4:14533/149‑1 (MQ=255)                                                                                                                                              atGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGt     <  2:176439/148‑1 (MQ=255)                                                                                                                                                gAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAt                        >  4:318207/1‑128 (MQ=255)                                                                                                                                                gAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAt                        <  3:318207/128‑1 (MQ=255)                                                                                                                                                 aGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGaa  <  1:56244/149‑1 (MQ=255)                                                                                                                                                 aGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGa   <  3:204261/148‑1 (MQ=255)                                                                                                                                                 aGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGa   <  6:232565/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      TTCGCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTC‑GTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAA  >  NC_002947/5237400‑5237692

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