Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,101,056	A→C	S83R (AGC→CGC)	sodA →	superoxide dismutase, Mn

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,101,056	0	A	C	100.0%	200.8 / NA	57	S83R (AGC→CGC)	sodA	superoxide dismutase, Mn
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (25/32);  total (25/32)

CCTACGTAAACAACGCCAACGCGGCGCTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACAGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATG  >  NC_000913/4100910‑4101201                                                                                                                                                   |                                                                                                                                                  ccTACGTAAACAACGCCAACGCGGCGCTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCt                                                                                                                                               <  2:291562/151‑1 (MQ=255)   tACGTAAACAACGCCAACGCGGCGCTGGAAAGCCTGCCAGAATTTGTCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGt                                                                                                                                             >  1:470075/1‑151 (MQ=255)    aCGTAAACAACGCCAACGCGGCGCTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGt                                                                                                                                             >  1:484116/1‑150 (MQ=255)      gTAAACAACGCCAACGCGGCGCTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCt                                                                                                                                          <  2:211074/151‑1 (MQ=255)                ccaacgcGGCGCTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTTACGCTAACCCCCGCCCGTTCTGGAATGGTc                                                                                                                                 >  1:393820/1‑150 (MQ=255)                ccaacgcGGCGCTGGAAAGCCGGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCt                                                                                                                                <  1:299845/151‑1 (MQ=255)                  aacgcGGCGCTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTg                                                                                                                               >  1:532921/1‑150 (MQ=255)                  aacgcGGCGCTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGa                                                                                                                              <  2:118553/151‑1 (MQ=255)                     gcggcgCTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGaaa                                                                                                                            >  1:85370/1‑150 (MQ=255)                          gcTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAgg                                                                                                                       <  2:719911/150‑1 (MQ=255)                          gcTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGc                                                                                                                      <  1:130456/151‑1 (MQ=255)                             ggAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGcac                                                                                                                    <  2:81441/150‑1 (MQ=255)                             ggAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGca                                                                                                                     <  2:423523/149‑1 (MQ=255)                                   ccTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTg                                                                                                             >  2:35437/1‑151 (MQ=255)                                          gAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTg                                                                                                      >  2:589452/1‑151 (MQ=255)                                           aaTTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTg                                                                                                      >  2:444775/1‑150 (MQ=255)                                            aTTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGAc                                                                                                    >  2:159051/1‑151 (MQ=255)                                                 ccAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGaa                                                                                               <  2:22488/151‑1 (MQ=255)                                                  cAACCTGCCGGGTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCCCGCTAACCACCGCCTGTTATGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCGGa                                                                                                >  1:334492/1‑149 (MQ=255)                                                      cTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCgg                                                                                          >  1:651663/1‑151 (MQ=255)                                                           ggTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATc                                                                                     >  1:266817/1‑151 (MQ=255)                                                            gTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCg                                                                                    >  2:339830/1‑151 (MQ=255)                                                                 agaaCTGATCACCACGCTGGACCAGCTGCCAGCAGACAAGAAACCCGTGCTGCGCAGCAACGCTGGCGTCCACGCTAAACACCGGCTGTTCTGGAAAGGTCTGAAAAAAGGCCCCACCCTGCAGGGTGACCTGAAAGCGGCTATCGaa                                                                                  <  1:809265/148‑1 (MQ=255)                                                                          cgccAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTGCTGCGCAACAACGCTGGCGGTGACGCTACCCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCgg                                                                       <  1:804727/148‑1 (MQ=255)                                                                             cAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGATCTGGAAAGGTCTGTAAAAAGGCACCACCCTGCAGGGGGACCTTAAAGCGGCTATCGAACGTGACTTCGGCttc                                                                   >  1:873227/1‑149 (MQ=255)                                                                              aaaCTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCg                                                                  <  2:846288/151‑1 (MQ=255)                                                                               aaCTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGt                                                                 <  2:178090/151‑1 (MQ=255)                                                                               aaCTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGt                                                                 <  1:721361/151‑1 (MQ=255)                                                                                       cAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACt                                                         <  2:122082/151‑1 (MQ=255)                                                                                             ccAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCaa                                                     >  2:590823/1‑149 (MQ=255)                                                                                              cagcagACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGc                                                  <  2:561472/151‑1 (MQ=255)                                                                                               agcagACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCa                                                 <  2:659128/151‑1 (MQ=255)                                                                                               agcagACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCa                                                 <  2:441412/151‑1 (MQ=255)                                                                                                gcagACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGGGACTTCGGCTCCGTTGATAACTTCAAAGc                                                  >  1:664939/1‑149 (MQ=255)                                                                                                  agACAAGAAACCGGGACTGCGCAAAAACGCAGGGGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACGCTGCATGGTGGCCTGAAAGCGGCTATCGAACGTGGCTTCGGCTCCGTTGATAACTTCAAAGCAg                                                <  1:419674/149‑1 (MQ=255)                                                                                                  agACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGa                                               >  2:81462/1‑150 (MQ=255)                                                                                                    aCAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAAtt                                            <  2:578066/151‑1 (MQ=255)                                                                                                    aCAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAAtt                                            >  2:24402/1‑151 (MQ=255)                                                                                                          aaaCCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGaaaa                                      <  2:574283/151‑1 (MQ=255)                                                                                                          aaaCCGTACTGCGCAACAACGCTGGCGGCCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGaaaa                                      <  2:616873/151‑1 (MQ=255)                                                                                                           aaCCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGaaaaa                                     <  1:212621/151‑1 (MQ=255)                                                                                                           aaCCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGaaaaa                                     <  2:311107/151‑1 (MQ=255)                                                                                                           aaCCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGaaaa                                      >  1:422157/1‑150 (MQ=255)                                                                                                                 aCTGCGCAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCa                               <  2:403286/151‑1 (MQ=255)                                                                                                                    gcgcAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCt                            <  2:845029/151‑1 (MQ=255)                                                                                                                     cgcAACAACGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCtt                           <  2:154621/151‑1 (MQ=255)                                                                                                                           aacGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCt                     >  2:126864/1‑151 (MQ=255)                                                                                                                            acGCTGGCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCtt                    >  1:889091/1‑151 (MQ=255)                                                                                                                              gCTGGCGGTCACGCTAACCACCGCCTGGTCTGGGAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGGACGCGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCttt                   >  1:214611/1‑150 (MQ=255)                                                                                                                                  gCGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTc              <  1:851929/151‑1 (MQ=255)                                                                                                                                   cGGTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTcc             <  1:206973/151‑1 (MQ=255)                                                                                                                                    ggTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCg            >  2:307691/1‑151 (MQ=255)                                                                                                                                    ggTCACGCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCg            >  1:721527/1‑151 (MQ=255)                                                                                                                                         ctcTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTgg       >  2:743599/3‑151 (MQ=255)                                                                                                                                         cGCTAACCACCGCCTGTTCTGGAAATGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCt         <  2:734023/149‑1 (MQ=255)                                                                                                                                          gCTAACCACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTggg      <  1:15884/151‑1 (MQ=255)                                                                                                                                               ccACCGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATg  <  1:104506/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CCTACGTAAACAACGCCAACGCGGCGCTGGAAAGCCTGCCAGAATTTGCCAACCTGCCGGTTGAAGAGCTGATCACCAAACTGGACCAGCTGCCAGCAGACAAGAAAACCGTACTGCGCAACAACGCTGGCGGTCACGCTAACCACAGCCTGTTCTGGAAAGGTCTGAAAAAAGGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATG  >  NC_000913/4100910‑4101201

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