Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	401,421	A→G	8.5%	Q336R (CAG→CGG)	PP_0334 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	401,421	0	A	G	8.5%	145.7 / 3.7	58	Q336R (CAG→CGG)	PP_0334	transposase
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base G (3/2);  total (31/27)
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

CGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAGGTAAGAC  >  NC_002947/401275‑401544                                                                                                                                                     |                                                                                                                            cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTTCAGAGCGATCAg                                                                                                                            <  1:88314/148‑1 (MQ=17)          ccAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAAc                                                                                                                  >  4:343011/1‑149 (MQ=32)             gCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGc                                                                                                                >  7:198875/1‑148 (MQ=32)                 aGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                                                                                                            <  1:218175/148‑1 (MQ=32)                     cTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGGGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                                                                                                       >  1:291803/1‑149 (MQ=17)                      tGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                                                                                                       >  6:97348/1‑148 (MQ=32)                       gATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGCGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACCGAGCGATCAGGTTTACCACCGCAAACcgcacg                                                                                                      >  7:171337/1‑148 (MQ=255)                        aTTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                                                                                                       <  5:97348/146‑1 (MQ=32)                            cGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAAc                                                                                                <  8:188152/149‑1 (MQ=32)                             gATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAAcc                                                                                               >  4:331229/1‑149 (MQ=32)                                 aGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                                            <  3:331229/148‑1 (MQ=32)                                  gCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAgg                                                                                          >  4:45515/1‑149 (MQ=32)                                   ccAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTAGAGAGCGATCAGGTTTACCAAGGCAAACCGAACGAAAAACCGGAggc                                                                                         >  3:23737/1‑149 (MQ=17)                                     aGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg                                                                                       >  1:89203/1‑149 (MQ=32)                                     aGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg                                                                                       >  8:290020/1‑149 (MQ=32)                                      gatgatCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCGGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcggc                                                                                      >  3:102063/1‑149 (MQ=17)                                        tgatCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTa                                                                                    <  8:235802/149‑1 (MQ=32)                                             gcggccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGc                                                                               <  7:70817/146‑1 (MQ=17)                                              ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGc                                                                               >  2:330253/1‑148 (MQ=32)                                               cgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATGCAGGCTGGCTAATCGAAATCTCCG‑GATATTCTGGGCCGTGTTAAAGAGGGAAAGGTTTTACCAACGCAACCCGAACGAAAAACCGGTGGCGCCTATGAGcg                                                                              >  1:281025/1‑148 (MQ=11)                                                         gcTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAAcc                                                                   <  2:291803/149‑1 (MQ=32)                                                             gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAAc                                                                                                >  4:188943/1‑116 (MQ=32)                                                             gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAAc                                                                                                <  3:188943/116‑1 (MQ=32)                                                                 tCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATg                                                           <  3:45515/149‑1 (MQ=32)                                                                 tCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTCCCTAATCGAAATCTGCG‑GATTATCTCGGCCCTGTAACAGAGCGATCATGTTTACCAACGCAAACCGAACGAAACACCGGAGGCGGCTATGAGCCTCGGACCAACCGGGTAATg                                                           >  8:196650/1‑149 (MQ=11)                                                                  cccGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGAG‑GATTATCTCGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACAGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAtt                                                           >  6:286105/1‑147 (MQ=12)                                                                   ccGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCGCGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGc                                                          >  1:96466/1‑149 (MQ=17)                                                                   ccGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGc                                                          <  3:343011/148‑1 (MQ=32)                                                                     gCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCAc                                                       >  6:218772/1‑149 (MQ=32)                                                                      ccTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCGGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCAcc                                                      >  7:47221/1‑149 (MQ=17)                                                                       cTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCg                                                     >  8:19709/1‑149 (MQ=32)                                                                           ggCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTCTGAGCCTCTGACCAACCGGGTAATGCCACCgggg                                                  <  6:332711/148‑1 (MQ=17)                                                                              gcAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCACCGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGgcgc                                               >  7:176112/1‑148 (MQ=17)                                                                                 aggaggCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCGGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAgg                                           <  8:47221/149‑1 (MQ=17)                                                                                      gcggcgAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTc                                      <  2:276355/149‑1 (MQ=32)                                                                                                 ggctggctAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGc                            <  2:89203/149‑1 (MQ=32)                                                                                                  gctggctAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCt                              >  4:231418/1‑146 (MQ=32)                                                                                                  gctggctAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCt                              <  3:231418/146‑1 (MQ=32)                                                                                                    tggctAATCGAAATCTGCG‑GATTATCTGGTCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAg                         <  2:376425/149‑1 (MQ=17)                                                                                                    tggctAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAg                         >  1:83334/1‑149 (MQ=32)                                                                                                      gctAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCACCCCCTGCTAGt                        >  7:346286/1‑148 (MQ=17)                                                                                                      gctAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGt                        >  5:67362/1‑148 (MQ=32)                                                                                                           tCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGTCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACAt                  <  5:218772/149‑1 (MQ=17)                                                                                                             gAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACAtt                 >  1:262822/1‑148 (MQ=32)                                                                                                              aaaTCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGa               >  3:19774/1‑149 (MQ=32)                                                                                                              aaaTCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGa               <  1:330253/149‑1 (MQ=32)                                                                                                               aagcTGCG‑GATTATCTGGGCCCTTGTACAGAGCGATCAGGCTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAc              <  8:176112/146‑1 (MQ=255)                                                                                                               aaTCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAc              <  7:19709/149‑1 (MQ=32)                                                                                                               aaTCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAc              <  2:83334/149‑1 (MQ=32)                                                                                                                  cTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCt           <  2:262822/149‑1 (MQ=32)                                                                                                                     cG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAccc            <  7:290020/145‑1 (MQ=32)                                                                                                                     cG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAgg        >  1:301910/1‑149 (MQ=32)                                                                                                                      cgtaTTATCTGGGCCCTGTTATAGAGCGATCAGGTTTACCAACGCAAACCGAACGTAAAACCTGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAg         <  5:286105/145‑1 (MQ=255)                                                                                                                          ttATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAGGt       >  7:123377/1‑146 (MQ=32)                                                                                                                          ttATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAGGt       <  8:123377/146‑1 (MQ=32)                                                                                                                            aTCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAGGTAAGAc  >  5:117645/1‑149 (MQ=32)                                                                                                                             tCTGGGCCCTGTTACAGAGCGATCGGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGt                        >  2:94607/1‑126 (MQ=14)                                                                                                                             tCTGGGCCCTGTTACAGAGCGATCGGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGt                        <  1:94607/126‑1 (MQ=14)                                                                                                                                                     |                                                                                                                            CGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGAT‑CAGGCTGGCTAATCGAAATCTGCG‑GATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAGGTAAGAC  >  NC_002947/401275‑401544

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