Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	2,177,051	T→A	8.9%	Q43H (CAA→CAT)	U712_10890 ←	Methylglyoxal synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,177,051	0	T	A	8.9%	152.0 / 2.3	62	Q43H (CAA→CAT)	U712_10890	Methylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base T (21/33);  new base A (4/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCGGTCAGCGGGTCGCGCAAAAAGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAACCATATCCTGTTTTTTCTTGTCATGCGCGATCAAAGCAATTTTCATCATTTATCCCC  >  CP006881/2176923‑2177187                                                                                                                                 |                                                                                                                                         tGCGCGGTCAGCGGGTCGCGCAAAAAGATGACAAGATCGAGTGCATTGGCCGCGCTCAGTGCCCCGATGTGCTGGTCTCCCCCGAAAGGGCCGGATTGAAAACGGTCAATTTGAAGACCTGTCGCGACGTGAAatttg                                                                                                                                 >  4:40952/1‑137 (MQ=255) tGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATCGAGTGTATTGGCAGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCATGAATTTTc                                                                                                                                 >  3:87994/1‑138 (MQ=255)   cgcgGTCAGCGGGTCGCGCAAAAAGATGACAAGATCGAGTGGATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAAc                                                                                                                              >  5:349706/1‑139 (MQ=255)       gTCAGCGGGTCGCGCAAAAAAATGACAAGATCGAGTGTATTGGCAGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCATGAATTTTCAACCCtg                                                                                                                          <  4:87994/139‑1 (MQ=255)             gggTCGCGCAAAAAGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGc                                                                                                                    <  5:216021/139‑1 (MQ=255)             gggTCGCGCAAAAAGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGc                                                                                                                    <  4:84858/139‑1 (MQ=255)                tCGCGCAAAAAGATGCCAAGATCGAGTGCAGTGGCCGCGATCAGTGCACCGATTTGCTCGTCTCCCCCTAAACGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCgg                                                                                                                 <  2:142480/139‑1 (MQ=255)                tCGCGCAAAAAGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCgg                                                                                                                 <  7:113104/139‑1 (MQ=255)                  gcgcAAAAAGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTGAACCCTGTTGTGCCgggt                                                                                                               >  1:37785/1‑137 (MQ=255)                  gcgcAAAAAGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGtt                                                                                                               <  6:155408/139‑1 (MQ=255)                  gcgcAAAAAGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGtt                                                                                                               <  3:340816/139‑1 (MQ=255)                       aaaaGCTGACAAGATCGAGTGCCTTGGCCGCGCTCAGTGCACCGAGTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGta                                                                                                          <  6:47179/139‑1 (MQ=255)                       aaaaGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCTGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGta                                                                                                          <  6:95320/139‑1 (MQ=255)                       aaaaGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGta                                                                                                          <  4:314644/139‑1 (MQ=255)                       aaaaGATGACAAGATCGAGTGCATTGGCCGAGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGta                                                                                                          <  1:128958/139‑1 (MQ=255)                         aaGATGACAAGATCGAGTGCATTGGCCGAGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGtata                                                                                                        <  6:349706/139‑1 (MQ=255)                                     tCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGcc                                                                                            <  5:348080/139‑1 (MQ=255)                                      cGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAAcc                                                                                                                             >  3:21678/1‑105 (MQ=255)                                      cGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAAcc                                                                                                                             <  4:21678/105‑1 (MQ=255)                                      cGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCGCGGTAAATTTACACCCTTTTGGTCCGGTTGCGGTAAGCTCAAGATGccg                                                                                           >  2:176650/1‑138 (MQ=255)                                       gAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGcc                                                                                            <  2:348951/137‑1 (MQ=255)                                        aGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCg                                                                                                                  >  7:286998/1‑114 (MQ=255)                                        aGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCg                                                                                                                  <  8:286998/114‑1 (MQ=255)                                             cttGGCCGCTATCAGGGCACCGATTTGCTGGTCGCCCCCTAAAGGGCCGCATTGAACACGTTCAATTTGAAGACCTGGCGACTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAata                                                                                    <  2:54454/138‑1 (MQ=255)                                             aTTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAata                                                                                    <  5:294812/139‑1 (MQ=255)                                             aTTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAata                                                                                    <  8:258522/139‑1 (MQ=255)                                               tGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGCCCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATc                                                                                  >  8:224357/1‑139 (MQ=255)                                               tGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTGAACCCAGTTGTGCCGGTTGCGGAATGATCATTGTGCCACAAAatata                                                                                  >  4:2029/1‑138 (MQ=255)                                               tGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATc                                                                                  >  3:139385/1‑139 (MQ=255)                                                       aTCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAgg                                                                          <  7:20029/139‑1 (MQ=255)                                                       aTCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAgg                                                                          <  7:224357/139‑1 (MQ=255)                                                        tCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAgg                                                                          <  2:160136/138‑1 (MQ=255)                                                           gTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCATGAATTTTCAACCCTGTGGTTCCGGTTGCGTATAGATCATGATTCTTTAAAATATCCCGATAGGCAga                                                                      >  5:113686/1‑138 (MQ=255)                                                            tGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTc                                                                     >  1:101155/1‑139 (MQ=255)                                                            tGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTc                                                                     <  3:322355/139‑1 (MQ=255)                                                                 cGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGtata                                                                                                        <  1:129306/99‑1 (MQ=255)                                                                 cGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGtata                                                                                                        >  2:129306/1‑99 (MQ=255)                                                                 cGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGt                                                                                                           <  8:192377/96‑1 (MQ=255)                                                                 cGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGt                                                                                                           >  7:192377/1‑96 (MQ=255)                                                                   aTTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAAtt                                                              <  5:114185/139‑1 (MQ=255)                                                                   aTTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAAtt                                                              <  6:79927/139‑1 (MQ=255)                                                                      tGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCGTCTTGAATTTTAAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGaa                                                           >  3:83487/1‑139 (MQ=255)                                                                      tGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTAAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGaa                                                           >  8:241318/1‑139 (MQ=255)                                                                      tGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTAAACCCTGTTGTGCCGGTTGCGTATAGATAATGATTCCTCAAAATATCCCGATAGGCAGTCGTAAAtttgt                                                           >  7:145661/1‑136 (MQ=255)                                                                       gCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCATGAATTTTCAACCCTGTGGTTCCGGTTGCGTATAGATCATGATTCTTTAAAATATCCCGATAGGCAGTCGTAAATTGaa                                                           <  6:113686/138‑1 (MQ=255)                                                                            tctcCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAACCatat                                                     <  3:40952/139‑1 (MQ=255)                                                                             ctcCCCCTCAAGGGCCGGATTGAAAACGTTCAATTAGAAGACCTGGCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAACCATATc                                                    <  6:61249/139‑1 (MQ=255)                                                                             ctcCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTTAACCataaa                                                    >  4:140415/1‑137 (MQ=255)                                                                             ctcCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAACCATATc                                                    >  1:31358/1‑139 (MQ=255)                                                                              tcCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCATGAATTTTCAACCCTGTGGTTGCGGTTGCGGATAGATCATGATTCGTTTAAATATCCCGGTAGGCAGTGGGTAATTGAACCAGATcc                                                   >  3:214434/1‑139 (MQ=255)                                                                                ccccTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAACCATATc                                                    <  7:241318/136‑1 (MQ=255)                                                                                     aaaGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAAtt                                                              >  6:246478/1‑121 (MQ=255)                                                                                     aaaGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAAtt                                                              <  5:246478/121‑1 (MQ=255)                                                                                          gCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCGGCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTGGTAAATTGAACCATATCCTGGTGTttgttg                                       >  4:43708/1‑139 (MQ=255)                                                                                          gCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAAcc                                                         <  8:288884/121‑1 (MQ=255)                                                                                          gCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAAcc                                                         >  7:288884/1‑121 (MQ=255)                                                                                          gCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTAAACCCTGTTGTGCGCGTTGCGTATTGATCATGATTGCTCAAAATATCCCGTGAGGGACTCCTTAAATGCACCAAAACGAGTttttttttt                                       >  2:93347/1‑138 (MQ=255)                                                                                            cGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATATATCATGATGCCTCAAAATATACCGATAGGCAGTCGTAAATTGAACCATATACTGTTTTTtgttgt                                      >  6:60304/1‑138 (MQ=255)                                                                                            cGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCATGGATTTTGACCCCTTTTGTGCCGGTTGCGTTTAGATCATTATGGCTCAAAATATCCCGATAGGAAGTCGTAAATTGAACCATATGCTGTTTTTTCATGt                                      >  8:17811/1‑138 (MQ=255)                                                                                                                  aaGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATc                                                                                  <  3:151539/72‑1 (MQ=255)                                                                                                                  aaGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATc                                                                                  >  4:151539/1‑72 (MQ=255)                                                                                                                    gACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAACCATATCCTGTTTTTTCTTGTCATGCGCGATCAAAGCAATTTtcat             <  2:120781/139‑1 (MQ=255)                                                                                                                          tagcatctTGAATTTTCAACCCTGTTGTGCCGGTTGAGTATAGATCATGATGCATCAAAATATCCCGATAGGCAGTCGTAAATTGAACCATATCCTGTTTTTTCTTGTCATGCGCGATCAAAGCAATTTTCATCATTTa       <  7:201002/134‑1 (MQ=255)                                                                                                                                 ttGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAACCATATCCTGTTTTTTCTTGTCATGCGCGATCAAAGCAATTTTCATCATTTATcccc  <  1:93347/137‑1 (MQ=255)                                                                                                                                 |                                                                                                                                         TGCGCGGTCAGCGGGTCGCGCAAAAAGATGACAAGATCGAGTGCATTGGCCGCGATCAGTGCACCGATTTGCTGGTCTCCCCCTAAAGGGCCGGATTGAAAACGTTCAATTTGAAGACCTGTCGCCTCTTGAATTTTCAACCCTGTTGTGCCGGTTGCGTATAGATCATGATGCCTCAAAATATCCCGATAGGCAGTCGTAAATTGAACCATATCCTGTTTTTTCTTGTCATGCGCGATCAAAGCAATTTTCATCATTTATCCCC  >  CP006881/2176923‑2177187

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑