Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,713,788	T→C	6.8%	W674R (TGG→CGG)	PP_4171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,713,788	0	T	C	6.8%	207.8 / 4.5	73	W674R (TGG→CGG)	PP_4171	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base C (2/3);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

GTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGC  >  NC_002947/4713640‑4713931                                                                                                                                                     |                                                                                                                                                gTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCACCCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCtgt                                                                                                                                                 <  4:312886/149‑1 (MQ=255)     aTGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCtgt                                                                                                                                                 <  1:11775/145‑1 (MQ=255)     aTGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCtgt                                                                                                                                                 >  2:11775/1‑145 (MQ=255)     aTGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTgg                                                                                                                                               >  1:167933/1‑147 (MQ=255)       gcgcGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgcagc                                                                                                                                           >  1:301491/1‑149 (MQ=255)          cgACTGGGGCGCCGGCGGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCagcagc                                                                                                                                        <  8:242469/149‑1 (MQ=255)              tGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCt                                                                                                                                    <  2:179188/149‑1 (MQ=255)                gggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGCAGCCCTCAACCAGTACCGCTACGCGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAGCGccgc                                                                                                                                   <  7:122833/148‑3 (MQ=255)                gggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAGCGccgc                                                                                                                                   <  7:101538/148‑3 (MQ=255)                 ggTGCCTGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTTGa                                                                                                                                 <  3:5329/149‑1 (MQ=255)                 ggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGa                                                                                                                                 >  7:83171/1‑149 (MQ=255)                   tGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCACCCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGCCg                                                                                                                               >  4:209184/1‑149 (MQ=255)                        gCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGaa                                                                                                                          >  1:272805/1‑149 (MQ=255)                              cGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCAcc                                                                                                                    >  4:252905/1‑149 (MQ=255)                                  tgcgcTGTTCTATTACGGAGCCGAGCCCCTCACCCAGTACCGCTACGAGCTTCTGGACGTGCGCCTGGCCAAGCGTTTTCGGGTGTATTGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCgg                                                                                                                <  4:247944/148‑1 (MQ=255)                                  cgcgcTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCgg                                                                                                                >  2:300428/1‑149 (MQ=255)                                     gcTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGGAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGccc                                                                                                             >  8:187072/1‑149 (MQ=255)                                      cTGTTCTATTACGGTCACGTTGCCCTCTACCTGTAGCTCTACGTGCGTTTGGAGGTGCGCCTGGCCAAGCGTTTTCGGGTGTTTGGCAGCAACCTGGAGCTGGGGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCt                                                                                                            <  2:167933/149‑1 (MQ=255)                                       tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAgcgc                                                                                                                                      >  4:234399/1‑122 (MQ=255)                                       tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAgcgc                                                                                                                                      <  3:234399/122‑1 (MQ=255)                                         ttCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCa                                                                                                         <  3:227150/149‑1 (MQ=255)                                           cTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTTCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAg                                                                                                        >  3:202538/1‑148 (MQ=255)                                            tATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGTAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGaa                                                                                                      >  1:92864/1‑149 (MQ=255)                                              ttACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAAcc                                                                                                    <  1:300428/149‑1 (MQ=255)                                                   gTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGTAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTAc                                                                                               <  2:92864/149‑1 (MQ=255)                                                   gTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTAc                                                                                               <  2:272805/149‑1 (MQ=255)                                                   gTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGGTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCGGTGGCAGCAGCGCCGGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTAc                                                                                               <  8:58816/149‑1 (MQ=255)                                                       cGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACa                                                                                           >  3:136172/1‑149 (MQ=255)                                                          tGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGcc                                                                                        <  7:187072/149‑1 (MQ=255)                                                           gCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCg                                                                                       >  2:231610/1‑149 (MQ=255)                                                             ccTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGt                                                                                      >  2:4196/1‑148 (MQ=255)                                                             ccTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGt                                                                                      >  2:242718/1‑148 (MQ=255)                                                              cTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCa                                                                                    >  4:182598/1‑149 (MQ=255)                                                                   ccAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGAt                                                                                                                             >  2:116273/1‑103 (MQ=255)                                                                   ccAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGAt                                                                                                                             <  1:116273/103‑1 (MQ=255)                                                                    cAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCt                                                                                                                                    >  6:331793/1‑95 (MQ=255)                                                                           gCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtg                                                                       <  3:182598/149‑1 (MQ=255)                                                                            cTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtga                                                                      <  1:322321/149‑1 (MQ=255)                                                                            cTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgt                                                                        >  6:131570/1‑147 (MQ=255)                                                                               cGATCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGTGTTTATTGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACTATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtgag                                                                     <  4:202538/147‑1 (MQ=255)                                                                                      cTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtga                                                                            <  4:153256/133‑1 (MQ=255)                                                                                      cTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtga                                                                            >  3:153256/1‑133 (MQ=255)                                                                                      cTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCt                                                            >  3:127218/1‑149 (MQ=255)                                                                                          aCCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAg                                                        <  1:130164/149‑1 (MQ=255)                                                                                              gcgcCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCt                                                    >  6:127831/1‑149 (MQ=255)                                                                                                 ccTGGCCAAGCGTGTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCGGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGAt                                                 <  3:209184/149‑1 (MQ=255)                                                                                                  cTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAg                                                        <  4:292875/141‑1 (MQ=255)                                                                                                  cTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAg                                                        >  3:292875/1‑141 (MQ=255)                                                                                                    ggCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGc                                              >  3:174345/1‑149 (MQ=255)                                                                                                        aaGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctg                                          >  5:252011/1‑149 (MQ=255)                                                                                                         aGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgc                                         >  1:314376/1‑149 (MQ=255)                                                                                                               ttCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCtgctgctgc                                   >  2:179959/1‑149 (MQ=255)                                                                                                                 cGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTg                                                           >  6:248776/1‑123 (MQ=255)                                                                                                                 cGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTg                                                           <  5:248776/123‑1 (MQ=255)                                                                                                                  ggggGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTgcgc                                 <  4:136172/148‑1 (MQ=255)                                                                                                                  gggTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTgcgc                                 <  3:252905/148‑1 (MQ=255)                                                                                                                     tgtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGt                             <  4:127218/149‑1 (MQ=255)                                                                                                                      gtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGaa                                                                                                      >  4:314583/1‑75 (MQ=255)                                                                                                                      gtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGaa                                                                                                      <  3:314583/75‑1 (MQ=255)                                                                                                                      gtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGtt                            >  2:206267/1‑149 (MQ=255)                                                                                                                       tATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtg                           <  1:206267/149‑1 (MQ=255)                                                                                                                        aTGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCtgctgctgcac                                 <  5:250466/142‑3 (MQ=255)                                                                                                                        aTGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCtgctgctgcac                                 >  6:250466/1‑140 (MQ=255)                                                                                                                            cagcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctg                                          >  5:244248/1‑129 (MQ=255)                                                                                                                            cagcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctg                                          <  6:244248/129‑1 (MQ=255)                                                                                                                             agcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCATCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGctc                     >  2:126116/1‑149 (MQ=255)                                                                                                                                 aCCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCAGCGCCTGTTGCTGCTCTGcc                 >  6:308179/1‑149 (MQ=255)                                                                                                                                   cTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCtg               >  6:218924/1‑149 (MQ=255)                                                                                                                                       aGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtgc                          <  6:154439/134‑1 (MQ=255)                                                                                                                                        gCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGcc          <  2:314376/149‑1 (MQ=255)                                                                                                                                           ggcggcGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGtt       <  2:301491/149‑1 (MQ=255)                                                                                                                                            gcggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGTAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGtt       <  1:207675/148‑1 (MQ=255)                                                                                                                                             cggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTgg     <  4:174345/149‑1 (MQ=255)                                                                                                                                              ggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGtt                                                      <  3:264269/99‑1 (MQ=255)                                                                                                                                              ggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGtt                                                      >  4:264269/1‑99 (MQ=255)                                                                                                                                                cgcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGc  <  8:83171/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                GTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGC  >  NC_002947/4713640‑4713931

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