Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,615,883 A→C intergenic (‑2/‑211) yfgO ← / → bepA putative UPF0118 family inner membrane permease/OM protein maintenance and assembly metalloprotease and chaperone, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,615,8830AC98.2% 210.9 / NA 57intergenic (‑2/‑211)yfgO/bepAputative UPF0118 family inner membrane permease/OM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (30/26);  total (31/26)
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

CCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATAAAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACTCAAAGGCGG  >  NC_000913/2615737‑2616027
                                                                                                                                                  |                                                                                                                                                 
ccAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGAc                                                                                                                                               <  2:642736/151‑1 (MQ=255)
 cAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGAcc                                                                                                                                              >  2:562529/1‑151 (MQ=255)
        aTAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑Tg                                                                                                                                      <  2:497876/151‑1 (MQ=255)
             caccaGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCt                                                                                                                                 >  2:399331/1‑151 (MQ=255)
                              gcaGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGa                                                                                                                <  2:149028/151‑1 (MQ=255)
                                aGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATc                                                                                                              >  1:66852/1‑151 (MQ=255)
                                aGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGCCGCCAAAACATTTGC‑TCAAAATTTTCGGCCTAAAGACCCTCCCTTAAACCCCAAAGGCGGGATTGcaacc                                                                                                              >  1:758542/1‑147 (MQ=255)
                                 ggCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATcc                                                                                                             >  1:743453/1‑151 (MQ=255)
                                     aCTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCa                                                                                                          >  2:555028/1‑150 (MQ=255)
                                      cTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAAt                                                                                                        >  2:890322/1‑151 (MQ=255)
                                          aGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTa                                                                                                    <  2:132852/151‑1 (MQ=255)
                                            aaaaaGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATc                                                                                                  <  1:334409/151‑1 (MQ=255)
                                             aaaaGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCg                                                                                                 <  2:347752/151‑1 (MQ=255)
                                             aaaaGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCg                                                                                                 >  1:30618/1‑151 (MQ=255)
                                                aGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGaaa                                                                                              <  2:556695/151‑1 (MQ=255)
                                                aGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGaaa                                                                                              <  1:738334/151‑1 (MQ=255)
                                                 gATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAAt                                                                                             <  1:541429/151‑1 (MQ=255)
                                                  aTAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATg                                                                                            <  1:279550/151‑1 (MQ=255)
                                                     aTGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTca                                                                                         <  1:374955/151‑1 (MQ=255)
                                                      tGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTcac                                                                                        >  1:540182/1‑151 (MQ=255)
                                                       gCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTcac                                                                                        >  1:873775/1‑150 (MQ=255)
                                                        ccAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTcaca                                                                                       >  1:807468/1‑150 (MQ=255)
                                                        ccAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTcaca                                                                                       >  1:391956/1‑150 (MQ=255)
                                                         cAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACaaa                                                                                     <  1:813211/151‑1 (MQ=255)
                                                           aaTCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACaaaaa                                                                                   <  2:107965/151‑1 (MQ=255)
                                                           aaTCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACaaaa                                                                                    >  2:734232/1‑150 (MQ=255)
                                                              ccGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAg                                                                                 <  1:855501/150‑1 (MQ=255)
                                                               cGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGAc                                                                               >  2:469443/1‑151 (MQ=255)
                                                                    aCTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGc                                                                          >  2:884636/1‑151 (MQ=255)
                                                                     cTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCt                                                                         >  1:712968/1‑151 (MQ=255)
                                                                               cAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACg                                                               <  1:800791/151‑1 (MQ=255)
                                                                                aGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACgg                                                              >  2:558888/1‑151 (MQ=255)
                                                                                aGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACgg                                                              >  1:320831/1‑151 (MQ=255)
                                                                                  cAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGa                                                             <  2:585907/150‑1 (MQ=255)
                                                                                        aaTCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAgg                                                       <  2:685890/150‑1 (MQ=255)
                                                                                           cGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAAc                                                   >  2:203005/1‑151 (MQ=255)
                                                                                                 cGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGaaa                                              >  2:461594/1‑150 (MQ=255)
                                                                                                  gggTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGaaa                                              >  1:711120/1‑149 (MQ=255)
                                                                                                   ggTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGt                                           <  2:163189/151‑1 (MQ=255)
                                                                                                         tAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCAttt                                      >  2:796664/1‑150 (MQ=255)
                                                                                                                       tACCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCtt                       >  2:664529/1‑151 (MQ=255)
                                                                                                                        aCCATTGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCtt                       <  1:362366/150‑1 (MQ=255)
                                                                                                                             tGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCg                  <  1:399317/150‑1 (MQ=255)
                                                                                                                             tGCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCg                  <  1:890301/150‑1 (MQ=255)
                                                                                                                              gCATCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGtt                <  1:283952/151‑1 (MQ=255)
                                                                                                                                aTCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTac              >  2:178479/1‑151 (MQ=255)
                                                                                                                                aTCAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTac              >  2:158596/1‑151 (MQ=255)
                                                                                                                                  cAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTacac            >  2:107264/1‑151 (MQ=255)
                                                                                                                                  cAACATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTacac            >  1:448523/1‑151 (MQ=255)
                                                                                                                                   aaCATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACt           <  2:608441/151‑1 (MQ=255)
                                                                                                                                     cATTTCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACt           >  1:523309/1‑149 (MQ=255)
                                                                                                                                       tttCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACTCaa        >  2:386920/1‑150 (MQ=255)
                                                                                                                                        ttCGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACTCaaa       <  1:790927/150‑1 (MQ=255)
                                                                                                                                          cGAGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACTCAAAggc    >  2:759694/1‑151 (MQ=255)
                                                                                                                                            aGCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACTCAAAggcgg  <  2:319625/151‑1 (MQ=255)
                                                                                                                                             gCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACTCAAAggcgg  >  2:342456/1‑150 (MQ=255)
                                                                                                                                             gCATACAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACGTCGCTTTTTATGACGGATGCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACTCAAAggcgg  <  2:835792/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCATAAAGACCTTCCC‑TGAACCTCAAGAGCGGGATTGCGATCCGCAATTGTATCGAAATGTCACAAAAAAGACTTCGCTTTTTATGACGGATTCAGGAAACTGAAAAGTCATTTGAGTGGGCTAATCTTCGCCGTTACACTCAAAGGCGG  >  NC_000913/2615737‑2616027

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