Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,597,313:1	(C)7→8	100%	pseudogene (675/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,306	1	.	C	97.1%	243.8 / 0.5	70	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (37/31);  total (38/32)
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.87e-01

AGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597207‑1597400                                                                                                     |                                                                                                   agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa                                                                                                     >  1:1466869/1‑100 (MQ=255)    gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑‑‑‑cc                                                                                               <  1:1964615/99‑1 (MQ=255)    gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑‑‑cc‑c                                                                                             <  1:2169390/100‑1 (MQ=255)      cTGCGTACCGTGATGACAGAGTCTTGATCCTGTACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑‑‑cc‑ccc                                                                                           <  1:2380777/100‑1 (MQ=255)            aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑‑‑cc‑cccccctcg                                                                                     >  1:512685/1‑96 (MQ=255)            aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑‑‑cc‑cccccctcg                                                                                     >  1:687562/1‑96 (MQ=255)              cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTcgcg                                                                                    <  2:343791/100‑1 (MQ=255)                tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGaa                                                                                  >  1:2274825/1‑100 (MQ=255)                 gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAAt                                                                                 <  2:1224190/100‑1 (MQ=255)                  atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATa                                                                                >  1:1922732/1‑100 (MQ=255)                   tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAg                                                                               <  2:462429/100‑1 (MQ=255)                    gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGt                                                                              >  1:1771229/1‑100 (MQ=255)                    gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGt                                                                              >  2:454535/1‑100 (MQ=255)                      cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCg                                                                            >  1:922308/1‑100 (MQ=255)                       agagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGc                                                                           <  1:1065487/100‑1 (MQ=255)                           tCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTcc                                                                       <  1:868850/100‑1 (MQ=255)                              tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACC‑‑CC‑CCCCCTCGCGAATAGTCGCCTcccc                                                                     >  2:1752551/1‑100 (MQ=255)                               gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAg                                                                   <  1:908513/100‑1 (MQ=255)                                 tCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACC‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGt                                                                  >  1:1574891/1‑100 (MQ=255)                                  ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTc                                                                >  2:1881724/1‑100 (MQ=255)                                   cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCCGTTcc                                                               >  1:1012461/1‑100 (MQ=255)                                   cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                               >  1:2266464/1‑100 (MQ=255)                                    tGCACATGCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑CCCCCTCGCGAATAATCGCCTCCCCAGtt                                                                 >  1:1248042/1‑100 (MQ=255)                                     gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                             >  1:1749271/1‑100 (MQ=255)                                      cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            <  1:1788373/100‑1 (MQ=255)                                      cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            <  1:307969/100‑1 (MQ=255)                                      caaaTTCAAGGTTCCGATACCGGTGGCATTGCCACCGAAAATCGTTTTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                              <  1:1349172/95‑1 (MQ=255)                                         aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         <  2:45038/100‑1 (MQ=255)                                         aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         <  1:79020/100‑1 (MQ=255)                                         aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         >  2:920618/1‑100 (MQ=255)                                          ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                        <  2:354215/100‑1 (MQ=255)                                               ggTTCCGATACCGGTGGCATTGCCACCGATAATCTTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                                                   <  1:1592188/100‑1 (MQ=255)                                               ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCCAt                                                   >  2:1935495/1‑100 (MQ=255)                                               ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTAGCCTCCCCAGGTCCTTGATTACCAAt                                                   >  1:193777/1‑100 (MQ=255)                                                tctCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                  <  1:1199671/98‑1 (MQ=255)                                                gTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                  <  1:964874/100‑1 (MQ=255)                                                 ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                 <  1:1004442/100‑1 (MQ=255)                                                 ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                 <  2:1749257/100‑1 (MQ=255)                                                 ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                 <  2:1537813/100‑1 (MQ=255)                                                  tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTTATTACCAATTTg                                                <  1:1562577/100‑1 (MQ=255)                                                   cgattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                                               >  1:2261635/8‑100 (MQ=255)                                                   ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGACTAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                                               >  1:435051/1‑100 (MQ=255)                                                   ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttaga                                               >  2:1199664/1‑97 (MQ=255)                                                   ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGACTCCCCAGTTCCTTGATTACCAATTTg                                                >  2:2120116/1‑99 (MQ=255)                                                    cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                              <  1:2145194/100‑1 (MQ=255)                                                    cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGACTACCAATTTGa                                               >  1:1031835/1‑99 (MQ=255)                                                     gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             <  1:1830357/100‑1 (MQ=255)                                                      aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                            <  2:1490454/100‑1 (MQ=255)                                                        tagcggtgCAATGCCACCGCTAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGCATAGTCGCCTCCACCGTTCCTTGACTACCAATTTGAAATTc                                          >  2:2146805/8‑100 (MQ=255)                                                         ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                         <  2:1238354/100‑1 (MQ=255)                                                         ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                         >  2:58579/1‑100 (MQ=255)                                                           ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                       >  2:1669019/1‑100 (MQ=255)                                                           ggTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAAGCAAC‑‑‑CC‑CCCCCTCGGGAATAGTCGCCTCCCCAGTTCCTTGATTAACAATTTGAAAATCAAt                                       >  1:635486/1‑100 (MQ=255)                                                           ggTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACGAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCGATTGGAAATTCagt                                       >  2:1182766/1‑98 (MQ=255)                                                           ggTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                         >  1:2199719/1‑98 (MQ=255)                                                            gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                      <  1:454537/100‑1 (MQ=255)                                                              ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa                                    >  2:2140344/1‑100 (MQ=255)                                                                cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGGATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTTAAATTCAACTGAgg                                  >  2:321562/1‑100 (MQ=255)                                                                cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                  >  1:620418/1‑100 (MQ=255)                                                                  ttGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa                                <  2:2197419/100‑1 (MQ=255)                                                                      cACCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                            >  1:1983442/1‑100 (MQ=255)                                                                       aCCGATAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                           <  1:2228844/100‑1 (MQ=255)                                                                            tAATCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                      <  2:7626/100‑1 (MQ=255)                                                                             aaTCGTATTTTCAGCCGTAACCAA‑‑‑‑CCTCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATGATTAttttt                     <  2:540638/100‑1 (MQ=255)                                                                              aTCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    <  1:1669032/100‑1 (MQ=255)                                                                              aTCGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    >  2:891853/1‑100 (MQ=255)                                                                               tCGTATTTTCAGGCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                     <  2:1975236/98‑1 (MQ=255)                                                                                cGTATTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                  <  1:1704290/100‑1 (MQ=255)                                                                                   aTTTTCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg               >  2:994258/1‑100 (MQ=255)                                                                                    ttttCAGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              <  2:587055/100‑1 (MQ=255)                                                                                         aGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc         >  1:2093076/1‑100 (MQ=255)                                                                                         aGCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc         >  2:577351/1‑100 (MQ=255)                                                                                          gCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc         >  1:166147/1‑99 (MQ=255)                                                                                          gCCGTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg        >  2:99098/1‑100 (MQ=255)                                                                                           ccGTAACCAA‑‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcc      >  1:2056198/1‑100 (MQ=255)                                                                                             gTAACCAAC‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAATTTCCTGATTACCAATTTGAAAATCAATTTACGAATTATTTTTTTTTATTAACCACTcccc      >  2:166438/1‑99 (MQ=255)                                                                                                accaac‑‑‑cC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  >  1:1966246/1‑100 (MQ=255)                                                                                                     |                                                                                                   AGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597207‑1597400

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