Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,457,986	A→T	6.3%	T164S (ACC→TCC)	paaF →	putative 2,3‑dehydroadipyl‑CoA hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,457,986	0	A	T	6.3%	184.6 / 3.2	64	T164S (ACC→TCC)	paaF	putative 2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (39/21);  new base T (2/2);  total (41/23)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

GCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCG  >  NC_000913/1457837‑1458132                                                                                                                                                      |                                                                                                                                                   gCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCa                                                                                                                                                    >  1:1154423/1‑150 (MQ=255)   gATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCAcc                                                                                                                                                  >  2:1012651/1‑150 (MQ=255)      gtggtgGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCt                                                                                                                                               <  2:69466/150‑1 (MQ=255)      gtggtgGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCt                                                                                                                                               >  2:79098/1‑150 (MQ=255)       tggtggTTGCCGGATAAAACGCTAGTATTGTGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTc                                                                                                                                              >  2:655237/1‑150 (MQ=255)        ggtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCa                                                                                                                                             <  2:1506486/150‑1 (MQ=255)        ggtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCa                                                                                                                                             <  1:1012651/150‑1 (MQ=255)        ggtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCa                                                                                                                                             >  1:1223903/1‑150 (MQ=255)         gtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCag                                                                                                                                            <  2:1154423/150‑1 (MQ=255)             ttGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGCACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaag                                                                                                                                        >  1:1333247/1‑150 (MQ=255)              tGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagc                                                                                                                                       >  1:1482799/1‑150 (MQ=255)              tGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagc                                                                                                                                       <  1:1420694/150‑1 (MQ=255)              tGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagc                                                                                                                                       >  2:1420694/1‑150 (MQ=255)              tGCCGGAGAGAACGCGCGTTTTGAGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagc                                                                                                                                       <  1:1104347/150‑1 (MQ=255)               gCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagca                                                                                                                                      >  2:267110/1‑150 (MQ=255)                ccGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGcac                                                                                                                                     >  2:463084/1‑150 (MQ=255)                   gagagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagc                                                                                                                                  <  2:714259/150‑1 (MQ=255)                    agagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCagc                                                                                                                                           >  1:406668/1‑140 (MQ=255)                    agagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCagc                                                                                                                                           <  2:406668/140‑1 (MQ=255)                          cgcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCgg                                                                                                                           >  2:1183278/1‑150 (MQ=255)                                     ggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCg                                                                                                                >  2:1393712/1‑150 (MQ=255)                                        tGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACg                                                                                                             <  2:1223903/150‑1 (MQ=255)                                                aTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTcccc                                                                                                     >  2:938250/1‑150 (MQ=255)                                                    ctctCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCg                                                                                                 >  2:25002/1‑150 (MQ=255)                                                      ctcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGATAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAt                                                                                               >  2:401929/1‑150 (MQ=255)                                                      ctcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAt                                                                                               >  2:137667/1‑150 (MQ=255)                                                          gCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTaa                                                                                           <  1:938250/150‑1 (MQ=255)                                                                tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                                     <  1:1233578/150‑1 (MQ=255)                                                                tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                                     <  2:402513/150‑1 (MQ=255)                                                                tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                                     >  2:160772/1‑150 (MQ=255)                                                                  ccTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGaa                                                                                   <  1:1279372/150‑1 (MQ=255)                                                                   cTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAAt                                                                                  <  1:79098/150‑1 (MQ=255)                                                                     ggCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATAc                                                                                >  1:1061844/1‑150 (MQ=255)                                                                           ggCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGAAGTCTTCCCCAGCGATTTACCCCTCGAATACGACTTa                                                                          >  2:1001241/1‑150 (MQ=255)                                                                            gCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTAc                                                                         <  1:655237/150‑1 (MQ=255)                                                                            gCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTAc                                                                         <  1:467106/150‑1 (MQ=255)                                                                                aaCGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCt                                                                     >  1:1218691/1‑150 (MQ=255)                                                                                 aCGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTg                                                                    <  1:991927/150‑1 (MQ=255)                                                                                 aCGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTg                                                                    >  2:485051/1‑150 (MQ=255)                                                                                  cGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTgg                                                                   <  1:1405046/150‑1 (MQ=255)                                                                                    cAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCa                                                                 >  2:1386368/1‑150 (MQ=255)                                                                                      aCGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATc                                                               >  1:513381/1‑150 (MQ=255)                                                                                       cGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCg                                                              >  2:941498/1‑150 (MQ=255)                                                                                       cGTTTAATCCGTAGTGTCGGAAAATAGTTAGCCAGCAAAATGGTCCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGTTTAGCGACGTCTTCCCCAGCGAGTTAACCCTCGAATACGCCTTACAGCTGGCATCg                                                              >  2:782654/1‑150 (MQ=255)                                                                                        gTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGa                                                             >  2:219005/1‑150 (MQ=255)                                                                                        gTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGa                                                             >  2:543644/1‑150 (MQ=255)                                                                                           tAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGaaaa                                                          >  1:1417805/1‑150 (MQ=255)                                                                                            aaTCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAAt                                                         >  2:928104/1‑150 (MQ=255)                                                                                                         ggTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcg                                            >  2:1264449/1‑150 (MQ=255)                                                                                                           tAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcgcc                                          >  1:1463989/1‑150 (MQ=255)                                                                                                              aTCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCt                                       >  1:1178823/1‑150 (MQ=255)                                                                                                              aTCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCt                                       >  2:447310/1‑150 (MQ=255)                                                                                                               tCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAAt                                                                                  >  1:384039/1‑106 (MQ=255)                                                                                                               tCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAAt                                                                                  <  2:384039/106‑1 (MQ=255)                                                                                                                  ttAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGcc                                   <  2:1417805/150‑1 (MQ=255)                                                                                                                   tatccaccaAAATGGTGCTGACCGGACAAACTATCACCGCCCATCAACCACACCCCGCCGGGCTGGTTAGCGAAGGCCTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCt                                  <  1:1184246/147‑1 (MQ=255)                                                                                                                       cAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACa                              >  1:884267/1‑150 (MQ=255)                                                                                                                       cAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACa                              <  1:941498/150‑1 (MQ=255)                                                                                                                               tGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCaa                      >  1:960388/1‑150 (MQ=255)                                                                                                                                    cTGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa                 >  1:73049/1‑150 (MQ=255)                                                                                                                                      gAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaagc               >  1:1306587/1‑150 (MQ=255)                                                                                                                                            agaAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTgcg         >  2:881977/1‑150 (MQ=255)                                                                                                                                            agaAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTgcg         <  2:1061844/150‑1 (MQ=255)                                                                                                                                                aGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAg     >  1:469408/1‑150 (MQ=255)                                                                                                                                                   aTCACCGCTCAGCAAGCACCGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCCGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGCAAATGGCACGTCCCTCGCCGCTGGCCTTACAAGCGGCAAAGCCAGCGCTGCGCCAGTCg  >  2:330238/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCG  >  NC_000913/1457837‑1458132

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑