Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,907,539 T→C Y18C (TAC→TGC)  yobF ← DUF2527 family heat‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,907,5390TC100.0% 278.3 / NA 76Y18C (TAC→TGC) yobFDUF2527 family heat‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/36);  total (40/36)

TTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATG  >  NC_000913/1907410‑1907669
                                                                                                                                 |                                                                                                                                  
ttAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATtcaacgt                                                                                                                           >  2:302249/1‑139 (MQ=255)
   aCCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATtcaacgtcaa                                                                                                                        <  1:599853/139‑1 (MQ=255)
         ccTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGtt                                                                                                                  >  1:85275/1‑139 (MQ=255)
               aTCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCa                                                                                                            >  1:1174809/1‑139 (MQ=255)
                    tGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGAct                                                                                                       >  1:237857/1‑139 (MQ=255)
                    tGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGAct                                                                                                       >  1:1508539/1‑139 (MQ=255)
                              aaaaTTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGa                                                                                             <  1:302249/139‑1 (MQ=255)
                                                   ttcGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGaa                                                                        <  1:1047829/139‑1 (MQ=255)
                                                        ccGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTa                                                                   >  1:188595/1‑139 (MQ=255)
                                                                      tAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCg                                                     <  2:1174809/139‑1 (MQ=255)
                                                                        gATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGc                                                                                       >  1:65949/1‑103 (MQ=255)
                                                                        gATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGc                                                                                       <  2:65949/103‑1 (MQ=255)
                                                                        gATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCcacac                                                                                  >  2:1276711/1‑108 (MQ=255)
                                                                        gATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCcacac                                                                                  <  1:1276711/108‑1 (MQ=255)
                                                                        gATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATaaa                                                                             >  1:613499/1‑113 (MQ=255)
                                                                        gATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATaaa                                                                             <  2:613499/113‑1 (MQ=255)
                                                                         aTAAAACTGACACGTTACTGCATAAGGCACCAATATAAGGCTCGGCAGAGAAGCTGCATTCAACGACAACGTGTTTACTTAGGACTTCTTTACTGAAAACGCCACACATAAACAGAACTGTACCTCGTTTAACCCGaaa                                                  <  2:994722/139‑1 (MQ=255)
                                                                         aTAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGaa                                                                        <  2:142217/117‑1 (MQ=255)
                                                                         aTAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGaa                                                                        >  1:142217/1‑117 (MQ=255)
                                                                             aaCTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTg                                              <  2:188595/139‑1 (MQ=255)
                                                                              aCTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTg                                              >  1:938063/1‑138 (MQ=255)
                                                                              aCTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTg                                              <  2:938063/138‑1 (MQ=255)
                                                                                  acacATTACTGCATGAGGCACCAATATAAGGCTCGGGAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCAACACATAAACAGAACTGTACGTCGTTTAACGCGAAATCTGTTATc                                         >  2:1409264/1‑139 (MQ=255)
                                                                                  acacATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATc                                         >  1:252874/1‑139 (MQ=255)
                                                                                     caTTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATcaca                                      <  2:85275/139‑1 (MQ=255)
                                                                                        tACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATAc                                   >  2:1458770/1‑139 (MQ=255)
                                                                                         aCTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGttt                                                           >  1:272484/1‑114 (MQ=255)
                                                                                         aCTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGttt                                                           <  2:272484/114‑1 (MQ=255)
                                                                                            gCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACg                               <  1:1336264/139‑1 (MQ=255)
                                                                                             cATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACt                                                                                                              >  1:1411979/1‑59 (MQ=255)
                                                                                             cATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACt                                                                                                              <  2:1411979/59‑1 (MQ=255)
                                                                                               tGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGc                                                                                       <  2:325175/80‑1 (MQ=255)
                                                                                               tGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGc                                                                                       >  1:325175/1‑80 (MQ=255)
                                                                                                gAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATtcaacgtcaa                                                                                                                        <  1:196830/46‑1 (MQ=255)
                                                                                                gAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATtcaacgtcaa                                                                                                                        >  2:196830/1‑46 (MQ=255)
                                                                                                gAGGCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTaa                           >  1:238684/1‑139 (MQ=255)
                                                                                                  ggCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGaaaa                                                                                          <  2:1425496/74‑1 (MQ=255)
                                                                                                  ggCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGaaaa                                                                                          >  1:1425496/1‑74 (MQ=255)
                                                                                                  ggCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGaa                                                                                            <  2:1403630/72‑1 (MQ=255)
                                                                                                  ggCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGaa                                                                                            >  1:1403630/1‑72 (MQ=255)
                                                                                                  ggCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACAt                                                                                >  1:1288168/1‑84 (MQ=255)
                                                                                                  ggCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACAt                                                                                <  2:1288168/84‑1 (MQ=255)
                                                                                                  ggCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATa                                    >  2:460151/1‑128 (MQ=255)
                                                                                                  ggCACCAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATa                                    <  1:460151/128‑1 (MQ=255)
                                                                                                       cAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCcacac                                                                                  >  2:687949/1‑77 (MQ=255)
                                                                                                       cAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCcacac                                                                                  <  1:687949/77‑1 (MQ=255)
                                                                                                       cAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATaaa                                                                             >  1:564124/1‑82 (MQ=255)
                                                                                                       cAATATAAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATaaa                                                                             <  2:564124/82‑1 (MQ=255)
                                                                                                           ataAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGc                <  1:1458770/139‑1 (MQ=255)
                                                                                                           ataAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGc                <  2:1508539/139‑1 (MQ=255)
                                                                                                            taAGGCTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGcc               >  2:278645/1‑139 (MQ=255)
                                                                                                               ggCTCGGCAGAGAAGCGGCATTCAAAGCCAACGTGTTTACCCAGGACTTCTTTACTGAAAATGCCCCACATAAAAAGAACTGTACTTCGATTAACCCGAAAACTGTTATAACATACAACATTAATTATGGCAGGCAAtt            >  2:741563/1‑139 (MQ=255)
                                                                                                                 cTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCcacac                                                                                  <  2:1184032/67‑1 (MQ=255)
                                                                                                                 cTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCcacac                                                                                  >  1:1184032/1‑67 (MQ=255)
                                                                                                                 cTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATcaca                                      <  2:1149865/111‑1 (MQ=255)
                                                                                                                 cTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATcaca                                      >  1:1149865/1‑111 (MQ=255)
                                                                                                                 cTCGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCAtttt          >  2:559903/1‑139 (MQ=255)
                                                                                                                   cGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGtt                                            >  1:1102473/1‑103 (MQ=255)
                                                                                                                   cGGCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGtt                                            <  2:1102473/103‑1 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGaa                                                                        <  1:133349/74‑1 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGaa                                                                        >  2:133349/1‑74 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGt                                                             >  2:1503965/1‑85 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGt                                                             <  1:1503965/85‑1 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTg                                              >  1:556387/1‑100 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTg                                              <  2:556387/100‑1 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTAt                                          >  1:466238/1‑104 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTAt                                          <  2:466238/104‑1 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACg                               >  2:1404527/1‑115 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACg                               <  1:1404527/115‑1 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGcc               <  2:1347059/131‑1 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGcc               >  1:1347059/1‑131 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTaa       >  1:216063/1‑139 (MQ=255)
                                                                                                                    ggCAGAGAAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTaa       >  1:1446622/1‑139 (MQ=255)
                                                                                                                         agaAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATa                                    <  2:1064046/105‑1 (MQ=255)
                                                                                                                         agaAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATa                                    >  1:1064046/1‑105 (MQ=255)
                                                                                                                         agaAGCGGCATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATg  <  2:216063/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                  
TTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATG  >  NC_000913/1907410‑1907669

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