Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,982,684	(C)8→9	noncoding (79/87 nt)	leuT →	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,982,676	1	.	C	91.7%	195.4 / 8.7	60	noncoding (71/87 nt)	leuT	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base C (27/28);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGA  >  NC_000913/3982544‑3982793                                                                                                                                      |                                                                                                                        tgATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑cccc‑cc                                                                                                                 <  2:1492081/139‑1 (MQ=255)  gATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑cccc‑ccc                                                                                                                <  1:1009326/139‑1 (MQ=255)      ccAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑cccc‑ccccctc                                                                                                            >  2:366304/1‑136 (MQ=255)         gTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGca                                                                                                          >  2:246490/1‑139 (MQ=255)              cGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACg                                                                                                     <  1:1878628/139‑1 (MQ=255)                    ttCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTa                                                                                               <  1:776906/139‑1 (MQ=255)                         aTCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGGACCACGACTTTAAAGaa                                                                                          <  2:1176315/139‑1 (MQ=255)                               ttAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAAc                                                                                    >  2:1795327/1‑139 (MQ=255)                                    cACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTaaaa                                                                               <  1:246490/139‑1 (MQ=255)                                       cccATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGa                                                                                                    <  1:255819/115‑1 (MQ=255)                                       cccATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGa                                                                                                    >  2:255819/1‑115 (MQ=255)                                       cccATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTc                                                                           >  1:1873259/1‑139 (MQ=255)                                        ccATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGcaccac                                                                                                      >  2:1273059/1‑112 (MQ=255)                                        ccATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTC‑CCCC‑CCCCTCGcaccac                                                                                                      <  1:1273059/112‑1 (MQ=255)                                          attattAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaa                                                                         >  2:750003/1‑139 (MQ=255)                                                     ttGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                                              >  1:356655/1‑139 (MQ=255)                                                      tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAg                                                                     >  1:30812/1‑131 (MQ=255)                                                      tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAg                                                                     <  2:30812/131‑1 (MQ=255)                                                      tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAttt                                                             >  2:938899/1‑139 (MQ=255)                                                      tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑cccc‑ccccct                                                                                                             >  1:262399/1‑88 (MQ=255)                                                      tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑cccc‑ccccct                                                                                                             <  2:262399/90‑3 (MQ=255)                                                      tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                                                            >  1:636195/1‑139 (MQ=255)                                                       gtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                                                            >  2:1531425/1‑139 (MQ=255)                                                        tgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAAt                                                                                         >  2:55156/1‑109 (MQ=255)                                                        tgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAAt                                                                                         <  1:55156/109‑1 (MQ=255)                                                        tgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg                                                           >  2:1290328/1‑139 (MQ=255)                                                        tgACAATGCGAAGGTGGCGAAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg                                                           >  2:1252197/1‑139 (MQ=255)                                                          aCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAg                                                                     <  2:888685/128‑1 (MQ=255)                                                          aCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAg                                                                     >  1:888685/1‑128 (MQ=255)                                                           cAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCg                                                        <  2:1304227/139‑1 (MQ=255)                                                               gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAg                                                   <  2:1873259/139‑1 (MQ=255)                                                                cGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGa                                                                                                    >  2:1463121/1‑90 (MQ=255)                                                                cGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGa                                                                                                    <  1:1463121/90‑1 (MQ=255)                                                                  aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGa                                                                                                    <  2:287737/88‑1 (MQ=255)                                                                  aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGa                                                                                                    >  1:287737/1‑88 (MQ=255)                                                                  aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGGCTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcg                                                 <  1:1531425/139‑1 (MQ=255)                                                                  aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcg                                                 <  1:366304/139‑1 (MQ=255)                                                                        gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCtt                                           >  2:1607127/1‑139 (MQ=255)                                                                        gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCtt                                           >  1:1367662/1‑139 (MQ=255)                                                                          ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTgg                                         <  1:1795327/139‑1 (MQ=255)                                                                             aTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAg                                      <  2:356655/139‑1 (MQ=255)                                                                                 gTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCaa                                 <  2:636195/139‑1 (MQ=255)                                                                                   aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAc                                <  1:1607127/139‑1 (MQ=255)                                                                                   aGACGATCTAGTTTCAGGTGTTAGTGTCCTTACGGTCGTGGGCGGTAAAGT‑‑CCCCACCCCTCGCACGAAGACTTTCAAGAATTGAACTAAAAATTCAAAAAGCATTATTACCGCGAGTAGCGCAGCTTGGTAGCGCAAc                                <  1:285018/139‑1 (MQ=25)                                                                                    gACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACt                               <  1:1000051/139‑1 (MQ=255)                                                                                      cgcgCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCa                                                                          >  2:334589/1‑94 (MQ=255)                                                                                      cgcgCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCa                                                                          <  1:334589/94‑1 (MQ=255)                                                                                        cgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGtt                           <  1:1290328/139‑1 (MQ=255)                                                                                          cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAg                                      <  1:201247/126‑1 (MQ=255)                                                                                          cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAg                                      >  2:201247/1‑126 (MQ=255)                                                                                           tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaa                                                                        <  2:1289747/91‑1 (MQ=255)                                                                                           tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaa                                                                        >  1:1289747/1‑91 (MQ=255)                                                                                           tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTgg                        <  1:1252197/139‑1 (MQ=255)                                                                                            aGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTggg                       >  1:404861/1‑139 (MQ=255)                                                                                                       ttAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACt                               <  1:1009110/120‑1 (MQ=255)                                                                                                       ttAGTGTCCTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACt                               >  2:1009110/1‑120 (MQ=255)                                                                                                        tAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTc            >  1:1302007/1‑139 (MQ=255)                                                                                                              ccTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGtt     <  1:938899/139‑1 (MQ=255)                                                                                                               cTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAg                                                   <  2:457562/92‑1 (MQ=255)                                                                                                               cTTACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAg                                                   >  1:457562/1‑92 (MQ=255)                                                                                                                 tACGGACGTGGGGGTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGa  >  2:602126/1‑139 (MQ=255)                                                                                                                         tGGGGGTTCAAGTCCCCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg                                                           >  2:512325/1‑75 (MQ=39)                                                                                                                         tGGGGGTTCAAGTCCCCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg                                                           <  1:512325/75‑1 (MQ=39)                                                                                                                           ggggTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaa                                                                        >  2:1189497/1‑59 (MQ=38)                                                                                                                           ggggTTCAAGTC‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaa                                                                        <  1:1189497/59‑1 (MQ=38)                                                                                                                                      |                                                                                                                        TGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCC‑CCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGA  >  NC_000913/3982544‑3982793

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