Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,648,458	T→G	12.4%	N1277T (AAC→ACC)	yfhM ←	alpha2‑macroglobulin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,648,458	0	T	G	12.4%	147.6 / 15.1	57	N1277T (AAC→ACC)	yfhM	alpha2‑macroglobulin
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base G (4/3);  total (29/28)
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.80e-01

CCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGG  >  NC_000913/2648314‑2648605                                                                                                                                                 |                                                                                                                                                     cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACAGGATCGACATCATGCCCGGATCTTTTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                                 <  2:1247574/150‑1 (MQ=255)  ccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATgag                                                                                                                                               <  2:847482/150‑3 (MQ=255)     gCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGt                                                                                                                                            >  1:1336927/1‑150 (MQ=255)      cGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTc                                                                                                                                           <  2:1176589/150‑1 (MQ=255)      cGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTc                                                                                                                                           <  1:857992/150‑1 (MQ=255)         aGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                        <  1:1095564/150‑2 (MQ=255)            cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCgg                                                                                                                                     <  1:429569/150‑1 (MQ=255)             ttttGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGc                                                                                                                                    <  1:395862/150‑1 (MQ=255)                    gACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGt                                                                                                                             >  2:120321/1‑150 (MQ=255)                             acacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGtt                                                                                                                    >  1:867116/1‑150 (MQ=255)                              cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                   >  1:1108459/1‑150 (MQ=255)                                ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                                 <  2:180557/150‑1 (MQ=255)                                      cagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTGGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTgcgc                                                                                                           <  1:756345/150‑1 (MQ=255)                                       agcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGcgcg                                                                                                          <  1:44749/150‑1 (MQ=255)                                       agcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGcgcg                                                                                                          <  1:1178545/150‑1 (MQ=255)                                        gcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCgg                                                                                                         >  2:304113/1‑150 (MQ=255)                                        gcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCgg                                                                                                         >  2:333680/1‑150 (MQ=255)                                            aaGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCa                                                                                                     <  1:798032/150‑1 (MQ=255)                                             aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                    <  2:1336927/150‑1 (MQ=255)                                                    aCGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATcc                                                                                             <  2:478322/150‑1 (MQ=255)                                                              tACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAg                                                                                   <  1:313969/150‑1 (MQ=255)                                                                 tGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCg                                                                                >  2:1212599/1‑150 (MQ=255)                                                                  ggCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGt                                                                               >  2:886803/1‑150 (MQ=255)                                                                    cTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACACCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCa                                                                             >  2:1279284/1‑150 (MQ=255)                                                                          gaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGt                                                                       <  1:574591/150‑1 (MQ=255)                                                                           aTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTa                                                                      <  1:686516/150‑1 (MQ=255)                                                                                 agcgTAACGGATAGACATAAAGCCAGGATCTTGTAAAAACCCCAGCAGACGCTAATAACCCCGAAAAAAG‑GCGTCTGTAGGCACGCTGTAACCCTGTTCGCCAGCGCGGACCAGGAAATACATCACGTAAGCCGTCAGCCAGTACTAtta                                                                >  2:609039/2‑149 (MQ=255)                                                                                       cGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAcc                                                          <  1:120321/150‑1 (MQ=255)                                                                                         ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                        >  2:1252962/1‑150 (MQ=255)                                                                                           aTCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCATCAGCCAGTACTCTTCGTCACCGttt                                                      >  1:892656/1‑150 (MQ=255)                                                                                            tCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGtttt                                                     <  2:1108459/150‑1 (MQ=255)                                                                                             cGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGttttt                                                    >  2:1280795/1‑150 (MQ=255)                                                                                              gACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTa                                                   >  1:465807/1‑150 (MQ=255)                                                                                              gACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTa                                                   >  2:524335/1‑150 (MQ=255)                                                                                              gACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTa                                                   >  1:804229/1‑150 (MQ=255)                                                                                                   catGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCca                                              >  2:270411/1‑150 (MQ=255)                                                                                                     tGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcaca                                            <  1:38319/150‑1 (MQ=255)                                                                                                       cccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTTGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgc                                          >  1:341678/1‑150 (MQ=255)                                                                                                             tCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAg                                    <  1:524335/150‑1 (MQ=255)                                                                                                             tCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAg                                    <  2:319939/150‑1 (MQ=255)                                                                                                              cTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgc                                   >  1:593166/1‑150 (MQ=255)                                                                                                               ttGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgcc                                  >  1:599695/1‑150 (MQ=255)                                                                                                                   aaaTAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGccgccg                              >  1:1201297/1‑150 (MQ=255)                                                                                                                       aGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTAt                          >  1:1268239/1‑150 (MQ=255)                                                                                                                        gcgcAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATc                         >  1:1290298/1‑150 (MQ=255)                                                                                                                        gcgcAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATc                         <  2:833017/150‑1 (MQ=255)                                                                                                                         gccgacagACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCa                        >  1:52092/6‑150 (MQ=38)                                                                                                                            gacagACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGt                     >  2:532841/3‑150 (MQ=255)                                                                                                                            gacagACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACCTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGt                     >  2:1387639/3‑150 (MQ=255)                                                                                                                              cagACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTg                   <  2:599695/150‑1 (MQ=255)                                                                                                                              cagACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTg                   >  2:1293285/1‑150 (MQ=255)                                                                                                                               agACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCCTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTga                  >  1:866815/1‑149 (MQ=255)                                                                                                                                gACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa                 <  2:777585/150‑1 (MQ=255)                                                                                                                                 gacgccATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa                 >  2:1002858/6‑150 (MQ=255)                                                                                                                                     tCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc            <  2:1268239/150‑1 (MQ=255)                                                                                                                                      cATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca           <  2:867116/150‑1 (MQ=255)                                                                                                                                               ggTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACggg  <  2:1001988/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGG  >  NC_000913/2648314‑2648605

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