Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,811,367	T→C	intergenic (+67/‑60)	mprA → / → emrA	DNA‑binding transcriptional repressor MprA/multidrug efflux pump membrane fusion protein EmrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,811,367	0	T	C	100.0%	218.4 / NA	59	intergenic (+67/‑60)	mprA/emrA	DNA‑binding transcriptional repressor MprA/multidrug efflux pump membrane fusion protein EmrA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/35);  total (24/35)

GCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCTGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTCCTCCT  >  NC_000913/2811243‑2811503                                                                                                                             |                                                                                                                                         gCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAg                                                                                                                           <  1:1015709/140‑1 (MQ=255) gCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAg                                                                                                                           <  2:857594/140‑1 (MQ=255)       ccGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTc                                                                                                                     >  1:718148/1‑140 (MQ=255)              gACCAGATGGAACAAGACGGTGTGGTTCTGGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAg                                                                                                              <  2:800751/140‑1 (MQ=255)               aCCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGc                                                                                                             >  1:620133/1‑140 (MQ=255)               aCCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGc                                                                                                             >  1:894043/1‑140 (MQ=255)               aCCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGc                                                                                                             >  1:829722/1‑140 (MQ=255)                   gATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGAt                                                                                                         <  1:652907/140‑1 (MQ=255)                    aTGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATg                                                                                                        >  1:323412/1‑140 (MQ=255)                    aTGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATg                                                                                                        <  1:504984/140‑1 (MQ=255)                          cAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTaa                                                                                                  <  1:371329/140‑1 (MQ=255)                          cAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTaa                                                                                                  <  2:551330/140‑1 (MQ=255)                              aCGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTaagaag                                                                                              >  2:44654/1‑140 (MQ=255)                              aCGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTaagaag                                                                                              >  1:238419/1‑140 (MQ=255)                               cGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAagaaga                                                                                             <  2:56810/140‑1 (MQ=255)                                  tgtgGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCg                                                                                          >  2:152310/1‑140 (MQ=255)                                  tgtgGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCg                                                                                          >  2:860241/1‑140 (MQ=255)                                  tgtgGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCg                                                                                          >  2:264329/1‑140 (MQ=255)                                  tgtgGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCg                                                                                          >  2:404087/1‑140 (MQ=255)                                  tgtgGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCg                                                                                          >  1:177966/1‑140 (MQ=255)                                  tgtgGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCg                                                                                          >  1:114276/1‑140 (MQ=255)                                     ggTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTgg                                                                                       <  2:477491/140‑1 (MQ=255)                                       ttCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGgag                                                                                     <  2:1024531/140‑1 (MQ=255)                                         ctcGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTATGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGaa                                                                                   <  2:934712/140‑1 (MQ=255)                                          tcGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAAc                                                                                  <  2:448335/140‑1 (MQ=255)                                            gAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCTGATGAGTTAAGAAGATCGTGGAGAACaa                                                                                >  1:973246/1‑140 (MQ=255)                                                  aTGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAg                                                                          <  1:744090/140‑1 (MQ=255)                                                    gAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAgcg                                                                        <  1:626948/140‑1 (MQ=255)                                                              gTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGga                                                              <  2:323412/140‑1 (MQ=255)                                                               tCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGgag                                                             >  2:247023/1‑140 (MQ=255)                                                                    tcGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCa                                                        <  2:797166/140‑1 (MQ=255)                                                                        tCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAAcc                                                    >  2:234976/1‑140 (MQ=255)                                                                             aaTCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCgca                                               <  2:305577/140‑1 (MQ=255)                                                                                ccAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCTGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCgcagca                                            <  2:973246/140‑1 (MQ=255)                                                                                ccAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCgcagca                                            <  2:718148/140‑1 (MQ=255)                                                                                ccAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCgcagca                                            <  1:404087/140‑1 (MQ=255)                                                                                ccAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCgcagca                                            <  2:944669/140‑1 (MQ=255)                                                                                  aGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAAc                                          <  1:10048/140‑1 (MQ=255)                                                                                  aGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAAc                                          <  1:215358/140‑1 (MQ=255)                                                                                  aGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAAc                                          <  2:416974/140‑1 (MQ=255)                                                                                         taAAAGAAAAAATGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTaaa                                   >  1:963001/1‑140 (MQ=255)                                                                                                aaaaaTGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGc                            >  1:737104/1‑140 (MQ=255)                                                                                                 aaaaTGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCg                           <  1:177257/140‑1 (MQ=255)                                                                                                   aaTGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGc                         >  1:1030180/1‑140 (MQ=255)                                                                                                    aTGACTGGCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCa                        <  2:801127/140‑1 (MQ=255)                                                                                                           gCCAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTaa                 >  2:915998/1‑140 (MQ=255)                                                                                                             cAGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGc               >  2:298967/1‑140 (MQ=255)                                                                                                              aGCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCg              <  2:829722/140‑1 (MQ=255)                                                                                                               gCATCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGt             <  2:540030/140‑1 (MQ=255)                                                                                                                 aTCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCt           <  2:703919/140‑1 (MQ=255)                                                                                                                  tCGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTg          <  1:44654/140‑1 (MQ=255)                                                                                                                   cGCAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGc         <  2:115058/140‑1 (MQ=255)                                                                                                                     cAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTc            <  2:1010494/135‑1 (MQ=255)                                                                                                                     cAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTc            >  1:1010494/1‑135 (MQ=255)                                                                                                                     cAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGctc       <  1:845200/140‑1 (MQ=255)                                                                                                                     cAACATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGctc       >  1:810596/1‑140 (MQ=255)                                                                                                                      aaCATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGctcc      <  2:177966/140‑1 (MQ=255)                                                                                                                      aaCATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGctcc      <  2:114276/140‑1 (MQ=255)                                                                                                                       aCATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGctcct     <  2:238409/140‑1 (MQ=255)                                                                                                                       aCATGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGctcct     <  2:126878/140‑1 (MQ=255)                                                                                                                          tGCCGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTcctcct  >  1:6757/1‑140 (MQ=255)                                                                                                                             |                                                                                                                                         GCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCTGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTCCTCCT  >  NC_000913/2811243‑2811503

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