Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,590,960	A→G	intergenic (‑94/‑287)	ypfM ← / → yffB	uncharacterized protein YpfM/putative reductase YffB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,590,960	0	A	G	100.0%	245.6 / NA	65	intergenic (‑94/‑287)	ypfM/yffB	uncharacterized protein YpfM/putative reductase YffB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/47);  total (18/47)

CTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATACAGACAGATAACAAGATGAATATTC  >  NC_000913/2590824‑2591096                                                                                                                                         |                                                                                                                                         cTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCaa                                                                                                                                       <  2:56799/140‑1 (MQ=255)     gATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCt                                                                                                                                   <  2:19361/140‑1 (MQ=255)      aTAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCtt                                                                                                                                  >  2:78575/1‑140 (MQ=255)        aaaGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTAt                                                                                                                                <  2:1018672/140‑1 (MQ=255)            tCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTaa                                                                                                                            >  1:1203973/1‑140 (MQ=255)             cTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTaaa                                                                                                                           <  2:487406/140‑1 (MQ=255)                  cAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTa                                                                                                                      <  2:957644/140‑1 (MQ=255)                  cAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTa                                                                                                                      <  2:566988/140‑1 (MQ=255)                  cAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTa                                                                                                                      <  1:937001/140‑1 (MQ=255)                  cAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTa                                                                                                                      <  1:853732/140‑1 (MQ=255)                  cAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTa                                                                                                                      <  1:1196738/140‑1 (MQ=255)                  cAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATACACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTa                                                                                                                      <  2:821822/140‑1 (MQ=255)                     ttCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATgg                                                                                                                   >  1:956273/1‑140 (MQ=255)                         ccAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCt                                                                                                               <  2:569806/140‑1 (MQ=255)                          cAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTg                                                                                                              <  1:1130920/140‑1 (MQ=255)                          cAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTg                                                                                                              <  2:797442/140‑1 (MQ=255)                            gTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCg                                                                                                            <  1:215163/140‑1 (MQ=255)                              tCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGt                                                                                                          <  1:692617/140‑1 (MQ=255)                               cACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGtg                                                                                                         <  1:120972/140‑1 (MQ=255)                               cACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGtg                                                                                                         <  1:1017204/140‑1 (MQ=255)                                aCGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGtga                                                                                                        <  2:128927/140‑1 (MQ=255)                                aCGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGtga                                                                                                        <  1:1480717/140‑1 (MQ=255)                                aCGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGtga                                                                                                        <  1:195412/140‑1 (MQ=255)                                 cGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGtgat                                                                                                       <  2:421552/140‑1 (MQ=255)                                 cGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGtgat                                                                                                       >  2:1175127/1‑140 (MQ=255)                                  gTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGtgatg                                                                                                      >  1:203905/1‑140 (MQ=255)                                  gTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGtgatg                                                                                                      >  2:1203320/1‑140 (MQ=255)                                    tCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTgatgat                                                                                                    <  1:1287645/140‑1 (MQ=255)                                      aaTCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGAttt                                                                                                  <  1:1433081/140‑1 (MQ=255)                                       aTCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTg                                                                                                 >  1:1138143/1‑140 (MQ=255)                                        tCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGt                                                                                                <  2:789741/140‑1 (MQ=255)                                        tCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGt                                                                                                <  1:2814/140‑1 (MQ=255)                                           tAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTAtg                                                                                             >  1:747829/1‑140 (MQ=255)                                           tAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTAtg                                                                                             <  2:48676/140‑1 (MQ=255)                                                acCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATa                                                                                        <  1:643006/140‑1 (MQ=255)                                                  ctctctTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATaca                                                                                      <  1:865254/140‑1 (MQ=255)                                                        tataATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCaa                                                                                <  1:1347630/140‑1 (MQ=255)                                                                tGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCga                                                                        <  1:97845/140‑1 (MQ=255)                                                                tGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCga                                                                        <  1:1163284/140‑1 (MQ=255)                                                                tGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCga                                                                        <  2:917197/140‑1 (MQ=255)                                                                   gTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGat                                                                     <  1:524147/140‑1 (MQ=255)                                                                          aCGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATa                                                              >  2:924751/1‑140 (MQ=255)                                                                           cGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATAc                                                             >  1:315343/1‑140 (MQ=255)                                                                            ggAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACg                                                            <  1:78575/140‑1 (MQ=255)                                                                            ggAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACg                                                            <  2:1138143/140‑1 (MQ=255)                                                                            ggAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACg                                                            <  1:1175127/140‑1 (MQ=255)                                                                            ggAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACg                                                            <  1:1203320/140‑1 (MQ=255)                                                                            ggAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACg                                                            <  1:260069/140‑1 (MQ=255)                                                                                    tataCCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACg                                                            <  2:928482/132‑1 (MQ=255)                                                                                    tataCCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACg                                                            >  1:928482/1‑132 (MQ=255)                                                                                       aCCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATAc                                                 <  2:203905/140‑1 (MQ=255)                                                                                       aCCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATAc                                                 <  2:956273/140‑1 (MQ=255)                                                                                         cGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATa                                                              >  2:652455/1‑125 (MQ=255)                                                                                         cGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATa                                                              <  1:652455/125‑1 (MQ=255)                                                                                          gtggtgAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGa                                              >  2:1065047/1‑140 (MQ=255)                                                                                             gtgAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGAttt                                           >  1:162829/1‑140 (MQ=255)                                                                                              tgAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTa                                          >  2:730048/1‑140 (MQ=255)                                                                                                 aGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAAtt                                       <  1:1464240/140‑1 (MQ=255)                                                                                                     aaTTTACTCGATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTgg                                   >  2:556063/1‑140 (MQ=255)                                                                                                              gATTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATac                          <  2:747829/140‑1 (MQ=255)                                                                                                               aTTGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATaca                         <  1:197900/140‑1 (MQ=255)                                                                                                                ttGCTGCGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATacag                        >  2:1298137/1‑140 (MQ=255)                                                                                                                     gcGATGACTACCTCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATAcagacaga                   >  1:529832/1‑140 (MQ=255)                                                                                                                                 tCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATACAGACAGATAACAAGATGaa       <  1:556063/140‑1 (MQ=255)                                                                                                                                 tCCGGGGGGCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATACAGACAGATAACAAGATGaa       <  1:1065047/140‑1 (MQ=255)                                                                                                                                      ggggCAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATACAGACAGATAACAAGATGAATATTc  <  2:1203973/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         CTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATACAGACAGATAACAAGATGAATATTC  >  NC_000913/2590824‑2591096

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