Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,783,691	C→T	7.4%	R1208Q (CGG→CAG)	pvdJ ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,783,691	0	C	T	7.4%	157.4 / 2.4	55	R1208Q (CGG→CAG)	pvdJ	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (29/26)
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.27e-01

GCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGG  >  NC_002947/4783545‑4783836                                                                                                                                                     |                                                                                                                                                  gcgcCCATTGGCTCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGt                                                                                                                                                 <  7:97120/149‑1 (MQ=255) gcgcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGt                                                                                                                                                 >  2:30336/1‑149 (MQ=255) gcgcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGt                                                                                                                                                 >  1:34715/1‑149 (MQ=255)   gcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTg                                                                                                                                               >  8:3420/1‑149 (MQ=255)     ccATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCt                                                                                                                                             >  6:132436/1‑149 (MQ=255)     ccATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCt                                                                                                                                             <  2:190433/149‑1 (MQ=255)        ttGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGAc                                                                                                                                          <  1:104878/149‑1 (MQ=255)        ttGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGAc                                                                                                                                          <  2:34715/149‑1 (MQ=255)           gcgcAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCACTAGGCTAGCCTGATGCTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCAt                                                                                                                                       >  2:38024/1‑149 (MQ=255)             gcAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCAACg                                                                                                                                     >  7:201659/1‑149 (MQ=255)                      ccGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCg                                                                                                                            <  2:7940/149‑1 (MQ=255)                       cGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCg                                                                                                                            >  3:106189/1‑148 (MQ=255)                           cTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACaa                                                                                                                       >  1:234696/1‑149 (MQ=255)                                    aGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCt                                                                                                              <  2:234696/149‑1 (MQ=255)                                       cTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGaa                                                                                                           <  2:40512/149‑1 (MQ=255)                                          aCCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAAc                                                                                                          <  2:21930/147‑1 (MQ=255)                                                  cccTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCAc                                                                                                >  5:265000/1‑149 (MQ=255)                                                   ccTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCTGGCCCCCAGTTCATGCACCGCCAGCAA‑CACCGGTTGATCACAAACGGCCATTTCGAACAATGACCCCCTGAACACCTCGcacc                                                                                                >  1:80990/1‑146 (MQ=255)                                                        gCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGc                                                                                            >  5:157339/1‑147 (MQ=255)                                                              ctgctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGa                                                                                    <  7:3420/149‑1 (MQ=255)                                                               tgctgGTAGGCCGTCTGCATGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCATCAACCAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGaa                                                                                                           <  1:191824/126‑1 (MQ=255)                                                               tgctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGaa                                                                                                           >  2:191824/1‑125 (MQ=255)                                                                 ctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAgg                                                                                  <  6:157339/148‑1 (MQ=255)                                                                  tgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCt                                                                                >  1:72870/1‑149 (MQ=255)                                                                    gTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCt                                                                                <  6:265000/147‑1 (MQ=255)                                                                    gTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTg                                                                               <  2:223524/148‑1 (MQ=255)                                                                               gCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTc                                                                    <  4:106189/148‑1 (MQ=255)                                                                                   gTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTgg                                                               <  5:132436/149‑1 (MQ=255)                                                                                   gTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGAAccaagtttc                                                               >  8:94141/1‑142 (MQ=255)                                                                                    tCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTggc                                                              <  2:72870/149‑1 (MQ=255)                                                                                           aCAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAg                                                                                         >  3:210553/1‑115 (MQ=255)                                                                                           aCAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAg                                                                                         <  4:210553/115‑1 (MQ=255)                                                                                                aTGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAg                                                  <  1:30336/149‑1 (MQ=255)                                                                                                    gcCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCt                                              >  1:182962/1‑149 (MQ=255)                                                                                                    gcCAAGACACCCCATC‑CACCACCAGGAGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCt                                              >  7:155879/1‑149 (MQ=255)                                                                                                     ccAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTg                                             >  2:64434/1‑149 (MQ=255)                                                                                                          acacCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCAACGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAgg                                        <  8:201659/149‑1 (MQ=255)                                                                                                            acCCCATC‑CACGACCAGGTGATGCACCGCCAGCAA‑CAGCCGGTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCGACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGGGGTTTCTAGCGCTGCAAGGCt                                      <  8:155879/149‑1 (MQ=255)                                                                                                            acCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCt                                      >  6:52809/1‑149 (MQ=255)                                                                                                              cccATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGt                                    >  7:153937/1‑149 (MQ=255)                                                                                                                 aTC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCa                                 <  1:64434/149‑1 (MQ=255)                                                                                                                     caccaccaGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCagag                              >  3:163096/1‑149 (MQ=255)                                                                                                                          ccaGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGAc                          >  1:36465/1‑148 (MQ=255)                                                                                                                          ccaGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACt                         <  2:36465/149‑1 (MQ=255)                                                                                                                            aGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTAc                       >  2:159890/1‑149 (MQ=255)                                                                                                                             ggTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAgc                                                 <  3:128257/122‑1 (MQ=255)                                                                                                                             ggTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAgc                                                 >  4:128257/1‑122 (MQ=255)                                                                                                                                gATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTg                    >  4:244273/1‑148 (MQ=255)                                                                                                                                         tcAGCAA‑CGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGgcggcg                                                          >  6:243603/2‑101 (MQ=255)                                                                                                                                         tcAGCAA‑CGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGgcggcg                                                          <  5:243603/100‑1 (MQ=255)                                                                                                                                         tcAGCAA‑CGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggagg             <  3:143722/145‑1 (MQ=255)                                                                                                                                         tcAGCAA‑CGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggagg             >  4:143722/2‑146 (MQ=255)                                                                                                                                           aGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTg        <  7:94141/149‑1 (MQ=255)                                                                                                                                            gCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTg        >  8:18872/1‑148 (MQ=255)                                                                                                                                                  aGCCGTTGCTGCCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGGTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCgg  >  6:141024/1‑149 (MQ=255)                                                                                                                                                    ccGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTAcc                      >  7:45832/1‑127 (MQ=255)                                                                                                                                                    ccGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTAcc                      <  8:45832/127‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGG  >  NC_002947/4783545‑4783836

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