Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,539,330	C→A	7.3%	E270D (GAG→GAT) ‡	PP_2231 ←	DMT superfamily permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,539,330	0	C	A	7.3%	175.6 / 3.3	55	E270D (GAG→GAT) ‡	PP_2231	DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base A (2/2);  total (28/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 = 8.70e-01

CGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGCCCACTGCACTGGCCA  >  NC_002947/2539185‑2539478                                                                                                                                                   |                                                                                                                                                     cGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCg                                                                                                                                                   >  6:64399/1‑149 (MQ=255)         aGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCgccagc                                                                                                                                              >  5:312341/1‑146 (MQ=255)         aGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCgccagc                                                                                                                                              >  8:85793/1‑146 (MQ=255)          gCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGa                                                                                                                                          >  6:158758/1‑149 (MQ=255)            cccGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTTGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACg                                                                                                                                        >  6:144696/1‑149 (MQ=255)              cGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGcc                                                                                                                                      <  6:63465/149‑1 (MQ=255)               gCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGccc                                                                                                                                     >  7:233838/1‑149 (MQ=255)                      gTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCTGGCGAGCgcctgcctggc                                                                                                                                         >  7:248717/1‑134 (MQ=255)                      gTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCTGGCGAGCgcctgcctggc                                                                                                                                         <  8:248717/138‑5 (MQ=255)                          gCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGATGCCCGACGATGCCgg                                                                                                                          >  3:275084/1‑149 (MQ=255)                               ggTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCca                                                                                                                      >  4:275189/1‑148 (MQ=255)                                    tGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCgccagc                                                                                                                                              >  1:134146/1‑119 (MQ=255)                                    tGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCgccagc                                                                                                                                              <  2:134146/119‑1 (MQ=255)                                       acaTCGCGTGGTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGCTCAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGtt                                                                                                              <  1:255065/149‑1 (MQ=255)                                        caTCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCa                                                                                                            <  4:245453/149‑1 (MQ=255)                                          tCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCAcc                                                                                                          >  7:95847/1‑149 (MQ=255)                                            gcgTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCTAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCgg                                                                                                        >  2:333163/1‑149 (MQ=255)                                            gcgTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCg                                                                                                         <  6:312341/148‑1 (MQ=255)                                             cgTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGc                                                                                                       <  1:167783/149‑1 (MQ=255)                                                gtTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCAt                                                                                                     >  2:917/1‑148 (MQ=255)                                                   tCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATaaa                                                                                                  <  3:275189/148‑1 (MQ=255)                                                   tCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATaaa                                                                                                  <  5:158758/148‑1 (MQ=255)                                                    cGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGATGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAAt                                                                                                 <  4:275084/148‑1 (MQ=255)                                                     gACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAAt                                                                                                 >  6:326399/1‑147 (MQ=255)                                                       cAGTTTCGCGAGGGCGCCGCACCTTGAGTTCATTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAgt                                                                                             <  7:137872/149‑1 (MQ=255)                                                         gtatCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAgtgt                                                                                           <  6:42058/146‑1 (MQ=255)                                                          tttCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGtgtg                                                                                          <  8:95847/149‑1 (MQ=255)                                                              gcgAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGgcgc                                                                                      >  5:217108/1‑149 (MQ=255)                                                              gcgAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGgcgc                                                                                      >  2:145596/1‑149 (MQ=255)                                                                gAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCTAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCgg                                                                                    <  1:333163/149‑1 (MQ=255)                                                                      gccgcACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGc                                                                              <  8:233838/149‑1 (MQ=255)                                                                           aCCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAg                                                                                              <  7:305520/128‑1 (MQ=255)                                                                           aCCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAg                                                                                              >  8:305520/1‑128 (MQ=255)                                                                           aCCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGc                                                                              >  5:202066/1‑144 (MQ=255)                                                                           aCCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGc                                                                              <  6:202066/144‑1 (MQ=255)                                                                                gAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCACGGGTAGGCGAGACCGa                                                                      <  2:118933/147‑1 (MQ=255)                                                                                   ttCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCt                                                                 <  5:144696/149‑1 (MQ=255)                                                                                    tCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCt                                                                 <  2:201994/148‑1 (MQ=255)                                                                                           ccaccaCGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAg                                                         >  5:19480/1‑149 (MQ=255)                                                                                              ccacGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAg                                                         >  5:346187/1‑146 (MQ=255)                                                                                              ccacGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAg                                                         <  6:346187/146‑1 (MQ=255)                                                                                                acGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCAt                                                    <  1:917/149‑1 (MQ=255)                                                                                                 cGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCAtt                                                   <  8:38804/149‑1 (MQ=255)                                                                                                  ggCGACCACGATCAGCAGGCCGCCGAGCAGGGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTg                                                  >  1:27645/1‑149 (MQ=255)                                                                                                                gcagGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACt                                    >  8:288049/1‑149 (MQ=255)                                                                                                                    gccgccGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCAc                                >  2:209844/1‑149 (MQ=255)                                                                                                                           agcagcGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTgc                         <  5:64399/149‑1 (MQ=255)                                                                                                                              agcGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTc                              >  1:315557/1‑141 (MQ=255)                                                                                                                              agcGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTc                              <  2:315557/141‑1 (MQ=255)                                                                                                                              agcGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGcc                      >  6:108613/1‑149 (MQ=255)                                                                                                                                     cgc‑caGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGcc                >  7:24900/1‑148 (MQ=255)                                                                                                                                          tggctggcGAGCGCCTGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGCCca              <  8:250495/145‑1 (MQ=255)                                                                                                                                          tggctggcGAGCGCCTGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACAGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGCCca              >  7:250495/2‑146 (MQ=255)                                                                                                                                          aggcaggcGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGCCcactg           >  5:41413/1‑149 (MQ=255)                                                                                                                                               ggcGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGCCcactgcactg      <  1:134555/149‑1 (MQ=255)                                                                                                                                                gcgcTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGCCCACTGCACTgg     >  4:50482/1‑149 (MQ=255)                                                                                                                                                   cTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGCCCACTGCACTGGCCa  <  7:288049/149‑1 (MQ=255)                                                                                                                                                   cTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGCCCACTGCACTGGCCa  <  1:145596/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGACCACCACGGCGACCACGATCAGCAGGCCGCCGAGCAGCGCCACGC‑CAGGCAGGCGCTCGCCAGCCAGACGCCCGACGATGCCGGCCCACACCGGTTCACCGGCATAAATGAGTGTGGCGCGGGTAGGCGAGACCGACTTCTGTGCCCAGTTCATTGCCACCTGGATCACTGCACTCATTGCGCCCAGGCCCACTGCACTGGCCA  >  NC_002947/2539185‑2539478

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