Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,927,062 +G coding (2/1284 nt)
coding (169/183 nt)
yieM ←
b4564 ←
hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,927,0621.G99.4% 420.1 / ‑3.7 154coding (2/1284 nt)
coding (169/183 nt)
yieM
b4564
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (72/81);  total (72/82)
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.48e-01

Only 100 of 166 total aligned reads displayed.

AGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCACGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAC  >  NC_000913/3926987‑3927131
                                                                            |                                                                     
agCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCa                                                                        <  1:1291687/76‑1 (MQ=255)
 gCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCag                                                                       >  2:2651374/1‑75 (MQ=255)
 gCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCag                                                                       <  1:167361/76‑2 (MQ=255)
  cTGGCTTTGCTGATAGCACATCGCCTCTACGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagc                                                                      <  1:1328861/76‑3 (MQ=255)
  cTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTTCGCagc                                                                      <  1:5121885/76‑3 (MQ=255)
  cTGGCTTTGCTGATAGCACATCACCTCGTCGGCGAGTTCCGGCGGGACTCGGGCAGCTTTCAGCCGACTGCGCagc                                                                      >  1:5612602/1‑74 (MQ=255)
   tGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGACATCTTTCAGCCGACTGCGCagcg                                                                     >  1:3324093/1‑73 (MQ=255)
      cTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCtt                                                                  <  1:3324594/76‑1 (MQ=255)
       tttGCTGATAGCTCATCACCTCTTCGGTGCGTTCCGTCGGGACTCGGGCATCTTTCAGCCGACTGCTCAGCGCTTc                                                                 <  1:1494766/76‑1 (MQ=255)
       tttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTc                                                                 <  1:481189/76‑1 (MQ=255)
       tttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTc                                                                 <  1:5593827/76‑1 (MQ=255)
       tttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTc                                                                 <  2:2656845/76‑1 (MQ=255)
        ttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCa                                                                >  1:4205736/1‑76 (MQ=255)
         tGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCAc                                                               <  1:5066046/76‑1 (MQ=255)
           cTGATAGCACATCCCCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGc                                                             <  1:239382/76‑1 (MQ=255)
            tGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGcc                                                            <  1:734773/76‑1 (MQ=255)
             gATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCa                                                           >  1:1917621/1‑76 (MQ=255)
               tAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGc                                                         >  1:4329585/1‑76 (MQ=255)
               tAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGc                                                         <  1:4498127/76‑1 (MQ=255)
               tAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGc                                                         <  2:5311436/76‑1 (MQ=255)
                 gCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgg                                                       <  1:4402130/76‑1 (MQ=255)
                   acaTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgggg                                                     >  1:1076499/1‑76 (MQ=255)
                     aTCCCCGCTTCGGTGAGTTCGGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGaa                                                   <  2:1917621/76‑1 (MQ=255)
                     aTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCCGCggggga                                                   >  1:322506/1‑74 (MQ=255)
                     aTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCCGCGGGGaa                                                   >  2:1008540/1‑76 (MQ=255)
                     aTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGaa                                                   <  1:2060202/76‑1 (MQ=255)
                     aTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGaa                                                   >  2:2189818/1‑76 (MQ=255)
                     aTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGAGCAGCGCTTCACGCCAGCGGGGaa                                                   >  1:2991178/1‑76 (MQ=255)
                      tctccTCTTCGGTGCGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCTCTTCACGCCAGCGGGGAAc                                                  <  1:738122/73‑1 (MQ=255)
                       cACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAca                                                 >  2:5306386/1‑76 (MQ=255)
                       cACCGCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAca                                                 <  2:73024/76‑1 (MQ=255)
                        cccTCTTTGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcat                                                <  2:3955646/75‑1 (MQ=255)
                        cccTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcat                                                <  2:4369299/75‑1 (MQ=255)
                        acgtctcCGGTGAGTGCCGGCGGGACTCGTGCATCTTTCAGCCGACTGCGCAGCGCGTCACGCCAGCGGGGAAcat                                                <  2:4406770/73‑1 (MQ=255)
                        aCCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcat                                                <  2:5612600/76‑1 (MQ=255)
                        aCCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcat                                                >  2:4540461/1‑76 (MQ=255)
                          ctctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCAAGCGGGGAAcatca                                              <  1:5195418/76‑1 (MQ=255)
                            ctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACatcatc                                            >  2:817930/1‑76 (MQ=255)
                            ctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACatcatc                                            <  1:1431519/76‑1 (MQ=255)
                             ttCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACAtcatca                                           >  2:4230540/1‑76 (MQ=255)
                              tCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAg                                          >  1:4953098/1‑76 (MQ=255)
                               cGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGt                                         >  1:5387208/1‑76 (MQ=255)
                                ggTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTg                                        >  1:4894078/1‑76 (MQ=255)
                                 gTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGa                                       >  1:1975129/1‑76 (MQ=255)
                                  tGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGAt                                      <  2:3060021/76‑1 (MQ=255)
                                   gAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGAtt                                     >  1:2906225/1‑76 (MQ=255)
                                     gTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgc                                   <  2:1026986/76‑1 (MQ=255)
                                       tCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgcca                                 <  2:2805454/76‑2 (MQ=255)
                                        ccGGGGGGACTCGGGGATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccgc                                <  2:5460878/76‑1 (MQ=255)
                                        ccGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccgc                                >  1:2782527/1‑76 (MQ=255)
                                         cggcggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccgcc                               <  1:201262/76‑1 (MQ=255)
                                         cggcggGACTCGCTCATCTTTCAGCCGGCTGCGCAGCGCTTCACGCGAGCGGGGAACATCATCAGTGATTgccgcc                               <  1:2824220/76‑1 (MQ=255)
                                          ggcggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCt                              <  1:1734700/76‑1 (MQ=255)
                                            cggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTc                            <  2:4323039/76‑1 (MQ=255)
                                            cggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTc                            >  1:3826581/1‑76 (MQ=255)
                                             gggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCa                           <  1:41452/76‑1 (MQ=255)
                                             gggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCa                           <  2:2906225/76‑1 (MQ=255)
                                             gggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCa                           >  1:1130805/1‑76 (MQ=255)
                                              ggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAg                          <  2:4630509/76‑1 (MQ=255)
                                                aCTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGc                        <  2:3693505/76‑1 (MQ=255)
                                                 cTCGGTCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCg                       <  2:4790411/76‑1 (MQ=255)
                                                  tCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGt                      <  2:3778477/76‑1 (MQ=255)
                                                   cGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTg                     >  1:1052710/1‑76 (MQ=255)
                                                   cGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTg                     >  1:5101329/1‑76 (MQ=255)
                                                    gggCATCTTTCAGCCGACGTCGCAGCGCTTCCCGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTgg                    <  2:291632/76‑1 (MQ=255)
                                                    gggCAGCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTgg                    <  2:5653076/76‑1 (MQ=255)
                                                     ggCATCTTTCAGCCGACTGCGCAGAGATTCACGCAAGCGGGGACCATCATCAGTGATTGCCGCCTTCCGGCGttgg                   <  1:4812927/76‑4 (MQ=255)
                                                      gCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGa                  <  1:1206598/76‑1 (MQ=255)
                                                      gCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGa                  >  2:4307212/1‑76 (MQ=255)
                                                      gCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGa                  <  2:495401/76‑1 (MQ=255)
                                                       cATCTTTCCGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGACTTCCGGCGTGGGaa                 >  1:1840060/1‑76 (MQ=255)
                                                       cATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGaa                 >  1:1836427/1‑76 (MQ=255)
                                                       cATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGaa                 >  1:2430810/1‑76 (MQ=255)
                                                       cATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATCGCCGCCTTCAGGTGTGGGaa                 >  2:606059/1‑76 (MQ=255)
                                                       cATCTTTCAGCCGACCGCGCAGCGCTTCACGCCAGCGGGGAACAGCATCAGTGATTGCCGCCTTCAGGCGTgggca                 >  2:3191409/1‑74 (MQ=255)
                                                       cATCTCCCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACAGCAGCAGTGAGTGCCGCCTTCCGGCGTgggca                 >  2:3167320/1‑74 (MQ=255)
                                                        gtCTTCCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAAt                <  1:91303/75‑1 (MQ=255)
                                                           tttCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAAtttt             >  1:41225/1‑76 (MQ=255)
                                                           tttCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAAtttt             <  1:4952449/76‑1 (MQ=255)
                                                            ttCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAAttttt            >  2:2510732/1‑76 (MQ=255)
                                                             tCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCCGGCGTGGGAATTTTTc           >  2:5028202/1‑76 (MQ=255)
                                                             tCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTc           >  2:1368161/1‑76 (MQ=255)
                                                              cAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCa          <  2:4205735/76‑1 (MQ=255)
                                                              cAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCa          >  1:3650437/1‑76 (MQ=255)
                                                               aGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCaa         >  1:4711861/1‑76 (MQ=255)
                                                               aGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCaa         >  2:724394/1‑76 (MQ=255)
                                                               aGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCaa         <  2:2430810/76‑1 (MQ=255)
                                                                gCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCaaa        <  1:5318811/76‑1 (MQ=255)
                                                                 ccGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaag       <  2:4953097/76‑1 (MQ=255)
                                                                 ccGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaag       <  1:1120494/76‑1 (MQ=255)
                                                                  cGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaaga      >  1:2134524/1‑76 (MQ=255)
                                                                  cGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaaga      >  2:2456571/1‑76 (MQ=255)
                                                                   gACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAAtt               >  1:4705552/1‑66 (MQ=255)
                                                                   gACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaagaa     >  1:4420660/1‑76 (MQ=255)
                                                                    gcTGCGCAGCGCTTCACGCCAGGGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaagaag    <  1:2249877/75‑1 (MQ=255)
                                                                    aCTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaagaag    <  1:3710353/76‑1 (MQ=255)
                                                                      tGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAc  >  1:341696/1‑76 (MQ=255)
                                                                      tGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAc  >  1:5722745/1‑76 (MQ=255)
                                                                      tGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAc  >  2:3583043/1‑76 (MQ=255)
                                                                      tGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAc  <  2:1262330/76‑1 (MQ=255)
                                                                            |                                                                     
AGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCACGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAC  >  NC_000913/3926987‑3927131

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