Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,075,612 C→T E1091K (GAA→AAA)  putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,075,6120CT95.1% 257.5 / 2.6 82E1091K (GAA→AAA) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (38/40);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

GCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATG  >  NC_000913/1075515‑1075706
                                                                                                 |                                                                                              
gcgAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTtgc                                                                                              <  1:1485688/100‑1 (MQ=255)
gcgAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTtgc                                                                                              <  2:1504572/100‑1 (MQ=255)
    gCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTgccgcc                                                                                          <  2:323641/100‑1 (MQ=255)
    gCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTgccgcc                                                                                          <  2:324189/100‑1 (MQ=255)
     cGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGccgccg                                                                                         <  1:1057646/100‑1 (MQ=255)
     cGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGccgccg                                                                                         >  1:1111592/1‑100 (MQ=255)
        cACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCGGACGGTAGAGAGAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGCACCACCCTTTGccgccg                                                                                         >  2:1775174/1‑97 (MQ=255)
         aCTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACg                                                                                     <  1:1319743/100‑1 (MQ=255)
         aCTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACg                                                                                     <  1:292255/100‑1 (MQ=255)
           tGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGc                                                                                   <  2:1665980/100‑1 (MQ=255)
           tGCCAGCGCACTTTCCGGGCGACTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGc                                                                                   >  2:1186063/1‑100 (MQ=255)
              cAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGc                                                                                >  2:1489630/1‑100 (MQ=255)
               aGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGca                                                                               >  1:1310360/1‑100 (MQ=255)
                 cgcACTTTCCGGGCGATTCGCCAGCAGACGGTCGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGcaca                                                                             <  2:1435946/100‑1 (MQ=255)
                 cgcACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGcaca                                                                             <  2:1017273/100‑1 (MQ=255)
                   cACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACacc                                                                           >  1:523629/1‑100 (MQ=255)
                       ttCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACaccaacc                                                                       >  1:509487/1‑100 (MQ=255)
                            ggCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTgc                                                                  >  2:1535022/1‑100 (MQ=255)
                                   cGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCccacc                                                           >  2:1113398/1‑100 (MQ=255)
                                   cGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCccacc                                                           <  2:1997182/100‑1 (MQ=255)
                                    gCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCccacca                                                          >  2:734195/1‑100 (MQ=255)
                                    gCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCccacca                                                          >  2:504528/1‑100 (MQ=255)
                                      cagcagACGGTAGAGATAGAGCGGGCCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCccac                                                            <  2:579481/96‑1 (MQ=255)
                                       agcagACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCatat                                                       >  1:703936/1‑100 (MQ=255)
                                       agcagACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCatat                                                       <  2:391874/100‑1 (MQ=255)
                                       agcagACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCatat                                                       <  2:523628/100‑1 (MQ=255)
                                           gACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACg                                                   <  1:867539/100‑1 (MQ=255)
                                            aCGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACgg                                                  <  1:1460480/100‑1 (MQ=255)
                                             cGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACgg                                                  >  2:1661900/1‑99 (MQ=255)
                                              ggTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGtt                                                >  1:1364797/1‑100 (MQ=255)
                                                tAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTaa                                              >  2:1372275/1‑100 (MQ=255)
                                                tAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTaa                                              <  2:134502/100‑1 (MQ=255)
                                                tAGAGATAGAGCGGACCGCCCGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTaa                                              <  2:266353/100‑1 (MQ=255)
                                                 agagATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAAc                                             <  1:1046105/100‑1 (MQ=255)
                                                 agagATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAAc                                             <  2:1122285/100‑1 (MQ=255)
                                                 agagATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAAc                                             <  2:229268/100‑1 (MQ=255)
                                                  gagaTAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACa                                            <  2:584139/100‑1 (MQ=255)
                                                   agaTAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACAt                                           <  2:192291/100‑1 (MQ=255)
                                                   agaTAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCGTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACAt                                           <  2:599630/100‑1 (MQ=255)
                                                    gaTAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATa                                          >  2:179571/1‑100 (MQ=255)
                                                     aTAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATAc                                         <  1:734198/100‑1 (MQ=255)
                                                     aTAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATAc                                         >  2:446967/1‑100 (MQ=255)
                                                      tAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACa                                        >  2:12573/1‑100 (MQ=255)
                                                        gagCGGACCGCCTGCTTTCGGCCCGGTACCGGACAGCCCTTTGCCGCCGAACGGCTGCACACCAACCACTGGGCCCACCATATTACGGTTAACATACAgg                                      <  2:103047/100‑1 (MQ=255)
                                                        gagCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAgg                                      <  1:1922666/100‑1 (MQ=255)
                                                        gagCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAgg                                      <  1:727977/100‑1 (MQ=255)
                                                        gagCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAgg                                      <  1:1828623/100‑1 (MQ=255)
                                                         agCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGt                                     >  2:719161/1‑100 (MQ=255)
                                                           cGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTa                                   <  1:719163/100‑1 (MQ=255)
                                                            ggACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTAc                                  >  2:914794/1‑100 (MQ=255)
                                                            ggACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTAc                                  >  1:1880767/1‑100 (MQ=255)
                                                               ccGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGCACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCaa                               >  2:264245/1‑100 (MQ=255)
                                                                 gCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCa                                >  1:240380/1‑97 (MQ=255)
                                                                 gCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACa                             >  2:1474419/1‑100 (MQ=255)
                                                                 gCCTGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACa                             >  1:438664/1‑100 (MQ=255)
                                                                 gCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACa                             >  1:1423190/1‑100 (MQ=255)
                                                                    tGCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATgg                          >  2:96419/1‑100 (MQ=255)
                                                                     gCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATggg                         >  1:482690/1‑100 (MQ=255)
                                                                     gCTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATggg                         <  1:240344/100‑1 (MQ=255)
                                                                      cTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGc                        >  1:1889502/1‑100 (MQ=255)
                                                                      cTTTCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGc                        >  1:1226781/1‑100 (MQ=255)
                                                                       tttCGGCCCGGTATCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGcc                       <  1:1190373/100‑1 (MQ=255)
                                                                       tttCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGcc                       <  2:1761852/100‑1 (MQ=255)
                                                                       tttCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGcc                       <  1:1372284/100‑1 (MQ=255)
                                                                        ttCGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGcc                       <  1:1113403/99‑1 (MQ=255)
                                                                          cGGCCCGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAg                    <  2:393248/100‑1 (MQ=255)
                                                                             cccGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCa                 <  2:644385/100‑1 (MQ=255)
                                                                              ccGGTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAg                <  1:522671/100‑1 (MQ=255)
                                                                                ggTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTg              <  2:703934/100‑1 (MQ=255)
                                                                                ggTACCGGACAACCCCTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGtt              <  2:1972254/100‑2 (MQ=255)
                                                                                 gTACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGa             >  1:1515762/1‑100 (MQ=255)
                                                                                  tACCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGAc            >  1:649124/1‑100 (MQ=255)
                                                                                   aCCGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGAcc           >  1:760345/1‑100 (MQ=255)
                                                                                    ccGGACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCt          <  1:914798/100‑1 (MQ=255)
                                                                                      ggACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGTGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTgg        >  1:1580314/1‑100 (MQ=255)
                                                                                      ggACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTgg        >  1:1986261/1‑100 (MQ=255)
                                                                                      ggACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTgg        >  2:381016/1‑100 (MQ=255)
                                                                                      ggACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTgg        <  1:1186069/100‑1 (MQ=255)
                                                                                      ggACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTgg        >  2:1315774/1‑100 (MQ=255)
                                                                                       gACAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTggg       <  1:1139580/100‑1 (MQ=255)
                                                                                         cAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCg     >  2:1813541/1‑100 (MQ=255)
                                                                                         cAACCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCg     >  1:864140/1‑100 (MQ=255)
                                                                                           aCCCTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCgct   >  1:233464/1‑98 (MQ=255)
                                                                                            cccTTTGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATg  >  1:746868/1‑100 (MQ=255)
                                                                                                 |                                                                                              
GCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATG  >  NC_000913/1075515‑1075706

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