Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	824,956	C→G	5.8%	I109M (ATC→ATG) ‡	ybhQ →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	824,956	0	C	G	5.8%	311.4 / 3.2	93	I109M (ATC→ATG) ‡	ybhQ	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (46/41);  new base G (3/3);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCA  >  NC_000913/824809‑825104                                                                                                                                                     |                                                                                                                                                     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:1732139/150‑1 (MQ=255)                    cGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAgtg                                                                                                                                 <  2:1440431/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:827687/1‑150 (MQ=255)                                            gagCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaagaag                                                                                                         >  1:614448/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:685222/1‑150 (MQ=255)                                                                                                                      acTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCcac                                >  1:406889/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)                                                                                                                                               ggCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGc                                   >  1:316521/1‑122 (MQ=255)                                                                                                                                               ggCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGc                                   <  2:316521/122‑1 (MQ=255)                                                                                                                                                gCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACt      <  2:171188/150‑1 (MQ=255)                                                                                                                                                 cAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTc     >  1:946011/1‑150 (MQ=255)                                                                                                                                                  aaTCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCg    <  1:820143/150‑1 (MQ=255)                                                                                                                                                    tCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCa  >  2:971306/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTA‑TTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCA  >  NC_000913/824809‑825104

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