Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,184,110	A→C	5.4%	M83L (ATG→CTG)	PP_4564 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,184,110	0	A	C	5.4%	286.4 / 6.3	92	M83L (ATG→CTG)	PP_4564	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (43/44);  new base C (3/2);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

CCGCCTGGTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGAT  >  NC_002947/5183962‑5184255                                                                                                                                                     |                                                                                                                                                   ccGCCTGGTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACa                                                                                                                                                    >  8:2617/1‑149 (MQ=255)  cGCCTGGTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGCCCGGCTACGTGACAt                                                                                                                                                   >  5:240598/1‑149 (MQ=255)  cGCCTGGTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACAt                                                                                                                                                   >  5:234554/1‑149 (MQ=255)    ccTGGTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATgg                                                                                                                                                 >  5:159571/1‑149 (MQ=255)       ggTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCg                                                                                                                                              >  4:21023/1‑149 (MQ=255)        gTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGc                                                                                                                                             <  1:14552/149‑1 (MQ=255)        gTGCCGTGGGCGCTGAGCGACTTCCTGTACATGAACCGCGAGCTAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGc                                                                                                                                             <  4:344629/149‑1 (MQ=255)                 tCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCa                                                                                                                                    <  5:287601/149‑1 (MQ=255)                 tCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCa                                                                                                                                    <  5:120388/149‑1 (MQ=255)                   gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTATCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAg                                                                                                                                   <  1:307518/148‑1 (MQ=255)                   gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAgg                                                                                                                                  >  2:410300/1‑149 (MQ=255)                    cTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                 >  6:115422/1‑149 (MQ=255)                       aGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTc                                                                                                                                <  5:103607/147‑1 (MQ=255)                       aGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCa                                                                                                                               <  3:21023/148‑1 (MQ=255)                        gCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCAcc                                                                                                                             <  5:199509/149‑1 (MQ=255)                           aCTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCTGCTACGTGACATGGTCGCCGACGACCACGTCACCgcc                                                                                                                          >  1:333904/1‑147 (MQ=255)                           aCTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑c                                                                                                                         <  2:54186/149‑1 (MQ=255)                               ccTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCa                                                                                                                       <  4:63635/147‑1 (MQ=255)                                   gACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGccc                                                                                                                  <  1:112659/148‑1 (MQ=255)                                    aCATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTg                                                                                                                <  1:111365/149‑1 (MQ=255)                                     cATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGCCCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTg                                                                                                                >  8:311079/1‑148 (MQ=255)                                      aTGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTg                                                                                                                >  6:439505/1‑147 (MQ=255)                                      aTGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTg                                                                                                                >  2:92227/1‑147 (MQ=255)                                       tGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAgg                                                                                                                     >  7:129294/1‑141 (MQ=255)                                       tGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAgg                                                                                                                     <  8:129294/141‑1 (MQ=255)                                        gAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGACTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTgcag                                                                                                             >  5:95098/1‑148 (MQ=255)                                         aaCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTgcagc                                                                                                            >  8:23099/1‑148 (MQ=255)                                          aCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCaa                                                                                                          <  5:115422/149‑1 (MQ=255)                                           ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                 >  2:294258/1‑126 (MQ=255)                                           ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                 <  1:294258/126‑1 (MQ=255)                                           ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCGACGTGACATGGGCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCaaa                                                                                                         <  6:117039/149‑1 (MQ=255)                                            cgcgAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAg                                                                                                        >  7:169995/1‑149 (MQ=255)                                              cgtgctaaGGGCGCTTCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTGGCGGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAATGCa                                                                                                      <  2:333904/143‑1 (MQ=255)                                                     aGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATgg                                                                                                                                                 <  5:13235/100‑1 (MQ=255)                                                     aGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATgg                                                                                                                                                 >  6:13235/1‑100 (MQ=255)                                                            tgctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGt                                                                                        >  6:377288/1‑149 (MQ=255)                                                              ctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTg                                                                                       <  2:14498/148‑1 (MQ=255)                                                              ctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGa                                                                                      <  6:159571/149‑1 (MQ=255)                                                               tgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGa                                                                                      <  1:290632/149‑1 (MQ=255)                                                                ggATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAgg                                                                                    >  2:385769/1‑149 (MQ=255)                                                                ggATGATCGGTTGCGCGAGCACCTGACCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAgg                                                                                    >  5:203104/1‑149 (MQ=255)                                                                 gatgatCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTgca                                                                                                              >  2:307170/1‑122 (MQ=255)                                                                 gatgatCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTgca                                                                                                              <  1:307170/122‑1 (MQ=255)                                                                 gatgatCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAgg                                                                                    <  1:385769/148‑1 (MQ=255)                                                                        ggTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGCCGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCCAGGTCCGTGAGGCAATTGGt                                                                            >  8:348885/1‑149 (MQ=255)                                                                        ggTGGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGt                                                                            <  4:300735/149‑1 (MQ=255)                                                                         gTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGGCCGGGTAACTAGAGTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTGACCGA‑CCAGGCCGTGCAGCAAAGCACCCGCGAGGTCGGTGAGGCAATTGGTc                                                                           >  7:272786/1‑149 (MQ=255)                                                                         gTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTc                                                                           >  1:188588/1‑149 (MQ=255)                                                                         gTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGACGCAATTGGTc                                                                           >  7:270559/1‑149 (MQ=255)                                                                             cgcgAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                       <  5:109355/149‑1 (MQ=255)                                                                             cgcgAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                       <  7:23099/149‑1 (MQ=255)                                                                             cgcgAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                       <  4:178382/149‑1 (MQ=255)                                                                              gcgAGCACCTGGCCTGGCATTGCAAGACTCCACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGGGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTg                                                                      <  4:242945/149‑1 (MQ=255)                                                                                     cctggcctggCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAgg                                                               <  2:301517/149‑1 (MQ=255)                                                                                      ctggcctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAgg                                                               >  8:84277/1‑148 (MQ=255)                                                                                         gcctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGACTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGaaa                                                            <  6:95098/148‑1 (MQ=255)                                                                                           ctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCAACGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCa                                                         >  2:189837/1‑149 (MQ=255)                                                                                            tggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAc                                                        >  7:19058/1‑149 (MQ=255)                                                                                              gcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAccc                                                      >  5:127037/1‑149 (MQ=255)                                                                                              gcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAcc                                                       <  1:410300/148‑1 (MQ=255)                                                                                                   gCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCg                                                 >  8:371213/1‑149 (MQ=255)                                                                                                    cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGTCATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGa                                                                 >  2:44089/1‑132 (MQ=255)                                                                                                    cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGTCATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGa                                                                 <  1:44089/132‑1 (MQ=255)                                                                                                    cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGGCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCgg                                                <  7:128058/149‑1 (MQ=255)                                                                                                     aaGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGt                                                                                        >  6:167241/1‑108 (MQ=255)                                                                                                         cTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAcc                                           >  5:366026/1‑149 (MQ=255)                                                                                                            aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGaa                                        >  3:331454/1‑149 (MQ=255)                                                                                                             aCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAAc                                       <  6:240598/149‑1 (MQ=255)                                                                                                             aCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAAc                                       <  6:234554/149‑1 (MQ=255)                                                                                                               tGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCAAGGAAATCACCCCTTCGGCCACCGAActg                                     >  1:61530/1‑149 (MQ=255)                                                                                                                  ccGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGACctgctg                                  >  2:87434/1‑149 (MQ=255)                                                                                                                  ccGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAActgctg                                  >  3:21963/1‑149 (MQ=255)                                                                                                                   cGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACtgctgc                                 >  1:328344/1‑149 (MQ=255)                                                                                                                   cGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAAACACCCCTTCGGCCACCGAACtgctgc                                 >  1:65327/1‑149 (MQ=255)                                                                                                                          cTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCAt                          <  4:21963/149‑1 (MQ=255)                                                                                                                           tCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTCCTGCGCGGCAt                          >  4:315705/1‑148 (MQ=255)                                                                                                                              aCTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCACGAGc                      <  5:377288/149‑1 (MQ=255)                                                                                                                                  gCTCGACCGGCTACGTGACAGGGTCGCCGACGCCCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGa                    <  8:307155/147‑1 (MQ=255)                                                                                                                                    tCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAAc                >  7:375136/1‑149 (MQ=255)                                                                                                                                      gACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAAcc               >  6:364971/1‑148 (MQ=255)                                                                                                                                      gACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAAcc               <  5:364971/148‑1 (MQ=255)                                                                                                                                         ctggtcggTGACCTGGTCGTCGGCGACCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGt           >  1:433575/8‑149 (MQ=39)                                                                                                                                            cttcGTTACATGGTCGACGACGACATGGTCAGCGA‑GCAGGGCCTGCAGCATACCAGCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGATATGACCGCTTCCGCCACGGAACTGCTGCGCGGCATGAGCGACAAGCAGGTAgc        <  8:270559/146‑1 (MQ=255)                                                                                                                                            cTACGTGACATGGTCGCCGACGTCCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCTATTGTCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgc        <  8:272786/149‑1 (MQ=255)                                                                                                                                               cgggACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCa     <  2:179340/146‑1 (MQ=255)                                                                                                                                               cGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgcg       <  5:439505/147‑1 (MQ=255)                                                                                                                                                gTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCACCAAAGCCCCCGCGACGTCCGTGAGGCAATTCGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTCCGCGGCATCAGAGACAACCAGGAAgcc       >  4:361579/1‑145 (MQ=255)                                                                                                                                                 tGACCTGGTCGTCGGCGACCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgcg       >  1:185646/1‑145 (MQ=255)                                                                                                                                                 tGACCTGGTCGTCGGCGACCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAg         >  7:299233/1‑143 (MQ=255)                                                                                                                                                 tGACCTGGTCGTCGGCGACCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAg         <  8:299233/143‑1 (MQ=255)                                                                                                                                                 tGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgcg       >  3:182411/1‑145 (MQ=255)                                                                                                                                                 tGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGa   <  6:366026/149‑1 (MQ=255)                                                                                                                                                   aCATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGAt  <  8:375136/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CCGCCTGGTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGAT  >  NC_002947/5183962‑5184255

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