Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,977,720	G→A	E66K (GAA→AAA)	rffT →	4‑acetamido‑4,6‑dideoxy‑D‑galactose transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,977,720	0	G	A	100.0%	205.9 / NA	59	E66K (GAA→AAA)	rffT	4‑acetamido‑4,6‑dideoxy‑D‑galactose transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/31);  total (28/31)

GACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGT  >  NC_000913/3977597‑3977834                                                                                                                            |                                                                                                                   gACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAg                                                                                                  <  1:801679/140‑1 (MQ=255)  aCGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGc                                                                                                 >  2:1251249/1‑140 (MQ=255)   cgcgCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCa                                                                                                <  1:780955/140‑1 (MQ=255)     cgcTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTGAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCaaa                                                                                              <  1:1282312/140‑1 (MQ=255)          gCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTa                                                                                         <  2:408689/140‑1 (MQ=255)              cgacgaGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCg                                                                                     >  2:691040/1‑140 (MQ=255)                acgaGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTc                                                                                   >  1:1023071/1‑140 (MQ=255)                  gagcgagcACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTcag                                                                                 >  2:589098/1‑140 (MQ=255)                   agcgagcACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTcagc                                                                                <  1:1290175/140‑1 (MQ=255)                         cACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGtt                                                                          <  1:245047/140‑1 (MQ=255)                         cACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGtt                                                                          >  2:536536/1‑140 (MQ=255)                          aCGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGttt                                                                         <  2:1023071/140‑1 (MQ=255)                             cgcgcgAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGtttttt                                                                      >  1:259668/1‑140 (MQ=255)                              gcgcgAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTttc                                                                     >  2:1088664/1‑140 (MQ=255)                                  gAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTcc                                                                 <  2:520724/140‑1 (MQ=255)                                       tATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGt                                                            <  2:197683/140‑1 (MQ=255)                                          ggTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAg                                                         <  1:1085298/140‑1 (MQ=255)                                            ttgttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGtt                                                       <  2:58346/140‑1 (MQ=255)                                            ttgttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGtt                                                       <  1:589098/140‑1 (MQ=255)                                            ttgttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGtt                                                       <  1:623145/140‑1 (MQ=255)                                              gttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCa                                                     >  2:547314/1‑140 (MQ=255)                                              gttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCa                                                     >  1:255693/1‑140 (MQ=255)                                              gttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCa                                                     <  1:1088664/140‑1 (MQ=255)                                              gttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCa                                                     >  2:594210/1‑140 (MQ=255)                                              gttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCa                                                     >  2:1057482/1‑140 (MQ=255)                                               ttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCaa                                                    >  2:1119145/1‑140 (MQ=255)                                               ttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCaa                                                    >  2:849288/1‑140 (MQ=255)                                               ttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCaa                                                    >  1:247131/1‑140 (MQ=255)                                                 ggCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATc                                                  <  1:849288/140‑1 (MQ=255)                                                  gCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGTTCAGTTCAATcc                                                 >  1:346268/1‑140 (MQ=255)                                                   cAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATccc                                                <  2:1225739/140‑1 (MQ=255)                                                      ggACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCcaca                                             >  2:233927/1‑140 (MQ=255)                                                             ggCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGtggc                                      <  2:1133883/140‑1 (MQ=255)                                                                ttAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGtggctgg                                   <  1:1119145/140‑1 (MQ=255)                                                                ttAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGtggctgg                                   <  2:1056753/140‑1 (MQ=255)                                                                ttAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGtggctgg                                   <  2:642150/140‑1 (MQ=255)                                                                  aaGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTggctggct                                 <  1:1251249/140‑1 (MQ=255)                                                                  aaGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTggctggct                                 <  2:410703/140‑1 (MQ=255)                                                                    gTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGctct                               >  1:831227/1‑140 (MQ=255)                                                                       aTAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctgct                            <  2:1008806/140‑1 (MQ=255)                                                                       aTAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctgct                            <  2:346268/140‑1 (MQ=255)                                                                       aTAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctgct                            >  1:637449/1‑140 (MQ=255)                                                                        tAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctgctg                           <  1:898384/140‑1 (MQ=255)                                                                         aGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACtgtg                                        >  1:1100936/1‑126 (MQ=255)                                                                         aGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACtgtg                                        <  2:1100936/126‑1 (MQ=255)                                                                         aGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGa                          >  1:566310/1‑140 (MQ=255)                                                                         aGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGa                          >  1:122383/1‑140 (MQ=255)                                                                             gTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACtgtg                                        >  1:1047081/1‑122 (MQ=255)                                                                             gTCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACtgtg                                        <  2:1047081/122‑1 (MQ=255)                                                                              tCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAgtggt                     >  2:495185/1‑140 (MQ=255)                                                                              tCCGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAgtggt                     >  2:588473/1‑140 (MQ=255)                                                                               ccGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAgtggtg                    <  1:1057482/140‑1 (MQ=255)                                                                               ccGGCGCTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAgtggtg                    >  1:1031406/1‑140 (MQ=255)                                                                                   cgcTTTCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCAt                <  2:566310/140‑1 (MQ=255)                                                                                        tCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGc           >  1:849505/1‑140 (MQ=255)                                                                                        tCTGTGCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGc           >  1:957503/1‑140 (MQ=255)                                                                                           gtgCAATTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCcca        <  2:259668/140‑1 (MQ=255)                                                                                                aTTTTTCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCccagccag   <  2:1141110/140‑1 (MQ=255)                                                                                                 tttttCCCTGGG‑AAAAATTGC‑ACAGAAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGt  >  2:295360/1‑140 (MQ=255)                                                                                                                            |                                                                                                                   GACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGT  >  NC_000913/3977597‑3977834

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