Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	824,949	T→C	5.4%	L107P (CTG→CCG)	ybhQ →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	824,949	0	T	C	5.4%	271.9 / 4.8	92	L107P (CTG→CCG)	ybhQ	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (46/41);  new base C (3/2);  total (49/43)
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.37e-01

GTTTCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCC  >  NC_000913/824800‑825095                                                                                                                                                       |                                                                                                                                                   gtTTCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCt                                                                                                                                                    >  2:374861/1‑150 (MQ=255)    tCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGc                                                                                                                                                 <  1:735697/150‑1 (MQ=255)    ccAGCGCCACCCCGAGCAACGCTTGCGTGAGTTGGCAGACGTGCTGGAAGATCTGACCAAGACCTGGTA‑TTTTGGCTCAGGTCTGATTGGGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGc                                                                                                                                                 <  1:881273/149‑1 (MQ=255)     cAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGCAGCGCTGCCCACGACCTGGTA‑TGTTGGCGCCGCGCTGATGGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACGTTTTGCTGGCCATTGCCGGGCTGGcc                                                                                                                                                >  1:80416/1‑149 (MQ=255)     cAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCa                                                                                                                                                >  1:766820/1‑150 (MQ=255)     cAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCa                                                                                                                                                >  1:1253142/1‑150 (MQ=255)          ccACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCt                                                                                                                                           <  1:663380/150‑1 (MQ=255)          ccACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCt                                                                                                                                           <  2:1676763/150‑1 (MQ=255)             ccccGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGc                                                                                                                                        >  1:833158/1‑150 (MQ=255)             ccccGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCT‑GCAATTGCAGGGCTGGCAATCCTTGcc                                                                                                                                       <  2:1188352/150‑1 (MQ=255)              cccGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGcc                                                                                                                                       >  1:171188/1‑150 (MQ=255)               ccGAGCAACGCTGGCGTGAGGGGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGccg                                                                                                                                      <  2:1291335/150‑1 (MQ=255)                  aGCAACGCTGCCGTGAGGTGGCAGAAGTGCTGGAAGAGCTCACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTCTGCTGTGGCTGTCGTCCCGCGTTCTGGAAAACAACTATTAGCTGCCAATTGCAGGGCTGTCAATCCTTGccggcc                                                                                                                                   >  2:685830/1‑150 (MQ=255)                   gCAACGCGTGCGTGATGTGGCAGACGTGCTGGAAGAGCTGTCCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCCGGCAAGGATTGccagccc                                                                                                                                  <  2:637177/150‑1 (MQ=255)                    ccaCGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCCGGCAAGGATTGccagccct                                                                                                                                 <  2:1143140/148‑2 (MQ=255)                         cTGGCGGGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGt                                                                                                                            <  1:1275043/150‑1 (MQ=255)                           ggCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAg                                                                                                                          >  1:206408/1‑150 (MQ=255)                           ggCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGACCCGGAGGTAg                                                                                                                          >  1:137604/1‑150 (MQ=255)                             cGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAgtg                                                                                                                        <  2:1440431/150‑1 (MQ=255)                             cGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAgtg                                                                                                                        <  2:1732139/150‑1 (MQ=255)                               tGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTc                                                                                                                      <  1:1669817/150‑1 (MQ=255)                                gAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCt                                                                                                                     >  2:945964/1‑150 (MQ=255)                                 aGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCt                                                                                                                     >  1:863413/1‑150 (MQ=255)                                  ggTGGCAGACGTGCTGGAAGAGCTGACCACGACCTTGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCttt                                                                                                                   <  2:1667712/150‑1 (MQ=255)                                   gTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTg                                                                                                                  >  1:1201585/1‑150 (MQ=255)                                    tGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGc                                                                                                                 >  2:1728932/1‑150 (MQ=255)                                     ggCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCt                                                                                                                >  1:327447/1‑150 (MQ=255)                                       cAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTgg                                                                                                              >  2:820143/1‑150 (MQ=255)                                         gACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCg                                                                                                            <  2:919734/150‑1 (MQ=255)                                         gACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCg                                                                                                            <  2:426272/150‑1 (MQ=255)                                           cGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGaa                                                                                                          <  2:1650449/150‑1 (MQ=255)                                                 ggAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCTGCGGTAGTGTCTTTGCTGGCGAAAGATaa                                                                                                    >  1:804820/1‑150 (MQ=255)                                                    agagCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaagaa                                                                                                 >  2:1475836/1‑150 (MQ=255)                                                     gagCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaagaag                                                                                                >  1:614448/1‑150 (MQ=255)                                                     gagCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaagaag                                                                                                >  1:827687/1‑150 (MQ=255)                                                      agCTGACAACGCCGTGGTA‑TTTTGGCCCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGt                                                                                               <  1:301259/150‑1 (MQ=255)                                                         tGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTAc                                                                                            >  1:1214449/1‑150 (MQ=255)                                                          gACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACa                                                                                           >  2:504579/1‑150 (MQ=255)                                                          gACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGATGGCGAAAGATAAGAAGTTACa                                                                                           >  1:1430805/1‑150 (MQ=255)                                                           atcACGACCTGGTA‑TTTTGGCGCTGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAACACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACAt                                                                                          <  2:250299/148‑1 (MQ=255)                                                            ccACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACAta                                                                                         >  1:122540/1‑149 (MQ=255)                                                            ccACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATc                                                                                         <  2:517666/150‑1 (MQ=255)                                                             cACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCa                                                                                        >  2:753662/1‑150 (MQ=255)                                                                gACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCt                                                                                     >  2:456301/1‑150 (MQ=255)                                                                 aCCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCtt                                                                                    >  2:618604/1‑150 (MQ=255)                                                                  ccTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTa                                                                                   <  2:1280329/150‑1 (MQ=255)                                                                    tGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACg                                                                                 <  2:206408/150‑1 (MQ=255)                                                                    tGGTA‑TATTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACg                                                                                 <  2:31846/150‑1 (MQ=255)                                                                     ggtatttttGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGGTCTGGAAAACAACTTTTTGCTGTCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACg                                                                                 <  2:863413/145‑1 (MQ=255)                                                                     ggTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGt                                                                                >  2:364382/1‑150 (MQ=255)                                                                            ttGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAAc                                                                          <  1:374861/150‑1 (MQ=255)                                                                              ggCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACAt                                                                        >  2:552151/1‑150 (MQ=255)                                                                               gcgcAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATc                                                                       <  1:1475836/150‑1 (MQ=255)                                                                                cgcAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATcg                                                                      >  2:1314102/1‑150 (MQ=255)                                                                                  cAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATcgcg                                                                    >  2:771141/1‑150 (MQ=255)                                                                                         gATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcgccg                                                             >  1:91755/1‑150 (MQ=255)                                                                                            tgtgCTGTGGCTGGTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTg                                                          <  2:1058210/150‑1 (MQ=255)                                                                                               gCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAcc                                                       <  2:1466421/150‑1 (MQ=255)                                                                                               gCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAcc                                                       <  2:1598867/150‑1 (MQ=255)                                                                                                gtgtgGCTGTTGTCCCGCGTTGTTGAAAACCACTTTTTGCTGGCAATTGCCGGGGTTGCAATCCTTGCCGGGCCGGCGGTAGTGGCTTTGCTTGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAccc                                                      <  2:1810/149‑1 (MQ=255)                                                                                                cTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATGCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGAGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAccc                                                      <  1:140446/150‑1 (MQ=255)                                                                                                cTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAccc                                                      <  1:945964/150‑1 (MQ=255)                                                                                                cTGTGGCTGATGTCCCGCGTTCTGGAAAACACCTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAccc                                                      <  1:1591950/150‑1 (MQ=255)                                                                                                      cTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGt                                                >  2:1237852/1‑150 (MQ=255)                                                                                                      cTGTTGTCCCGCGTTCTGGAAAACAACTTTTAGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACCTCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGt                                                >  1:1285493/1‑150 (MQ=255)                                                                                                       tgttgtCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTg                                               >  1:714166/1‑150 (MQ=255)                                                                                                        gttgtCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTgg                                              >  2:571232/1‑150 (MQ=255)                                                                                                          tgtCCCGCGTTCTGGTAAACCAGTTTTTGCTGGCAATTGCAGTGCTGGCAATCCTTGCCTGCCGGGAGGTAGTTTCTTTGCTGGCGATAGATTAGAAGTTACATCACCTTTCGTCTAAACATCGCGTACGCCGCTGACCCTGTtgtgtcg                                            <  2:1284339/150‑4 (MQ=255)                                                                                                           gtCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCg                                           <  2:1214449/150‑1 (MQ=255)                                                                                                             cccGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGtt                                         <  2:614448/150‑1 (MQ=255)                                                                                                              ccGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTa                                        >  1:350993/1‑150 (MQ=255)                                                                                                                 cgTTCTGGATAACAACTTTTTTCTGGCATTTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTATACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCa                                     <  1:704642/150‑1 (MQ=255)                                                                                                                   tgcTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCAcc                                   <  1:35193/148‑1 (MQ=255)                                                                                                                     cTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAg                                 >  2:477860/1‑150 (MQ=255)                                                                                                                      tGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTAACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATAACCAGt                                >  2:598983/1‑150 (MQ=255)                                                                                                                       ggAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGCCCCTGTCGTGGCCGTTATAACCAGTa                               >  2:166823/1‑150 (MQ=255)                                                                                                                        gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAg                              >  1:56906/1‑150 (MQ=255)                                                                                                                           aacaacTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCgg                           >  2:450306/1‑150 (MQ=255)                                                                                                                             caacTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGCCCCTGTCGTGGCCGTTATCACCAGTAGCGGcc                         >  2:422881/1‑150 (MQ=255)                                                                                                                               acTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCcac                       <  2:327447/150‑1 (MQ=255)                                                                                                                               acTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCcac                       >  1:406889/1‑150 (MQ=255)                                                                                                                               acTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCcac                       >  1:685222/1‑150 (MQ=255)                                                                                                                                       gtGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTAtt               <  2:833158/149‑1 (MQ=255)                                                                                                                                         ggCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTcc             >  1:56951/1‑150 (MQ=255)                                                                                                                                         ggCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACATTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTcc             <  2:1214851/150‑1 (MQ=255)                                                                                                                                            aaTTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcaca          <  1:1728932/150‑1 (MQ=255)                                                                                                                                             aTTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcacac         <  1:753662/150‑1 (MQ=255)                                                                                                                                              ccgccgggCCGGCAAGGATTGCCAGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcacaca        >  2:897147/6‑150 (MQ=255)                                                                                                                                               tGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCAACACaa       >  2:1449426/1‑150 (MQ=255)                                                                                                                                               cgccgggCCGGCAAGGATTGCCAGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACaa       >  1:1483644/5‑150 (MQ=255)                                                                                                                                                gCAGGGCTGGCAAGCCTAGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGCGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAat      <  1:103223/150‑1 (MQ=255)                                                                                                                                                   gggCCGGCAAGGATTGCCAGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGCTAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCCCAAACTATTCCACACCatatg   >  2:1373757/1‑149 (MQ=255)                                                                                                                                                    ggCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACGATTCCCCCCCatatgc  >  1:1573125/1‑148 (MQ=255)                                                                                                                                                       |                                                                                                                                                   GTTTCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCC  >  NC_000913/824800‑825095

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