Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,984,411	C→T	*394* (TAG→TAA)	yqeF ←	short chain acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,411	0	C	T	100.0%	240.4 / NA	76	*394* (TAG→TAA)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (39/37);  total (39/37)

AGTGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAG  >  NC_000913/2984275‑2984538                                                                                                                                         |                                                                                                                                aGTGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTAACGGTAAGGAGGAAATGCAGAGTTa                                                                                                                               >  2:1291999/1‑139 (MQ=255)             aTCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTGGTCACGTTc                                                                                                                   >  2:1173825/1‑139 (MQ=255)             aTCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTc                                                                                                                   >  2:679947/1‑139 (MQ=255)              tCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCa                                                                                                                  <  2:1525140/139‑1 (MQ=255)                   aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  1:1173825/139‑1 (MQ=255)                            tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCg                                                                                                                           <  2:268109/116‑1 (MQ=255)                            tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCg                                                                                                                           >  1:268109/1‑116 (MQ=255)                            tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCcac                                                                                                    >  1:1749589/1‑139 (MQ=255)                                aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                >  2:925788/1‑139 (MQ=255)                                  tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                              >  1:1537505/1‑139 (MQ=255)                                  tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                              >  1:491764/1‑139 (MQ=255)                                  tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                              >  1:1784468/1‑139 (MQ=255)                                         aTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccg                                                                                       <  1:925788/139‑1 (MQ=255)                                          tGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                      >  1:1403279/1‑139 (MQ=255)                                                          atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGc                                                                      <  2:1784468/139‑1 (MQ=255)                                                               tgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      >  1:130813/1‑102 (MQ=255)                                                               tgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                 >  1:867733/1‑139 (MQ=255)                                                               tgtCGACTTGCCAGCCAACATTATTTAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      <  2:130813/102‑1 (MQ=255)                                                                 tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                               >  2:511417/1‑139 (MQ=255)                                                                 tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                               >  2:342949/1‑139 (MQ=255)                                                                   gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCaa                                                                                                                 >  2:686877/1‑87 (MQ=255)                                                                   gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCaa                                                                                                                 <  1:686877/87‑1 (MQ=255)                                                                     cTTTCCAGCCATCATTATCCAACAGATTAATGACGGCGGATATTTATCGGTAAGGAGAAAATGAAGAGTTATTCGTCACCTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                           <  1:1291999/139‑1 (MQ=255)                                                                     cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                           <  1:342949/139‑1 (MQ=255)                                                                     cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                           <  1:511417/139‑1 (MQ=255)                                                                      ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      <  2:1535817/95‑1 (MQ=255)                                                                      ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      >  1:1535817/1‑95 (MQ=255)                                                                       tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                >  1:1228708/1‑100 (MQ=255)                                                                       tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                <  2:1228708/100‑1 (MQ=255)                                                                           agccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                     <  2:491764/139‑1 (MQ=255)                                                                             ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                           <  2:145680/83‑1 (MQ=255)                                                                             ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                           >  1:145680/1‑83 (MQ=255)                                                                               atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                 <  2:867733/139‑1 (MQ=255)                                                                                      ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACGCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                             <  2:971289/136‑1 (MQ=255)                                                                                      ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                             >  1:971289/1‑136 (MQ=255)                                                                                      ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                          >  1:608840/1‑139 (MQ=255)                                                                                      ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                          <  2:1777196/139‑1 (MQ=255)                                                                                       tCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                         <  2:1403279/139‑1 (MQ=255)                                                                                           cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                     >  2:1480194/1‑139 (MQ=255)                                                                                            aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                    <  2:1749589/139‑1 (MQ=255)                                                                                             gATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             >  2:118856/1‑65 (MQ=255)                                                                                             gATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  1:118856/65‑1 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                         >  2:294271/1‑52 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                         <  2:942388/52‑1 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                         <  1:294271/52‑1 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                         >  1:942388/1‑52 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             >  1:104258/1‑64 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  2:104258/64‑1 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      >  2:1645653/1‑71 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      <  1:1645653/71‑1 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                           <  2:1394236/98‑1 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                           >  1:1394236/1‑98 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                           >  2:1631835/1‑114 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                           <  1:1631835/114‑1 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                  >  2:174398/1‑139 (MQ=255)                                                                                              aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                  >  1:399363/1‑139 (MQ=255)                                                                                                tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      >  2:151821/1‑69 (MQ=255)                                                                                                tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      <  1:151821/69‑1 (MQ=255)                                                                                                tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                <  1:174398/139‑1 (MQ=255)                                                                                                       ggCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                         >  1:1245248/1‑139 (MQ=255)                                                                                                        gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                           >  2:819510/1‑56 (MQ=255)                                                                                                        gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                           <  1:819510/56‑1 (MQ=255)                                                                                                        gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                         <  2:1189646/122‑1 (MQ=255)                                                                                                        gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                         >  1:1189646/1‑122 (MQ=255)                                                                                                                    cGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                        >  1:1119831/1‑95 (MQ=255)                                                                                                                    cGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                        <  2:1119831/95‑1 (MQ=255)                                                                                                                       tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                         >  1:992783/1‑107 (MQ=255)                                                                                                                       tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa         <  2:608840/139‑1 (MQ=255)                                                                                                                       tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAAGGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                         <  2:992783/107‑1 (MQ=255)                                                                                                                         aggaggAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGt                                                                  <  2:1187173/80‑1 (MQ=255)                                                                                                                         aggaggAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGt                                                                  >  1:1187173/1‑80 (MQ=255)                                                                                                                          ggaggaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg      >  2:75005/1‑139 (MQ=255)                                                                                                                              gaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg  <  2:1537505/139‑1 (MQ=255)                                                                                                                                      gagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           >  2:469931/1‑42 (MQ=255)                                                                                                                                      gagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           <  1:469931/42‑1 (MQ=255)                                                                                                                                      gagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                     <  2:749545/80‑1 (MQ=255)                                                                                                                                         |                                                                                                                                AGTGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAG  >  NC_000913/2984275‑2984538

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