Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	101,803	A→T	10.3%	K208M (AAG→ATG) ‡	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	101,803	0	A	T	10.3%	146.7 / 7.5	58	K208M (AAG→ATG) ‡	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base A (25/26);  new base T (2/4);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

CTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTC  >  NC_002947/101664‑101948                                                                                                                                             |                                                                                                                                                   cTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTggg                                                                                                                                           <  2:52062/149‑2 (MQ=255)  tgaGCGAGCTTGGCTGCCGTTTCTCCAATCAGCTGGTGGACGCCACCGATTCCTGGAACAAGCTTGTTACCGACGAAGCTGCGCTGGCCGGACTGACCGATTCGGCCAAG‑GCGCAGATGGCAGCTGCAGTCCAGGACCTTGGCCTgggc                                                                                                                                          <  4:256847/147‑3 (MQ=255)  tCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAATCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGTAGACCATGCCAAGGGCCTGGAc                                                                                                                                          >  4:261991/1‑149 (MQ=255)  tCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAc                                                                                                                                          >  5:117497/1‑149 (MQ=255)  tCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAc                                                                                                                                          >  2:191507/1‑149 (MQ=255)    aGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTgggc                                                                                                                                          <  6:274974/147‑3 (MQ=255)    aGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTCGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCATATGGCCGCTGCAGCCCAGGCCAAGGGCCt                                                                                                                                              <  5:309965/143‑1 (MQ=255)    aGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTCGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCt                                                                                                                                              >  6:309965/1‑143 (MQ=255)        aGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTgg                                                                                                                                   <  8:28888/149‑1 (MQ=255)        aGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTgg                                                                                                                                   >  4:14834/1‑149 (MQ=255)         gcggggCAGCCGTTTCTCCAATCAGCGGCTGCACGCCACCCAGGCCTGCACCAAGCATGTTACCGACTATGCTGCGCTGGCCGGCCTGACCGTTTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGt                                                                                                                                  <  4:194829/146‑1 (MQ=255)             ggCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGc                                                                                                                                      >  2:113719/1‑141 (MQ=255)             ggCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACCAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGc                                                                                                                                      <  1:113719/141‑1 (MQ=255)                 gCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGTCCAGGCCCTTGGCCTgggctg                                                                                                                                        <  5:172975/136‑5 (MQ=255)                 gCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGTCCAGGCCCTTGGCCTgggctg                                                                                                                                        >  6:172975/1‑132 (MQ=255)                 gCCGTTTCTCCAATCAGCTGCTGGACGCCAACCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCAcc                                                                                                                          <  7:177791/149‑1 (MQ=255)                  ccGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGTCCAAGCATGTTACCGACGAAGGTGCGCTGGCCGGCCTGACCGATTCGGCCTAGTGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTAGTTGATCa                                                                                                                            <  1:89107/147‑1 (MQ=255)                  ccGTTTCTACAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGCCCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑TCCTAGATGATCa                                                                                                                            >  2:89107/1‑146 (MQ=255)                     tttCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCCGCCAAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTg                                                                                                                      <  6:157772/149‑1 (MQ=255)                       tctcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGa                                                                                                                     <  1:191507/148‑1 (MQ=255)                              tCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCa                                                                                                             <  6:117497/149‑1 (MQ=255)                              tCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCa                                                                                                             <  3:14834/149‑1 (MQ=255)                                cgctgcgGGCCGCCACCCCGGCCTGGACCACGCATGTTACCGACGCAGCGGCGCCGGCCGGCGTGACCGATTGGGCCAAG‑GCGCAGATGGCCGCTCCCGCCCCGGCCACGGCCCTGCAC‑GGCTGGTTGATCACCCTTGAATTCCCCa                                                                                                             <  2:331964/146‑1 (MQ=255)                                 gctgctGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAgg                                                                                                                                                  <  3:155019/110‑1 (MQ=255)                                 gctgctGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAgg                                                                                                                                                  >  4:155019/1‑110 (MQ=255)                                 gctgctGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCa                                                                                                             >  1:271461/1‑146 (MQ=255)                                 gctgcgGGACCCCACCCCGGCCTGGCCCAAGCATGTTACCTACTAAGCTGCGCTGGCCGGCCTTACCGATGCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCTGGCCAAGGGCGTGGAC‑GGCTGGTTTATCACCCTTTAATTCCCCa                                                                                                             <  2:271461/146‑1 (MQ=255)                                           gCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGt                                                                                                <  4:322236/149‑1 (MQ=255)                                              aCCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtgat                                                                                             >  6:357692/1‑149 (MQ=255)                                                ccAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtgatga                                                                                           >  2:29545/1‑149 (MQ=255)                                                    gCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGa                                                                                                                               >  4:314715/1‑109 (MQ=255)                                                    gCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGa                                                                                                                               <  3:314715/109‑1 (MQ=255)                                                         gACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGc                                                                                    >  3:237715/1‑147 (MQ=255)                                                           ccAAGCATGTTACCGACGAAGCTTCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGACCAGGCCAAGGGCCTGCAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGCTGATGACCTACGCCAGc                                                                                >  3:66559/1‑149 (MQ=255)                                                            cAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTg                                                                                                                                <  7:272524/100‑1 (MQ=255)                                                            cAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTg                                                                                                                                >  8:272524/1‑100 (MQ=255)                                                            cAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGCTCACCCTTGAATTCCCCAGCTACTACGCGGTGATGGCCTAGGCCAGCg                                                                               >  8:188554/1‑149 (MQ=255)                                                                 aTGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACcgc                                                                          <  1:58735/149‑1 (MQ=255)                                                                            gaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGACCACCCTTGAATTCCCCAGCTACTAcgcg                                                                                                  >  2:258349/1‑114 (MQ=255)                                                                            gaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGACCACCCTTGAATTCCCCAGCTACTAcgcg                                                                                                  <  1:258349/114‑1 (MQ=255)                                                                                    gcTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAagag                                                           <  8:321407/145‑1 (MQ=255)                                                                                    gcTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAagag                                                           >  7:321407/1‑145 (MQ=255)                                                                                       ggccggccTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATAcg                                                    >  2:100292/1‑149 (MQ=255)                                                                                           ggccTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTAGGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATAcgccgc                                                <  1:239734/149‑1 (MQ=255)                                                                                                 aCCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTAcgcg                                                                                                  >  4:300485/1‑93 (MQ=255)                                                                                                 aCCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTAcgcg                                                                                                  <  3:300485/93‑1 (MQ=255)                                                                                                 aCCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACt                                           <  4:237715/148‑1 (MQ=255)                                                                                                 aCCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTg                                          >  4:197173/1‑149 (MQ=255)                                                                                                   cGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCa                                        >  5:50233/1‑149 (MQ=255)                                                                                                   cGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCa                                        >  2:124303/1‑149 (MQ=255)                                                                                                             aaG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCt                               >  1:249960/1‑149 (MQ=255)                                                                                                              aG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTcc                             <  5:357692/149‑1 (MQ=255)                                                                                                                  cgcAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCCCCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGAc                          <  7:188554/149‑1 (MQ=255)                                                                                                                         ggCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGc                    <  6:50233/148‑1 (MQ=255)                                                                                                                                 cAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGcc           <  1:29545/149‑1 (MQ=255)                                                                                                                                   gtccAGGCCCTTGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCgg         >  3:301257/3‑149 (MQ=255)                                                                                                                                          ccAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTc  <  1:100292/149‑1 (MQ=255)                                                                                                                                           cAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTc  >  7:41989/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                                   CTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAG‑GCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAC‑GGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTC  >  NC_002947/101664‑101948

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