Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	787,474	G→A	A41V (GCA→GTA)	gpmA ←	phosphoglyceromutase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	787,474	0	G	A	98.2%	203.9 / NA	56	A41V (GCA→GTA)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/22);  minor base C (0/1);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTCGTTCAGTTTCCAGGATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTGCTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATATACTTACTCCTCAAATCATCTTTT  >  NC_000913/787326‑787620                                                                                                                                                     |                                                                                                                                                   aCGTTCGTTCAGTTTCCAGGATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACctact                                                                                                                                                  >  2:781200/1‑151 (MQ=255)   gttcgttcAGTTTCCAGGATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCtactgc                                                                                                                                                >  1:691015/1‑151 (MQ=255)      cgttcAGTTTCCAGGATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCttt                                                                                                                                             >  1:152871/1‑151 (MQ=255)      cgttcAGTTTCCAGGATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCttt                                                                                                                                             <  2:805946/151‑1 (MQ=255)                ccAGGATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCtt                                                                                                                                   <  1:479337/151‑1 (MQ=255)                ccAGGATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCtt                                                                                                                                   <  2:525870/151‑1 (MQ=255)                   ggATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTAc                                                                                                                                 >  2:685938/1‑150 (MQ=255)                    gATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGc                                                                                                                               <  2:532954/151‑1 (MQ=255)                      tttCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCt                                                                                                                             >  2:439190/1‑151 (MQ=255)                       ttCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCtt                                                                                                                            >  1:89789/1‑151 (MQ=255)                        tctcAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCttt                                                                                                                           >  1:385266/1‑151 (MQ=255)                               gggCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTcagac                                                                                                                    >  1:164358/1‑151 (MQ=255)                                 gCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTcagac                                                                                                                    >  1:536797/1‑149 (MQ=255)                                       aTGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCAc                                                                                                            >  1:511421/1‑151 (MQ=255)                                       aTGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCAc                                                                                                            <  2:295120/151‑1 (MQ=255)                                         gCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCAcgt                                                                                                          >  1:683472/1‑151 (MQ=255)                                         gCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCAACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCAcgt                                                                                                          >  2:577645/1‑151 (MQ=255)                                            tGATACAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCAcgtcgt                                                                                                       >  1:774197/1‑151 (MQ=255)                                             gATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCAcgtcgt                                                                                                       <  1:523152/150‑1 (MQ=255)                                             gATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCAcgtcgt                                                                                                       <  1:829597/150‑1 (MQ=255)                                              aTCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTAc                                                                                                     >  1:194923/1‑151 (MQ=255)                                                ccAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTTTAAGCAAAGACAAAGCTGTAACCTTCCTCTTTCAGAAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTAcc                                                                                                    >  1:232654/1‑150 (MQ=255)                                                           aGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAAc                                                                                        >  1:825004/1‑151 (MQ=255)                                                                 ttCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGtttt                                                                                  >  1:738673/1‑151 (MQ=255)                                                                   ccACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCt                                                                                >  1:398240/1‑151 (MQ=255)                                                                    cacaGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCtt                                                                               >  2:683256/1‑151 (MQ=255)                                                                     acaGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCttt                                                                              >  2:659266/1‑151 (MQ=255)                                                                        gggTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGtt                                                                           <  2:540036/151‑1 (MQ=255)                                                                           tATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTcc                                                                         >  2:681059/1‑150 (MQ=255)                                                                             tGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCAc                                                                       >  2:759509/1‑150 (MQ=255)                                                                              ggATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACt                                                                      <  2:645995/150‑1 (MQ=255)                                                                                 tAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACt                                                                  >  2:718218/1‑151 (MQ=255)                                                                                   gcgcGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACttt                                                                <  1:819324/151‑1 (MQ=255)                                                                                           cAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTCCGGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAGCCGGTGAAACGGTTTTATTTGTTCCACTGACTTTCGCCAt                                                          <  1:604207/149‑1 (MQ=255)                                                                                            aGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGAc                                                       <  1:720634/151‑1 (MQ=255)                                                                                              cacaGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACg                                                      >  2:550219/1‑150 (MQ=255)                                                                                               acaGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACga                                                     >  2:488922/1‑150 (MQ=255)                                                                                                   aaGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACgaaccag                                                <  2:548586/151‑1 (MQ=255)                                                                                                     gtgtAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACgaaccagaa                                              <  2:541449/151‑1 (MQ=255)                                                                                                     gtgtAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACgaaccaga                                               >  2:522036/1‑150 (MQ=255)                                                                                                            cAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCtt                                       <  2:330455/151‑1 (MQ=255)                                                                                                             agaGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCtt                                       >  2:324021/3‑150 (MQ=255)                                                                                                             aaaGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTa                                      <  1:292704/151‑1 (MQ=255)                                                                                                              aaGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAg                                     <  2:683515/151‑1 (MQ=255)                                                                                                                gTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGtt                                   >  2:131820/1‑151 (MQ=255)                                                                                                                  cAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTAc                                 >  1:138219/1‑151 (MQ=255)                                                                                                                  cAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTAc                                 <  2:272311/151‑1 (MQ=255)                                                                                                                     aGCTGTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGc                              >  2:402198/1‑151 (MQ=255)                                                                                                                         gTAACCTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCat                           >  2:797740/1‑150 (MQ=255)                                                                                                                             ccTTCCTCTTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATAtact                      >  2:423284/1‑151 (MQ=255)                                                                                                                                    ctTTCAGCAGCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATATACTTActcct               >  2:850094/1‑151 (MQ=255)                                                                                                                                        cagcagCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATATACTTACTCCTCaa            <  2:279407/150‑1 (MQ=255)                                                                                                                                        cagcagCTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATATACTTACTCCCCaa            <  2:396914/150‑1 (MQ=255)                                                                                                                                            agcTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATATACTTACTCCTCAAatcat       <  1:96863/151‑1 (MQ=255)                                                                                                                                            agcTTACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATATACTTACTCCTCAAatcat       >  1:805368/1‑151 (MQ=255)                                                                                                                                             gcTTCCCTACTGCTTTTGCTTCGCTTACGCCTTTCCCAGACAGATCCCCGGCGTACCAGCCGGTGCAACTGTGTTCTTTGTTCGCCTGCCTTTCGCCATTACGAACGAGAACCAGCTTAGTTACAGCCATATCCTTACTCCTCAAatcat       <  1:611908/150‑1 (MQ=255)                                                                                                                                               ttACCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATATACTTACTCCTCAAATCATCtt    >  2:42501/1‑151 (MQ=255)                                                                                                                                                 aCCTACTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATATACTTACTCCTCAAATCATCtttt  <  2:742505/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ACGTTCGTTCAGTTTCCAGGATTTCTCAACGGGCAGCCATGCCTGATCCAGTTCGTCCAGCACATTCCACAGGGTATGGATAGCGCGTTTCAGCACAGAAGTGTAAGCAAAGTCAAAGCTGTAACCTTCCTCTTTCAGCAGCTTACCTGCTGCTTTTGCTTCGCTTACGCCTTTCTCAGACAGATCCACGTCGTACCAACCGGTGAAACGGTTTTCTTTGTTCCACTGACTTTCGCCATGACGAACCAGAACCAGCTTAGTTACAGCCATATACTTACTCCTCAAATCATCTTTT  >  NC_000913/787326‑787620

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑