Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	237,539	G→C	6.3%	A69P (GCG→CCG) ‡	yafT →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	237,539	0	G	C	6.3%	213.9 / 3.7	63	A69P (GCG→CCG) ‡	yafT	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base C (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTGGTATCCCG  >  NC_000913/237399‑237676                                                                                                                                             |                                                                                                                                          aaTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGAt                                                                                                                                   <  4:201353/149‑1 (MQ=255)  aTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATg                                                                                                                                  <  2:441653/149‑1 (MQ=255)             atgaAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTc                                                                                                                       <  1:223460/149‑1 (MQ=255)                aaaaGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCcgc                                                                                                                    >  2:402561/1‑149 (MQ=255)                aaaaGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCcgc                                                                                                                    >  2:321220/1‑149 (MQ=255)                    gaagaaAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTc                                                                                                                >  3:216908/1‑149 (MQ=255)                      agaaAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                                                                                                              <  3:73657/149‑1 (MQ=255)                       gaaAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTcccc                                                                                                             >  3:80103/1‑149 (MQ=255)                         aaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCtt                                                                                                           >  3:55628/1‑149 (MQ=255)                          aaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                                                                                                              >  1:271954/1‑145 (MQ=255)                          aaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                                                                                                              <  2:271954/145‑1 (MQ=255)                          aaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCGCAGCTTGaata                                                                                                                                 >  4:341849/1‑123 (MQ=255)                          aaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCGCAGCTTGaata                                                                                                                                 <  3:341849/126‑4 (MQ=255)                                  aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCt                                                                                                  <  4:207441/149‑1 (MQ=255)                                  aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCGCAGCTTGaa                                                                                                                                   >  2:263536/1‑115 (MQ=255)                                  aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCGCAGCTTGaa                                                                                                                                   <  1:263536/116‑2 (MQ=255)                                        cacaAAGTAATATTAATAAAAACACTCCCAAGCAACTGACAGACAAAGATTTACTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATa                                                                                            <  2:269787/149‑1 (MQ=255)                                         acaAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATAc                                                                                           <  1:321220/149‑1 (MQ=255)                                         acaAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGAGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATAc                                                                                           >  3:299442/1‑149 (MQ=255)                                           aaaGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACtt                                                                                         >  3:63915/1‑149 (MQ=255)                                                 atatTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                                                                                                   >  3:456264/1‑133 (MQ=255)                                                 atatTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                                                                                                   <  4:456264/133‑1 (MQ=255)                                                    ttAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCcgc                                                                                                                    <  4:200119/113‑1 (MQ=255)                                                    ttAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCcgc                                                                                                                    >  3:200119/1‑113 (MQ=255)                                                        taaAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCa                                                                            >  1:410841/1‑149 (MQ=255)                                                          aaaaCACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTcc                                                                                                                      >  1:319102/1‑105 (MQ=255)                                                          aaaaCACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTcc                                                                                                                      <  2:319102/105‑1 (MQ=255)                                                           aaaCACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAAcc                                                                         >  3:11023/1‑149 (MQ=255)                                                              cacTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACcgcg                                                                      >  3:162346/1‑149 (MQ=255)                                                                   cccAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTaa                                                                                              >  3:219472/1‑120 (MQ=255)                                                                   cccAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTaa                                                                                              <  4:219472/120‑1 (MQ=255)                                                                     cAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGa                                                               >  2:209676/1‑149 (MQ=255)                                                                          aCTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCAc                                                                   <  2:335230/140‑1 (MQ=255)                                                                          aCTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCAc                                                                   >  1:335230/1‑140 (MQ=255)                                                                             gacagacaAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTa                                                       <  4:162346/149‑1 (MQ=255)                                                                              acagacaAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTAt                                                      <  4:55628/149‑1 (MQ=255)                                                                                 gacaAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCa                                                   >  4:323201/1‑149 (MQ=255)                                                                                   caAAGATTTATTCAGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCagg                                                 <  4:200853/149‑1 (MQ=255)                                                                                    aaaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCg                                                                               <  2:423149/118‑1 (MQ=255)                                                                                    aaaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCg                                                                               >  1:423149/1‑118 (MQ=255)                                                                                    aaaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCagga                                                >  4:282413/1‑149 (MQ=255)                                                                                     aaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                                                                                                   >  3:25684/1‑97 (MQ=255)                                                                                     aaGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                                                                                                   <  4:25684/97‑1 (MQ=255)                                                                                       gATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAggagga                                             <  4:80103/149‑1 (MQ=255)                                                                                         tttATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGgaga                                           >  2:375225/1‑149 (MQ=255)                                                                                          ttATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGAt                                          <  1:209676/149‑1 (MQ=255)                                                                                          ttATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATAAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTAAGCAGGAGGAGAt                                          <  4:299442/149‑1 (MQ=255)                                                                                             ttCGGTAATGAAACCACTCTGGCTGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGTAAAc                                                            >  1:449575/1‑128 (MQ=255)                                                                                             ttCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCg                                       >  1:147833/1‑149 (MQ=255)                                                                                             ttCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCCCACCGGAAAc                                                            <  2:449575/128‑1 (MQ=255)                                                                                                     tGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTGGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAAtatt                               >  1:416439/1‑149 (MQ=255)                                                                                                      gAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAAtatt                               >  4:22649/1‑148 (MQ=255)                                                                                                      gAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAAtatt                               <  3:22649/148‑1 (MQ=255)                                                                                                                  gCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCcac                  >  3:455174/1‑149 (MQ=255)                                                                                                                      tATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATtc              >  2:94951/1‑149 (MQ=255)                                                                                                                       aTCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATtct             >  1:214382/1‑149 (MQ=255)                                                                                                                        tCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATAGTGTTTCCACATtctc            <  4:455174/149‑1 (MQ=255)                                                                                                                          cGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTg          <  1:402561/149‑1 (MQ=255)                                                                                                                           gAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTgg         >  2:310364/1‑149 (MQ=255)                                                                                                                            aaGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTGGt        <  4:216908/149‑1 (MQ=255)                                                                                                                                 gATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGt                                                                                       <  4:211097/65‑1 (MQ=255)                                                                                                                                 gATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGt                                                                                       >  3:211097/1‑65 (MQ=255)                                                                                                                                  atatTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTGGTATCCCg  <  1:375225/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                          AATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTCAATCCGGCAACCGCGCACCGGAAACCATTATGCAGGAGGAGATGCGAAAATATTATACTGTTTCCACATTCTCTGGTATCCCG  >  NC_000913/237399‑237676

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