Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,579,709	G→A	7.1%	R359H (CGC→CAC)	PP_2260 →	glycerol‑phosphate ABC transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,579,709	0	G	A	7.1%	156.5 / 3.2	57	R359H (CGC→CAC)	PP_2260	glycerol‑phosphate ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base A (2/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

ATGTGCGCAACGAATACTGGCGCATGGTGCTGCACCTTCCCGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTGGGAACAGGGCGGTGCCTATGCGTTGCTCGGGCCGTC  >  NC_002947/2579562‑2579858                                                                                                                                                     |                                                                                                                                                      acgtgcgCAACGAATACTGGCGCATGGTGCTGCACCTTCCCGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGcc                                                                                                                                                     <  2:906856‑M2/148‑3 (MQ=255)   gtgCGCAACGAATACTGGCGCATGGTGCTGCACCTTCCCGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAg                                                                                                                                                   <  1:794478/150‑1 (MQ=255)       gcAACGAATACTGGCGCATGGTGCTGCACCTTCCCGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCga                                                                                                                                               <  1:1066202/150‑1 (MQ=255)              aTACTGGCGCATGGTGCTGCACCTTCCCGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAgg                                                                                                                                        <  1:654583/150‑1 (MQ=255)                   ggCGCATGGTGCTGCACCTTCCCGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGAt                                                                                                                                   <  1:358513/150‑1 (MQ=255)                      gcatcttgctgcACCTTCCCGGGGTTCACGAGTACCAGGTCGATACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCACGCCCCGGGGC‑TGCGCCAGGCGAGTAGAGGGTGATGGc                                                                                                                                <  2:156969/144‑1 (MQ=255)                              ctgcACCTTCCCGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGTTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATcc                                                                                                                        <  1:1461613/150‑1 (MQ=255)                              ctgcACCTTCCCGGGGTTCACGAGTACCAGCTGGACACGCCGATACGCGTGCTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCTCCAGGCGAGGAGAGGGTGATGGCCGAAATcc                                                                                                                        >  2:114740/1‑150 (MQ=255)                                         cGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAgg                                                                                                                                        <  2:103893/123‑1 (MQ=255)                                         cGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAgg                                                                                                                                        >  1:103893/1‑123 (MQ=255)                                              ttCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGaaa                                                                                                                           <  2:1223044/131‑1 (MQ=255)                                              ttCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGaaa                                                                                                                           >  1:1223044/1‑131 (MQ=255)                                                   gAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGG‑CCCGGGGC‑TGCGCCAGGCGAGGagag                                                                                                                                         <  2:329607/111‑1 (MQ=255)                                                   gAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGG‑CCCGGGGC‑TGCGCCAGGCGAGGagag                                                                                                                                         >  1:329607/1‑111 (MQ=255)                                                      tACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGtt                                                                                                >  2:1308939/1‑150 (MQ=255)                                                           ggTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGcc                                                                                                                               >  1:726172/1‑114 (MQ=255)                                                           ggTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGcc                                                                                                                               <  2:726172/114‑1 (MQ=255)                                                           ggTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAgc                                                                                           <  1:1420463/150‑1 (MQ=255)                                                            gTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAgcc                                                                                          >  2:134661/1‑150 (MQ=255)                                                                    gCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGccgccagccc                                                                                  >  1:901827/1‑150 (MQ=255)                                                                            gcgTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGa                                                                          >  2:584593/1‑150 (MQ=255)                                                                            gcgTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGa                                                                          >  1:675480/1‑150 (MQ=255)                                                                            gcgTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGCGCAG‑CCCCGGGGCCTGCAccaacagcccggcgctgtcgaaaa                                                                                                                               <  1:932346‑M2/97‑25 (MQ=255)                                                                            gcgTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGCGCAG‑CCCCGGGGCCTGCAccaacagcccggcgctgtcgaaaa                                                                                                                               >  2:932346‑M2/1‑73 (MQ=255)                                                                             cgTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAg                                                                         <  2:737995/150‑1 (MQ=255)                                                                                gtTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCa                                                                      <  1:1615443/150‑1 (MQ=255)                                                                                gtTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCa                                                                      >  1:1522459/1‑150 (MQ=255)                                                                                 ttCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAg                                                                     >  1:24356/1‑150 (MQ=255)                                                                                  tCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGa                                                                    <  1:584593/150‑1 (MQ=255)                                                                                    aTCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAggg                                                                                                                                       <  2:387286/81‑1 (MQ=255)                                                                                    aTCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAggg                                                                                                                                       >  1:387286/1‑81 (MQ=255)                                                                                       cccACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCGGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCCCCAGCCCGCCAGCGAGGCAGACTATg                                                               >  2:304438/1‑150 (MQ=255)                                                                                           cGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCt                                                           <  2:1315819/150‑1 (MQ=255)                                                                                             cacaAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTgc                                                         <  1:1308939/150‑1 (MQ=255)                                                                                              acaAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACa                                                                                             <  2:641391/113‑1 (MQ=255)                                                                                              acaAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACa                                                                                             >  1:641391/1‑113 (MQ=255)                                                                                              acaAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGCGCAG‑CCCCGGGGCCTGCAccaacagcccgg                                                                                                                                           <  1:827363‑M2/67‑13 (MQ=255)                                                                                              acaAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGCGCAG‑CCCCGGGGCCTGCAccaacagcccgg                                                                                                                                           >  2:827363‑M2/1‑55 (MQ=255)                                                                                               caAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTgcac                                                       >  1:119762/1‑150 (MQ=255)                                                                                                aaGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTgcacg                                                      >  1:871996/1‑150 (MQ=255)                                                                                                      ttCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCg                                                >  2:980704/1‑150 (MQ=255)                                                                                                         gTTTTCGACAGCGACGGGCTGTTGGTGCAGGCCCCGGGGG‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGAATGCGCCAGCTGGCGCACAGTTACAGCAGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCAAGCACTCGAAc                                             >  2:478972/1‑150 (MQ=255)                                                                                                           tttCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACAt                                           <  2:119762/150‑1 (MQ=255)                                                                                                               gACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACa                                                                                             <  1:952365/96‑1 (MQ=255)                                                                                                               gACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCATGCGCCAGCTGGCGCACAGTTACa                                                                                             >  2:952365/1‑96 (MQ=255)                                                                                                                aCAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACAtgtgtg                                      >  1:59736/1‑150 (MQ=255)                                                                                                                 cAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTgg                                     <  2:675480/150‑1 (MQ=255)                                                                                                                  aGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGcag                                                                                                    >  2:74575/1‑85 (MQ=255)                                                                                                                  aGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGcag                                                                                                    <  1:74575/86‑2 (MQ=255)                                                                                                                         ggCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGcaggcg                                                                                                 <  1:46457/82‑5 (MQ=255)                                                                                                                         ggCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGcaggcg                                                                                                 >  2:46457/1‑78 (MQ=255)                                                                                                                               tGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTGGGAACAGGGCGGTGc                       <  1:906856/150‑1 (MQ=255)                                                                                                                                   gCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTGGGAACAGGGCGGTGCCTAt                   <  2:1522459/150‑1 (MQ=255)                                                                                                                                      ggCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTGGGAACAGGGCGGTGCCTATGCg                <  1:1125435/150‑1 (MQ=255)                                                                                                                                       gCCCCGGGGC‑TGCGCCAGGCGAGGCGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGGTACAGCCGGCAGCCCACCAGCTAGGCGGACTATGCCCTGCAAGCACTCGAACATGTGTGGGAACAGGGCGGTGCATATGCGt               >  2:546226/1‑150 (MQ=255)                                                                                                                                          ccGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTGGGAACAGGGCGGTGCCTATGCGTTGc            >  2:1645706/1‑150 (MQ=255)                                                                                                                                                c‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTGGGAACAGGGCGGTGCCTATGCGTTGCTCGGGc      <  1:980704/150‑1 (MQ=255)                                                                                                                                                     gcCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTGGGAACAGGGCGGTGCCTAt                   <  2:564186/133‑1 (MQ=255)                                                                                                                                                     gcCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTGGGAACAGGGCGGTGCCTAt                   >  1:564186/1‑133 (MQ=255)                                                                                                                                                     gcCAGGCCATGAGAGGGTGATGGCCAAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCCCTCGAACATGTGTGGGAACAGGGCGGTGCCTATGCGTTGCTCGGGCCGTc  <  1:478972/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ATGTGCGCAACGAATACTGGCGCATGGTGCTGCACCTTCCCGGGGTTCACGAGTACCAGGTGGACACGCCGATACGCGTGTTCATCCCCACGCACAAGCTGTTCGTTTTCGACAGCGCCGGGCTGTTGGTGCAGGCCCCGGGGC‑TGCGCCAGGCGAGGAGAGGGTGATGGCCGAAATCCGCCTGCGCCAGCTGGCGCACAGTTACAGCCGCCAGCCCGCCAGCGAGGCAGACTATGCCCTGCACGCACTCGAACATGTGTGGGAACAGGGCGGTGCCTATGCGTTGCTCGGGCCGTC  >  NC_002947/2579562‑2579858

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑