Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,850,729	C→T	6.8%	Y478Y (TAC→TAT)	ilvD →	dihydroxy‑acid dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,850,729	0	C	T	6.8%	173.2 / 2.3	59	Y478Y (TAC→TAT)	ilvD	dihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base T (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CTGCGTGGTCAAAACCGCCGGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGCGGTACCTCGGGCCTGTCGATCGGC  >  NC_002947/5850588‑5850867                                                                                                                                               |                                                                                                                                           cTGCGTGGTCAAAACCGCCGGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGcc                                                                                                                                     >  6:119417/1‑149 (MQ=255) cTGCGGGGGCAAAACCGCCGGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGcc                                                                                                                                     <  8:293906/149‑1 (MQ=255)      tGGTCAAAACCGCCGGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGa                                                                                                                                  >  4:89772/1‑147 (MQ=255)       ggTCAAAACCGCCGGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAg                                                                                                                               >  6:333688/1‑149 (MQ=255)         tCAAAACCGCCGGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGAGTCTGCGAAAGCAAGGACAGCGCCGTGCGCGGCAACATCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGATACGAAGGCCCGAAAgg                                                                                                                              >  1:272696/1‑149 (MQ=255)         tCAAAACCGCCGGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCCTCATCCGCTACGAAGGCCCGAAAggc                                                                                                                             >  6:133409/1‑149 (MQ=255)            aaaCCGCCGGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACAGCGTCATCATCCGCTACGAAGGACCGAAAggcggc                                                                                                                          >  7:335879/1‑149 (MQ=255)               ccgccgGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCg                                                                                                                       >  2:347062/1‑149 (MQ=255)                cgccgGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCgg                                                                                                                      >  5:67550/1‑149 (MQ=255)                 gccgGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCggg                                                                                                                     <  5:133409/149‑1 (MQ=255)                 gccgGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCggg                                                                                                                     <  1:314765/149‑1 (MQ=255)                  ccgGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGc                                                                                                                    <  7:151905/149‑1 (MQ=255)                   cgGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCa                                                                                                                   >  5:324358/1‑149 (MQ=255)                     gtgtCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATg                                                                                                                 <  7:33996/149‑1 (MQ=255)                     gtgtCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATg                                                                                                                 <  3:89772/149‑1 (MQ=255)                        gcgacgaATCGATCCACGTGGTCGAGGGCACCGCGAAGA‑TCTTCGAAACCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCaa                                                                                                              <  8:335879/148‑1 (MQ=255)                          gacgaATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGGCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCGCGAAAGGCGGCCCGGGCATGCAAg                                                                                                             <  6:112155/148‑1 (MQ=255)                          gacgaATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTATGAAGGCCCGAAAGGCGGCCCGGGCATGCAAg                                                                                                             >  8:47600/1‑148 (MQ=255)                           acgaATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGa                                                                                                            >  3:108804/1‑148 (MQ=255)                           acgaATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGa                                                                                                            >  5:322351/1‑148 (MQ=255)                              aaTCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATg                                                                                                        >  3:287607/1‑149 (MQ=255)                              aaTCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATg                                                                                                        <  6:67550/149‑1 (MQ=255)                              aaTCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATg                                                                                                        <  2:272696/149‑1 (MQ=255)                                        cGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGa                                                                                              >  2:145909/1‑149 (MQ=255)                                               gAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTa                                                                                       >  2:108405/1‑149 (MQ=255)                                                aGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTAc                                                                                      <  5:119417/149‑1 (MQ=255)                                                aGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTAc                                                                                      >  4:287863/1‑149 (MQ=255)                                                   gCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTc                                                                                   <  4:103115/149‑1 (MQ=255)                                                     aCCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTc                                                                                   >  2:220624/1‑147 (MQ=255)                                                     aCCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCaa                                                                                 <  3:287863/149‑1 (MQ=255)                                                        gcgAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTc                                                                              <  4:108804/149‑1 (MQ=255)                                                               ttCTTCGAAACCCAGGACTCCGCCGGTCGCGGCTTCCTCGCCCCCGAGGTTAAAGCCGGCGACATCGTCATCATCCGCTACTAAGGCCCGAAAGGCGGCCCGGGCATTCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAg                                                                         <  8:226533/147‑1 (MQ=255)                                                                 cTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTATGAAGGCCCGAAAggcg                                                                                                                            >  5:243606/1‑95 (MQ=255)                                                                 cTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTATGAAGGCCCGAAAggcg                                                                                                                            <  6:243606/95‑1 (MQ=255)                                                                  ttCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCt                                                                     >  1:335541/1‑149 (MQ=255)                                                                    cGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCt                                                                     <  6:324358/147‑1 (MQ=255)                                                                    cGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCt                                                                     <  1:198917/147‑1 (MQ=255)                                                                    cGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTgg                                                                   >  7:223957/1‑149 (MQ=255)                                                                     gAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAATGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTgg                                                                   <  8:223957/148‑1 (MQ=255)                                                                        aGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCaa                                                               <  1:145909/149‑1 (MQ=255)                                                                             ggACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAgg                                                             <  4:206067/146‑1 (MQ=255)                                                                             ggACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAgg                                                             >  3:206067/1‑146 (MQ=255)                                                                                 aGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTgcgc                                                      <  5:333688/149‑1 (MQ=255)                                                                                     ccGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTg                                                  <  1:220624/149‑1 (MQ=255)                                                                                             ggCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGa                                          >  3:170493/1‑149 (MQ=255)                                                                                               cATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCAGCCCGGGCATGCAAg                                                                                                             <  7:97893/80‑1 (MQ=255)                                                                                               cATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCAGCCCGGGCATGCAAg                                                                                                             >  8:97893/1‑80 (MQ=255)                                                                                                   cTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGc                                                                                                       <  4:188500/82‑1 (MQ=255)                                                                                                   cTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGc                                                                                                       >  3:188500/1‑82 (MQ=255)                                                                                                             aGGTGAAAGCCGGCGACATCGTCATCATCCGCTATGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGc                          <  7:47600/149‑1 (MQ=255)                                                                                                               gTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGc                                                                       <  7:294131/102‑1 (MQ=255)                                                                                                               gTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGc                                                                       >  8:294131/1‑102 (MQ=255)                                                                                                                 gAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCTTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGCGGTa                      >  8:184165/1‑149 (MQ=255)                                                                                                                    aGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGCGGTACCt                   <  2:335541/149‑1 (MQ=255)                                                                                                                      ccGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGCGGTACCTCg                 <  4:287607/149‑1 (MQ=255)                                                                                                                      ccGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGCGGTACCTCg                 <  6:322351/149‑1 (MQ=255)                                                                                                                       cGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGCGGTACCTCgg                <  6:62975/149‑1 (MQ=255)                                                                                                                           gacatcgacatcATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGCGGTACCTCGGGCCt            <  6:174496/149‑1 (MQ=255)                                                                                                                                      catcCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGCGGTACCTCGGGCCTGTCGATCGGc  >  3:107323/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                           CTGCGTGGTCAAAACCGCCGGTGTCGACGAATCGATCCACGTGTTCGAGGGCACCGCGAAGA‑TCTTCGAAAGCCAGGACAGCGCCGTGCGCGGCATCCTCGCCGACGAGGTGAAAGCCGGCGACATCGTCATCATCCGCTACGAAGGCCCGAAAGGCGGCCCGGGCATGCAAGAAATGCTCTACCCGACCTCGTACCTCAAGTCCAAAGGCCTGGGCAAGGCCTGCGCCCTGCTCACCGACGGTCGCTTCTCGGGCGGTACCTCGGGCCTGTCGATCGGC  >  NC_002947/5850588‑5850867

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