Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,358,163 G→A A228V (GCG→GTG)  rhmA ← 2‑keto‑3‑deoxy‑L‑rhamnonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,358,1630GA96.9% 190.4 / ‑2.3 64A228V (GCG→GTG) rhmA2‑keto‑3‑deoxy‑L‑rhamnonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (29/33);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

GGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCGGCGTTATCCGGGTAGCCCAACG  >  NC_000913/2358030‑2358298
                                                                                                                                     |                                                                                                                                       
gggCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAgg                                                                                                                                    <  1:1178459/139‑1 (MQ=255)
      ttGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgc                                                                                                                              <  2:1356625/139‑1 (MQ=255)
       tGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCCAGCAAtgc                                                                                                                              <  1:732839/138‑1 (MQ=255)
       tGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgct                                                                                                                             >  2:1713503/1‑139 (MQ=255)
         atatCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgctgc                                                                                                                           <  2:332136/139‑1 (MQ=255)
                tAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCATTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCatat                                                                                                                    >  2:1271036/1‑139 (MQ=255)
                   cTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTgcg                                                                                                                          >  2:71697/1‑130 (MQ=255)
                   cTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTgcg                                                                                                                          <  1:71697/130‑1 (MQ=255)
                   cTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAg                                                                                                                 <  1:498170/139‑1 (MQ=255)
                     aCCTTTTATGCGTAGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGcc                                                                                                                        >  2:468588/1‑130 (MQ=255)
                     aCCTTTTATGCGTAGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGcc                                                                                                                        <  1:468588/130‑1 (MQ=255)
                       cTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGc                                                                                                             <  1:1271036/139‑1 (MQ=255)
                          ttATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAACGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCtac                                                                                                          >  1:1243061/1‑139 (MQ=255)
                           tATGCGTGGCACATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCcaca                                                                                                         <  2:1409949/139‑1 (MQ=255)
                             tGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGc                                                                                                       >  1:952841/1‑139 (MQ=255)
                                gTGGCCCATTTTTGCCTGATTTAAACATCGCCTGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAg                                                                                                    >  1:1116216/1‑139 (MQ=255)
                                gTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAg                                                                                                    >  2:1509094/1‑139 (MQ=255)
                                 tGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCAtt                                                                                                                                >  1:1673587/1‑110 (MQ=255)
                                 tGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCAtt                                                                                                                                <  2:1673587/110‑1 (MQ=255)
                                 tGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAg                                                                                                                 >  2:1305088/1‑125 (MQ=255)
                                 tGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAg                                                                                                                 <  1:1305088/125‑1 (MQ=255)
                                 tGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGa                                                                                                   >  1:1847855/1‑139 (MQ=255)
                                  ggCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGaa                                                                                                  >  1:907519/1‑139 (MQ=255)
                                        tttttGCCCGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGATAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAg                                                                                            <  2:1116216/139‑1 (MQ=255)
                                         ttttGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGc                                                                                           >  2:161966/1‑139 (MQ=255)
                                               cTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCttt                                                                                     <  1:1713503/139‑1 (MQ=255)
                                                  aTTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTAcc                                                                                  >  2:220134/1‑139 (MQ=255)
                                                  aTTTAAACATCCCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTAcc                                                                                  <  2:1223696/139‑1 (MQ=255)
                                                    ttAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCg                                                                                <  2:413659/139‑1 (MQ=255)
                                                    ttAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCg                                                                                <  1:1429378/139‑1 (MQ=255)
                                                            cGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTAcc                                                                                  >  1:461589/1‑129 (MQ=255)
                                                            cGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTAcc                                                                                  <  2:461589/129‑1 (MQ=255)
                                                            cGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTgcg                                                                                                                          >  2:1022182/1‑89 (MQ=255)
                                                            cGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTgcg                                                                                                                          <  1:1022182/89‑1 (MQ=255)
                                                            cGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACgg                                                                        >  1:1579414/1‑139 (MQ=255)
                                                              ccAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGAt                                                                      <  2:1579414/139‑1 (MQ=255)
                                                                 gTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCGGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATccg                                                                   >  2:1036324/1‑139 (MQ=255)
                                                                   cGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATccgcc                                                                 <  2:907519/139‑1 (MQ=255)
                                                                        aTCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCGGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATa                                                            <  1:1036324/139‑1 (MQ=255)
                                                                        aTCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATa                                                            <  2:952841/139‑1 (MQ=255)
                                                                            aGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTgg                                                        >  1:892286/1‑139 (MQ=255)
                                                                             gggCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATa                                                            <  1:1115529/134‑1 (MQ=255)
                                                                             gggCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATa                                                            >  2:1115529/1‑134 (MQ=255)
                                                                              ggCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGTCAGAAAACCAGCCGCTTTACCCg                                                                                >  2:435892/1‑113 (MQ=255)
                                                                              ggCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGTCAGAAAACCAGCCGCTTTACCCg                                                                                <  1:435892/113‑1 (MQ=255)
                                                                              ggCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGtt                                                      >  1:643130/1‑139 (MQ=255)
                                                                                cATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGt                                                       <  1:1492997/136‑1 (MQ=255)
                                                                                cATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGt                                                       >  2:1492997/1‑136 (MQ=255)
                                                                                                cGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCGAGAAAACCAGCCGCTTTACCCGCAGCACGGAt                                                                      <  1:1118750/105‑1 (MQ=255)
                                                                                                cGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGCAGCCACAGCCAGAAAACCAGCCGCTTTAACCGCAGCACGGAt                                                                      >  2:1118750/1‑105 (MQ=255)
                                                                                                  tgtCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCgcagc                                                                            >  2:1805310/1‑97 (MQ=255)
                                                                                                  tgtCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCgcagc                                                                            <  1:1805310/97‑1 (MQ=255)
                                                                                                  tgtCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTc                                  >  2:1120136/1‑139 (MQ=255)
                                                                                                          ccAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGcc                          <  2:1847855/139‑1 (MQ=255)
                                                                                                             aCAGCGACAAAGTTCGCTCCCCGCACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCAAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCgg                       <  2:1452380/139‑1 (MQ=255)
                                                                                                             aCAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCcaca                                                                                                         <  1:1178536/57‑1 (MQ=255)
                                                                                                             aCAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCcaca                                                                                                         >  2:1178536/1‑57 (MQ=255)
                                                                                                                 cGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGc                           <  1:626731/131‑1 (MQ=255)
                                                                                                                 cGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGc                           >  2:626731/1‑131 (MQ=255)
                                                                                                                   aCAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCg                                           <  1:51267/113‑1 (MQ=255)
                                                                                                                   aCAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCg                                           >  2:51267/1‑113 (MQ=255)
                                                                                                                   aCAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCGGCGTTAt                 <  1:1509094/139‑1 (MQ=255)
                                                                                                                        gTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATAGTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCGGCGTTATCCggg            <  1:161966/139‑1 (MQ=255)
                                                                                                                         ttCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCGGCGTTATCCGGGt           >  2:1195668/1‑139 (MQ=255)
                                                                                                                                  cacaCCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCGGCGTTATCCGGGTAGCCCAacg  <  2:120531/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
GGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCGGCGTTATCCGGGTAGCCCAACG  >  NC_000913/2358030‑2358298

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