Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	864,667	A→T	7.3%	F173Y (TTC→TAC)	uraA ←	uracil permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	864,667	0	A	T	7.3%	162.1 / 2.8	55	F173Y (TTC→TAC)	uraA	uracil permease
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base T (2/2);  total (28/27)
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.75e-01

TTGAACGCAGGGGCGGTGAAATTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAGGCCGATGGAGATGATCAC  >  NC_002947/864520‑864815                                                                                                                                                     |                                                                                                                                                      ttGAACGCAGGGGCGGTGAAATTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGaa                                                                                                                                                      <  4:335514/149‑1 (MQ=255) ttGAACGCAGGGGCGGTGAAATTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGaa                                                                                                                                                      >  3:335514/1‑149 (MQ=255)    aaCGCAGGGGCGGTGAAATTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATAc                                                                                                                                                   >  6:339330/1‑149 (MQ=255)       gCAGGGGCGGTGAAATTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGAGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACgg                                                                                                                                                 >  4:306568/1‑148 (MQ=255)           gggCGGTGAAATTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAAc                                                                                                                                             >  8:269235/1‑148 (MQ=255)                gTGAAATTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCaa                                                                                                                                       <  4:3054/149‑1 (MQ=255)                 tGAAATTGGGAATTTCCAGCCACCCGGCGGC‑AGCGATCCGGGCGGTGTCGACCACGCCAACGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAg                                                                                                                                      <  8:344993/149‑1 (MQ=255)                   aaaTTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGTATTCGGCaagg                                                                                                                                            <  6:97453/141‑3 (MQ=255)                   aaaTTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGTATTCGGCaagg                                                                                                                                            >  5:97453/1‑139 (MQ=255)                     aTTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTg                                                                                                                                  <  5:76603/149‑1 (MQ=255)                      ttGGGAATTTCCAGCCACGGTGCGGCGAGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTg                                                                                                                                  >  1:213365/1‑149 (MQ=255)                      ttGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGt                                                                                                                                 <  7:195018/149‑1 (MQ=255)                       tGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGt                                                                                                                                 >  4:179702/1‑148 (MQ=255)                         ggAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGAt                                                                                                                                          <  1:338832/137‑1 (MQ=255)                         ggAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGAt                                                                                                                                          >  2:338832/1‑137 (MQ=255)                         ggAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTca                                                                                                                               <  8:291984/148‑1 (MQ=255)                                 cAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑Tg                                                                                                                     <  2:213365/149‑1 (MQ=255)                                  aGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑Tgg                                                                                                                    >  7:23711/1‑149 (MQ=255)                                       cGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCg                                                                                                               <  7:269235/149‑1 (MQ=255)                                        ggTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCga                                                                                                              <  7:191509/149‑1 (MQ=255)                                             ggC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCa                                                                                                         <  8:346397/149‑1 (MQ=255)                                               c‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGc                                                                                                       <  7:117066/149‑1 (MQ=255)                                                   cGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGAGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATg                                                                                                    <  3:306568/149‑1 (MQ=255)                                                    gATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGTATTCGGCaaggg                                                                                                                                           >  1:256383/1‑107 (MQ=255)                                                    gATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGTATTCGGCaaggg                                                                                                                                           <  2:256383/110‑4 (MQ=255)                                                    gATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATgg                                                                                                   <  8:23711/149‑1 (MQ=255)                                                                    aCCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCAt                                                                                   <  5:339330/149‑1 (MQ=255)                                                                        cGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGt                                                                               >  3:381299/1‑149 (MQ=255)                                                                        cGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCCTGTACGGCATCAGt                                                                               <  7:353800/149‑1 (MQ=255)                                                                         gCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTg                                                                              <  7:90129/149‑1 (MQ=255)                                                                           cAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGAAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGcc                                                                            >  5:188422/1‑149 (MQ=255)                                                                            aaaGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACAAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGAATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGcag                                                                           >  4:222820/1‑147 (MQ=255)                                                                                         aGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCAcc                                                              >  8:34383/1‑149 (MQ=255)                                                                                          gggCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCt                                                             >  3:214494/1‑149 (MQ=255)                                                                                           ggCAAAGCGCCCGAGCGCGCCGGAGATGATGGGTGCCCGGCGGAAAGTGCCCTTGCCGAATGCGGCGACGATCAGGGTTGTCAGCAGCGC‑TGGCGTCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCGGTCCCGTCgccgccgt                                                            >  2:2786/1‑147 (MQ=255)                                                                                           ggCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGTAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCtt                                                            >  8:21307/1‑149 (MQ=255)                                                                                            gCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCttt                                                           >  5:172100/1‑149 (MQ=255)                                                                                            gCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCttt                                                           >  3:56244/1‑149 (MQ=255)                                                                                            gCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCttt                                                           >  1:348278/1‑149 (MQ=255)                                                                                              aaaGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGc                                                         <  4:381299/149‑1 (MQ=255)                                                                                              aaaGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGc                                                         <  3:179702/149‑1 (MQ=255)                                                                                               aaGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGcc                                                        >  8:196626/1‑149 (MQ=255)                                                                                                   ccaccaGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATc                                                    <  8:308745/149‑1 (MQ=255)                                                                                                    caccaGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCg                                                   >  3:65593/1‑149 (MQ=255)                                                                                                     accaGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACgg                                                                                      >  4:390067/1‑113 (MQ=255)                                                                                                     accaGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACgg                                                                                      <  3:390067/113‑1 (MQ=255)                                                                                                           aCGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGtt                                            <  4:333646/149‑1 (MQ=255)                                                                                                                 gagaTGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTggcggc                                      >  5:126355/1‑149 (MQ=255)                                                                                                                     tgatgGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGc                                                         >  7:81262/1‑126 (MQ=255)                                                                                                                     tgatgGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGc                                                         <  8:81262/126‑1 (MQ=255)                                                                                                                     tgatgGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATc                                  <  6:172100/149‑1 (MQ=255)                                                                                                                               cAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCATGGCATCGAGATCAGAATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCAATCGCCTTGTTGGCGGCGATCGGTGCCATg                         >  1:107945/1‑149 (MQ=255)                                                                                                                               cAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATgg                        >  3:263671/1‑149 (MQ=255)                                                                                                                                 ggCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGc                       >  8:220293/1‑148 (MQ=255)                                                                                                                                  gCGGAAAATGCCCTTGCCGAATCCGGCAACGATCAAGGTTCTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCa                     <  1:2786/149‑1 (MQ=255)                                                                                                                                    ggAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAgg                   >  6:335360/1‑149 (MQ=255)                                                                                                                                                  ccGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAGGCCGATGGAgatgat     <  4:263671/149‑1 (MQ=255)                                                                                                                                                     aaTACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAGGCCGATGGAGATGATCAc  <  4:65593/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TTGAACGCAGGGGCGGTGAAATTGGGAATTTCCAGCCACGGTGCGGC‑AGCGATCCTGGCGGTGTCGACCACGCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGC‑TGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAGGCCGATGGAGATGATCAC  >  NC_002947/864520‑864815

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