Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,380,499	T→C	T49A (ACA→GCA)	menF ←	isochorismate synthase MenF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,380,499	0	T	C	100.0%	230.3 / NA	61	T49A (ACA→GCA)	menF	isochorismate synthase MenF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/27);  total (34/27)

CGTAAGTCGGCGTGTTCCGGGTGCTGGCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGTTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGATTTTCCAGCGCCGTAGTA  >  NC_000913/2380366‑2380632                                                                                                                                      |                                                                                                                                      cGTAAGTCGGCGTGTTCCGGGTGCTGGCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTg                                                                                                                                 >  1:868076/1‑140 (MQ=255) cGTAAGTCGGCGTGTTCCGGGTGCTGGCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTg                                                                                                                                 >  1:52763/1‑140 (MQ=255) cGTAAGTCGGCGTGTTCCGGGTGCTGGCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTg                                                                                                                                 <  1:459570/140‑1 (MQ=255)           cGTGTTCCGGGTGCTGGCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGgccagc                                                                                                                       <  2:455714/140‑1 (MQ=255)                    ggTGCTGGCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAg                                                                                                              >  2:788113/1‑140 (MQ=255)                    ggTGCTGGCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAg                                                                                                              <  2:585520/140‑1 (MQ=255)                    ggTGCTGGCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAg                                                                                                              <  2:24458/140‑1 (MQ=255)                          ggCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCa                                                                                                        >  2:83226/1‑140 (MQ=255)                          ggCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCa                                                                                                        >  1:317948/1‑140 (MQ=255)                            cGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCaaa                                                                                                      >  2:231779/1‑140 (MQ=255)                             gaagaaAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCaaaa                                                                                                     <  1:512175/140‑1 (MQ=255)                               agaaAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGc                                                                                                   <  1:243459/140‑1 (MQ=255)                                gaaAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCg                                                                                                  >  2:695925/1‑140 (MQ=255)                                gaaAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCg                                                                                                  >  2:649214/1‑140 (MQ=255)                                   aaCGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCt                                                                                               >  2:619908/1‑140 (MQ=255)                                     cGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCttt                                                                                             >  2:684519/1‑140 (MQ=255)                                     cGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCttt                                                                                             <  2:608749/140‑1 (MQ=255)                                       ttGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTga                                                                                           <  2:145147/140‑1 (MQ=255)                                         gtgCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTgaga                                                                                         <  2:52763/140‑1 (MQ=255)                                           gCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAgg                                                                                       <  2:22389/140‑1 (MQ=255)                                           gCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAgg                                                                                       <  1:562592/140‑1 (MQ=255)                                              tGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGTCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGaa                                                                                    >  2:640064/1‑140 (MQ=255)                                              tGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGaa                                                                                    <  1:649214/140‑1 (MQ=255)                                              tGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGaa                                                                                    <  1:695925/140‑1 (MQ=255)                                                      cgacgTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAatat                                                                            >  2:348434/1‑140 (MQ=255)                                                          gTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGAt                                                                        >  1:162187/1‑140 (MQ=255)                                                           tAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATa                                                                       >  2:827271/1‑140 (MQ=255)                                                              aaCGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAAcc                                                                    >  2:537379/1‑140 (MQ=255)                                                                  ggTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGaa                                                                >  2:724729/1‑140 (MQ=255)                                                                        cgcgCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCgg                                                          >  2:340767/1‑140 (MQ=255)                                                                          cgcCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGgt                                                        >  2:494360/1‑140 (MQ=255)                                                                            cccAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGgtgt                                                      <  1:401017/140‑1 (MQ=255)                                                                             ccAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGgtg                                                       >  2:255924/1‑138 (MQ=255)                                                                             ccAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGgtg                                                       <  1:255924/138‑1 (MQ=255)                                                                               aGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGc                                                   >  1:369347/1‑140 (MQ=255)                                                                                      cAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAAtc                                            <  1:494360/140‑1 (MQ=255)                                                                                           tcttcATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCg                                                    <  2:766016/127‑1 (MQ=255)                                                                                           tcttcATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCg                                                    >  1:766016/1‑127 (MQ=255)                                                                                           tcttcATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTg                                       >  2:657647/1‑140 (MQ=255)                                                                                           tcttcATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTg                                       >  1:246588/1‑140 (MQ=255)                                                                                             ttcATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCg                                     <  1:684519/140‑1 (MQ=255)                                                                                                    ccATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATg                              <  1:537379/140‑1 (MQ=255)                                                                                                    ccATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATg                              <  1:231779/140‑1 (MQ=255)                                                                                                     cATTACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATgg                             >  2:576876/1‑140 (MQ=255)                                                                                                        tACGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGgcgc                          <  2:317948/140‑1 (MQ=255)                                                                                                         aCGTTGTTGCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCa                         >  2:768893/1‑140 (MQ=255)                                                                                                               ttgCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCa                         >  2:331721/1‑134 (MQ=255)                                                                                                               ttgCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCa                         <  1:331721/134‑1 (MQ=255)                                                                                                               ttgCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGAt                   >  2:123350/1‑140 (MQ=255)                                                                                                               ttgCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGAt                   >  2:175500/1‑140 (MQ=255)                                                                                                                 gCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGAttt                 >  2:740867/1‑140 (MQ=255)                                                                                                                 gCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGAttt                 >  1:538533/1‑140 (MQ=255)                                                                                                                 gCCAGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGAttt                 >  1:748022/1‑140 (MQ=255)                                                                                                                    aGTAGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGATTTTcc              >  2:443540/1‑140 (MQ=255)                                                                                                                       aGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGATTTTCCAgc           <  1:281668/140‑1 (MQ=255)                                                                                                                       aGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGATTTTCCAgc           <  1:340767/140‑1 (MQ=255)                                                                                                                       aGAATTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGATTTTCCAgc           <  2:246588/140‑1 (MQ=255)                                                                                                                         aaTTGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGATTTTCCAgcgc         <  1:123350/140‑1 (MQ=255)                                                                                                                           ttGCGGGTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTATATTTTCCAGCGCCg       >  2:732960/1‑140 (MQ=255)                                                                                                                               gggTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGATTTTCCAGCGCCgtagt   <  2:748022/140‑1 (MQ=255)                                                                                                                                ggTATGCTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGATTTTCCAGCGCCgtagta  >  2:201520/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                      CGTAAGTCGGCGTGTTCCGGGTGCTGGCGAAGAAAACGTTGTGCCTGGTCCAACGACGTAAAACGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTATGTTTGCTGACTGGCCAGCCAGCTCAAGGCATCAAAAGCGTCTTTGAGAGGGAAAGGAATATCGATAACCCGAATGCCGGGTGTCGCCGGAATCTCTTGCGACAAATGGCGCAGTAGATTTTCCAGCGCCGTAGTA  >  NC_000913/2380366‑2380632

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