Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,409,489	T→C	G140G (GGT→GGC)	prmA →	methyltransferase for 50S ribosomal subunit protein L11

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,409,489	0	T	C	100.0%	207.9 / NA	61	G140G (GGT→GGC)	prmA	methyltransferase for 50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/29);  total (32/29)

GGAGCGCGAATGGATGGATAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATT  >  NC_000913/3409354‑3409624                                                                                                                                        |                                                                                                                                        ggAGCGCGAATGGATGGATAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCAccc                                                                                                                                     <  2:30390/140‑1 (MQ=255) ggAGCGCGAATGGATGGATAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCAccc                                                                                                                                     >  1:715116/1‑140 (MQ=255)   aGCGCGAATGGATGGATAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCAt                                                                                                                                   <  2:131661/140‑1 (MQ=255)       cgAATGGATGGATAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCaa                                                                                                                               >  2:291022/1‑140 (MQ=255)        gAATGGATGGATAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAac                                                                                                                              <  2:695566/140‑1 (MQ=255)        gAATGGATGGATAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAac                                                                                                                              <  2:710510/140‑1 (MQ=255)           tggatggaTAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAaccac                                                                                                                           >  1:705971/1‑140 (MQ=255)                 gatAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACctctct                                                                                                                     <  1:644973/140‑1 (MQ=255)                  atAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCtg                                                                                                                    >  1:432627/1‑140 (MQ=255)                       ttCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAaccac                                                                                                                           <  2:905116/128‑1 (MQ=255)                       ttCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAaccac                                                                                                                           >  1:905116/1‑128 (MQ=255)                            cccGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAaccac                                                                                                                           >  2:955846/1‑123 (MQ=255)                            cccGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAaccac                                                                                                                           <  1:955846/123‑1 (MQ=255)                             ccGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCGAGGGCTGGCGTTTGGTACGGGCACCCATCCAaccac                                                                                                                           >  2:180002/1‑122 (MQ=255)                             ccGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAaccac                                                                                                                           <  1:180002/122‑1 (MQ=255)                              cGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATgg                                                                                                        <  1:291022/140‑1 (MQ=255)                                              aCGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGAtt                                                                                        >  2:271881/1‑140 (MQ=255)                                                  cTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGAGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAAc                                                                                    >  1:167358/1‑140 (MQ=255)                                                     tgGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCgg                                                                                 >  2:659686/1‑140 (MQ=255)                                                          cTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTaaaa                                                                            >  2:345599/1‑140 (MQ=255)                                                            gCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACa                                                                          >  1:837955/1‑140 (MQ=255)                                                            gCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACa                                                                          >  1:877986/1‑140 (MQ=255)                                                             cccTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAg                                                                         >  1:178768/1‑140 (MQ=255)                                                                       cGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCCCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACtt                                                               >  1:328278/1‑140 (MQ=255)                                                                           aTGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGc                                                           >  2:120884/1‑140 (MQ=255)                                                                            tgtgCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTaaaa                                                                            <  1:974758/122‑1 (MQ=255)                                                                            tgtgCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTaaaa                                                                            >  2:974758/1‑122 (MQ=255)                                                                            tgtgCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCt                                                          >  1:165043/1‑140 (MQ=255)                                                                                ccGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCa                                                                      <  2:525412/124‑1 (MQ=255)                                                                                ccGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCa                                                                      >  1:525412/1‑124 (MQ=255)                                                                                   gACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTc                                                   >  2:554413/1‑140 (MQ=255)                                                                                   gACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTc                                                   <  2:715116/140‑1 (MQ=255)                                                                                     cGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCg                                                 <  1:120884/140‑1 (MQ=255)                                                                                        aaaCGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCa                                              >  1:623128/1‑140 (MQ=255)                                                                                        aaaCGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCa                                              >  2:321916/1‑140 (MQ=255)                                                                                         aaCGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTc                                                   <  2:669351/134‑1 (MQ=255)                                                                                         aaCGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTc                                                   >  1:669351/1‑134 (MQ=255)                                                                                                     gTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCCCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATcgc                                 <  2:328278/140‑1 (MQ=255)                                                                                                         tGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCgcggcg                             <  2:837955/140‑1 (MQ=255)                                                                                                         tGTTAGATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCgcggcg                             <  1:377177/140‑1 (MQ=255)                                                                                                              gATCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGaa                        >  2:101165/1‑140 (MQ=255)                                                                                                               aTCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGaaa                       <  1:101165/140‑1 (MQ=255)                                                                                                               aTCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGaaa                       >  1:439187/1‑140 (MQ=255)                                                                                                               aTCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGaaa                       <  2:429636/140‑1 (MQ=255)                                                                                                                tCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGaaa                       >  1:145021/1‑139 (MQ=255)                                                                                                                tCCAGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGaaa                       <  2:145021/139‑1 (MQ=255)                                                                                                                   aGGGCTGGCGTTTGGTACGGGCACCCATCCAACTACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTgg                   >  2:668728/1‑140 (MQ=255)                                                                                                                   aGGGCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTgg                   >  2:552997/1‑140 (MQ=255)                                                                                                                     ggCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATcgc                                 >  1:744420/1‑124 (MQ=255)                                                                                                                     ggCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATcgc                                 <  2:744420/124‑1 (MQ=255)                                                                                                                     ggCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGt                 >  2:928586/1‑140 (MQ=255)                                                                                                                     ggCTGGCGTTGGGAACGGCAACCATTCCAACCACTTCTTTGGGCCGGCAAGGCTTCGCCACCTCCGTTTTAACCGGAAAAACATCAATCGACTTTGGCTGTGGTTCGGGCATTCGGGGGTTCGCGGGCCGGAAACTgggg                 >  2:386708/1‑139 (MQ=255)                                                                                                                      gCTGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTg                <  1:659686/140‑1 (MQ=255)                                                                                                                        tGGCGTTTGGTACGGGCACCCATCCAACCACCTCTCTCTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTgca              >  1:371024/1‑140 (MQ=255)                                                                                                                          gCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTgcagc            >  2:114133/1‑140 (MQ=255)                                                                                                                          gCGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTgcagc            >  1:808360/1‑140 (MQ=255)                                                                                                                           cGTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTgcagca           <  1:928586/140‑1 (MQ=255)                                                                                                                            gTTTGGTACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCaa          <  2:432627/140‑1 (MQ=255)                                                                                                                                  tACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCa    <  1:321916/140‑1 (MQ=255)                                                                                                                                  tACGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCa    <  1:163910/140‑1 (MQ=255)                                                                                                                                   aCGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCAt   <  2:178768/140‑1 (MQ=255)                                                                                                                                    cGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCAtt  <  2:557348/140‑1 (MQ=255)                                                                                                                                    cGGGCACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCAtt  <  1:554413/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        GGAGCGCGAATGGATGGATAATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCCAGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATT  >  NC_000913/3409354‑3409624

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: ‑