Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,609,561:1	+C	6.3%	coding (1120/1452 nt)	uxaB ←	altronate oxidoreductase, NAD‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,609,561	1	.	C	6.3%	304.2 / 16.6	79	coding (1120/1452 nt)	uxaB	altronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base . (30/44);  new base C (2/3);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTT  >  NC_000913/1609414‑1609710                                                                                                                                                     |                                                                                                                                                      gCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                                                                                     <  1:1181260/150‑1 (MQ=255)   tGATAACGTTCCAGCCAGTGTGCATCGTCCTGTAGCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCTCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCATTCTGACCTTCCAGCAGCTG‑TGGc                                                                                                                                                   >  2:917878/1‑150 (MQ=255)    gATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCa                                                                                                                                                  <  2:610568/150‑1 (MQ=255)    gATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCa                                                                                                                                                  >  1:159878/1‑150 (MQ=255)    gATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCa                                                                                                                                                  <  1:522925/150‑1 (MQ=255)      tAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAgg                                                                                                                                                >  2:811695/1‑150 (MQ=255)      tAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgg                                                                                                                                                <  1:418909/150‑1 (MQ=255)       aaCGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                                                                               <  2:670759/150‑2 (MQ=255)        aCGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGAt                                                                                                                                              >  1:1097452/1‑150 (MQ=255)        aCGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggc                                                                                                                                              <  1:645432/150‑3 (MQ=255)         cGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACCCTCACCACGATAGAACGCAATCAATGCCGCTAATCCGAAAGTCAGGCGCGCCGGAAGTGAGCCGTTTGCCTTCTGCCCTGCCAGCAGCTC‑TGGCAGGATg                                                                                                                                             >  1:339784/1‑150 (MQ=255)           ttCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCTGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCg                                                                                                                                           <  2:815569/150‑1 (MQ=255)            tCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGTCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCgg                                                                                                                                          <  2:106055/150‑1 (MQ=255)              cagccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATCCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGAGTGCCGTTTGCCTTCTGCCCAGCCAGCAGCTG‑TGGCAGGATGCGGGt                                                                                                                                        >  1:623050/1‑150 (MQ=255)                   agTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGGGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGa                                                                                                                                   <  1:342762/150‑1 (MQ=255)                      gtgCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACt                                                                                                                                >  1:1252771/1‑150 (MQ=255)                      gtgCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACt                                                                                                                                <  2:205555/150‑1 (MQ=255)                         cATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                                                                             <  1:74159/150‑1 (MQ=255)                         cATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                                                                             <  1:1466147/150‑1 (MQ=255)                         cATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                                                                             <  2:1475688/150‑1 (MQ=255)                           tCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTc                                                                                                                           <  2:515825/150‑1 (MQ=255)                              tCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGTTAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGCTGCGGGTGCGGACCTTGGTCATa                                                                                                                        >  2:1171395/1‑150 (MQ=255)                               ccTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGGGCCGTTTGCCTTCTGCCCTGCCAGTAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAc                                                                                                                       <  2:792852/150‑1 (MQ=255)                                cTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCACTCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                                                                                      <  2:1183071/150‑1 (MQ=255)                                cTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                                                                                      <  2:1356557/150‑1 (MQ=255)                                      cGGATATGTTGCGCCGTTACGCTCACCCCGATCGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGGCATACCGTTGAg                                                                                                                <  1:1565122/150‑1 (MQ=255)                                       ggATATGTTTCGCCGTTACGCTCACCACGATAGAACTCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTCCGGAACTTGGTCATACCGTTGAgc                                                                                                               >  2:260441/1‑150 (MQ=255)                                         atatGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcgc                                                                                                             <  1:52742/150‑1 (MQ=255)                                          tatGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCTCTAAGGCGAACGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                                                                            <  1:340940/150‑1 (MQ=255)                                          tatGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                                                                            <  1:60131/150‑1 (MQ=255)                                             gTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATa                                                                                                         <  1:1004941/150‑1 (MQ=255)                                             gTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATa                                                                                                         <  2:1252771/150‑1 (MQ=255)                                              tttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGAGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                                                                        >  1:1148426/1‑150 (MQ=255)                                                 cGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                                                                                     <  2:154798/150‑1 (MQ=255)                                                          gCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGGGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGAt                                                                                            <  2:159878/150‑1 (MQ=255)                                                           cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGGGCCGGAAGTGTGGCGTTTGCCTTCTGCCCTTGCAGCAGGTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                                                                           <  1:917878/150‑1 (MQ=255)                                                               ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTa                                                                                       >  1:1385735/1‑150 (MQ=255)                                                               ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGACGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGGGCGGAACTTGGTCATACCGTTGAGCGCGATAGACCGCAGCTGCTGCTTa                                                                                       >  2:855344/1‑150 (MQ=255)                                                                cacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                                                                                      >  2:467289/1‑150 (MQ=255)                                                                  cGATCGAACGCCATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATg                                                                                    <  1:1572033/150‑1 (MQ=255)                                                                  cGATAGATCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCTGCAGCTG‑TGGGAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATg                                                                                    <  2:339784/150‑1 (MQ=255)                                                                       gAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACgg                                                                               >  2:1189885/1‑150 (MQ=255)                                                                       gAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACgg                                                                               >  2:1183650/1‑150 (MQ=255)                                                                        aaCGCTATCTATGCCGCTAATGCGAAAGTAAGGCGCGCCGGACGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTT‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACggg                                                                              <  2:1315553/150‑1 (MQ=255)                                                                             aatcactGCCGCTAATGCGATAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                                                                         <  2:39487/150‑1 (MQ=255)                                                                                caatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaa                                                                      <  2:320579/150‑1 (MQ=255)                                                                                   tGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg                                                                   >  2:1316999/1‑150 (MQ=255)                                                                                   tGCCGCTAATGCGAAAGTAAGGCGCGCAGGAAGTGTGCCGTTTGCCTTCTGCCGTGGCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAATCg                                                                   <  2:623050/150‑1 (MQ=255)                                                                                        cTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCgg                                                              >  1:1009268/1‑150 (MQ=255)                                                                                         tAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGt                                                             >  2:938099/1‑150 (MQ=255)                                                                                              cGAAAGTAAGGCGCGCCGGAAGGGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGTCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTg                                                        <  1:1223212/150‑1 (MQ=255)                                                                                                aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCa                                                      >  2:796404/1‑150 (MQ=255)                                                                                                 aaGTCAGGCGCGCCGGAAGGGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATCGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCAc                                                     <  1:224683/150‑1 (MQ=255)                                                                                                      aGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAACAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCg                                                >  2:62210/1‑150 (MQ=255)                                                                                                          gcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAg                                            <  1:12508/150‑1 (MQ=255)                                                                                                          gcgcCGGAAGGGTGCCGTTTGCCTTCTGCCCTGCCCGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAg                                            <  1:1471694/150‑1 (MQ=255)                                                                                                             ccGGAAGTGTGCCGTTTGCCTTCTGCCCAGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCCGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGAtt                                         >  1:704411/1‑150 (MQ=255)                                                                                                               ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTcc                                       <  1:1700234/150‑1 (MQ=255)                                                                                                               ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCGG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACTCCCGGTGACTGCACTGGCGAAAGATTcc                                       <  1:383867/150‑1 (MQ=255)                                                                                                                   gtgtGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGtt                                   >  1:1434870/1‑150 (MQ=255)                                                                                                                       gCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTtcatc                               <  2:1210195/150‑1 (MQ=255)                                                                                                                        ccGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTtcatca                              <  2:1097452/150‑1 (MQ=255)                                                                                                                        ccGTTTGCCTTCTGCCCTGCCAGCAGCGG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTtcatca                              <  1:169805/150‑1 (MQ=255)                                                                                                                             tGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAgg                         <  2:1385735/150‑1 (MQ=255)                                                                                                                              gCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGc                        >  1:982000/1‑150 (MQ=255)                                                                                                                               ccTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCa                       >  2:1433717/1‑150 (MQ=255)                                                                                                                               ccTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGGGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGGTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCa                       <  1:854245/150‑1 (MQ=255)                                                                                                                                cTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaa                      >  2:975619/1‑150 (MQ=255)                                                                                                                                cTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaa                      >  2:1668420/1‑150 (MQ=255)                                                                                                                                 ttCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaaa                     >  2:1569987/1‑150 (MQ=255)                                                                                                                                 ccgcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaaa                     >  1:1434328/7‑150 (MQ=255)                                                                                                                                    tgccctgccAGCAGCTG‑GGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATcc                  <  1:1555688/150‑1 (MQ=255)                                                                                                                                       cctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGt               >  1:583347/1‑150 (MQ=255)                                                                                                                                        ctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTa              >  2:495615/1‑150 (MQ=255)                                                                                                                                          gccAGCAGCGG‑GTGCAGGATGCGGGTGCGGTACTTGGTCATACCGTTGAGCGCTAAAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAcc            <  2:618737/150‑1 (MQ=255)                                                                                                                                           ccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAc             >  2:1055454/1‑148 (MQ=255)                                                                                                                                           ccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAc             <  1:1055454/148‑1 (MQ=255)                                                                                                                                            cagcagCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCgg          <  2:1504311/150‑1 (MQ=255)                                                                                                                                            cagcagCTG‑TGGCAGGATGCGGGGTCGTAACTTGGTCATACGGTTGAGCGCGATAGACATCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCTGTTCATCACGAGGCAAATCCAGTACCgg          <  2:357186/150‑1 (MQ=255)                                                                                                                                              gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaa        >  2:1112080/1‑150 (MQ=255)                                                                                                                                                   tg‑tgGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAtt   <  1:467289/150‑1 (MQ=255)                                                                                                                                                    g‑tgGCAGGACGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAttt  >  2:381860/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTT  >  NC_000913/1609414‑1609710

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