Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,347,069 G→A *164* (TAG→TAA ptsN → sugar‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,347,0690GA100.0% 246.8 / NA 76*164* (TAG→TAAptsNsugar‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/39);  total (37/39)

TGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAGTTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTGCGTGGAT  >  NC_000913/3346939‑3347204
                                                                                                                                  |                                                                                                                                       
tGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCgg                                                                                                                                 <  1:198863/139‑1 (MQ=255)
      ctctGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           <  2:955951/139‑1 (MQ=255)
      ctctGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           <  2:1460542/139‑1 (MQ=255)
         tGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGtctc                                                                                                                        >  2:918078/1‑139 (MQ=255)
            cGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           >  2:963627/1‑133 (MQ=255)
            cGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           <  1:963627/133‑1 (MQ=255)
             gCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCtttt                                                                                                                    >  1:1228799/1‑139 (MQ=255)
              cTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTa                                                                                                                   >  2:1169812/1‑139 (MQ=255)
               tggtggCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAg                                                                                                                  <  1:1169812/139‑1 (MQ=255)
               tggtggCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAg                                                                                                                  <  2:1228799/139‑1 (MQ=255)
                 gtggCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGGTGAAGCGTAATTATTCGGTAATGTCTCTTTTAGAc                                                                                                                >  2:1010278/1‑139 (MQ=255)
                     cGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGt                                                                                                                                <  1:526607/119‑1 (MQ=255)
                     cGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGt                                                                                                                                >  2:526607/1‑119 (MQ=255)
                                    aCAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTc                                                                                                                                   >  1:1340096/1‑101 (MQ=255)
                                    aCAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTc                                                                                                                                   <  2:1340096/101‑1 (MQ=255)
                                      aaaaCCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                           <  1:918078/139‑1 (MQ=255)
                                                  cgccgTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           >  1:1638436/1‑95 (MQ=255)
                                                  cgccgTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           <  2:1638436/95‑1 (MQ=255)
                                                  cgccgTTTACGCACAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           <  2:942651/95‑1 (MQ=255)
                                                  cgccgTTTACGCACAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           >  1:942651/1‑95 (MQ=255)
                                                    ccgTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTAt                                                                                                                                     <  2:751383/83‑1 (MQ=255)
                                                    ccgTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTAt                                                                                                                                     >  1:751383/1‑83 (MQ=255)
                                                             cgcAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGCTGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGa                                                                    <  1:1010278/139‑1 (MQ=255)
                                                                aGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgt                                                                 >  1:1306439/1‑139 (MQ=255)
                                                                  cccAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           <  2:606373/79‑1 (MQ=255)
                                                                  cccAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                           >  1:606373/1‑79 (MQ=255)
                                                                            tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTa                                                                                                                               >  2:1376917/1‑65 (MQ=255)
                                                                            tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTa                                                                                                                               <  1:1376917/65‑1 (MQ=255)
                                                                            tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTgagg                                                                                                       >  1:1513649/1‑89 (MQ=255)
                                                                            tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTgagg                                                                                                       <  2:1513649/89‑1 (MQ=255)
                                                                            tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                           <  1:156134/101‑1 (MQ=255)
                                                                            tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                           >  2:156134/1‑101 (MQ=255)
                                                                            tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTtcaggttcaggt                                                     >  2:1513651/1‑139 (MQ=255)
                                                                            tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTtcaggttcaggt                                                     >  1:125272/1‑139 (MQ=255)
                                                                                      tATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGAc                                                                                                                >  1:718010/1‑70 (MQ=255)
                                                                                      tATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGAc                                                                                                                <  2:718010/70‑1 (MQ=255)
                                                                                        tCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                             <  1:523585/87‑1 (MQ=255)
                                                                                        tCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                             >  2:523585/1‑87 (MQ=255)
                                                                                        tCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGcc                                         >  2:135747/1‑139 (MQ=255)
                                                                                               aTTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgcgtg                                  <  1:135747/139‑1 (MQ=255)
                                                                                                ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGCAACAGTACATGGTACTGATGATCGTCAGCGGa                                                                    <  1:416944/104‑1 (MQ=255)
                                                                                                ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgcgtgc                                 >  2:710351/1‑139 (MQ=255)
                                                                                                ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATCTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGa                                                                    >  2:416944/1‑104 (MQ=255)
                                                                                                     agtACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTgg                            <  2:296195/137‑1 (MQ=255)
                                                                                                      aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                           <  1:192386/75‑1 (MQ=255)
                                                                                                      aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                           >  2:192386/1‑75 (MQ=255)
                                                                                                      aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGa                           <  2:376367/139‑1 (MQ=255)
                                                                                                       tACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGaa                          <  1:1513651/139‑1 (MQ=255)
                                                                                                        aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttc                                                               <  2:765236/101‑1 (MQ=255)
                                                                                                        aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttc                                                               >  1:765236/1‑101 (MQ=255)
                                                                                                        aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgc                                     <  1:1325487/127‑1 (MQ=255)
                                                                                                        aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgc                                     >  2:1325487/1‑127 (MQ=255)
                                                                                                          cGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGcc                                         <  1:324707/121‑1 (MQ=255)
                                                                                                          cGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGcc                                         >  2:324707/1‑121 (MQ=255)
                                                                                                           gAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGata                      <  2:1306439/139‑1 (MQ=255)
                                                                                                                  cTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGttt               >  1:508032/1‑139 (MQ=255)
                                                                                                                  cTCCGGATGAAGCGTAATTATTCGGTAATGTCTCGTTTAGACGTTGCGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCTGACGGTCAGGTTCAGGCAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGtct               >  1:296195/1‑137 (MQ=255)
                                                                                                                   tCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTAc                                                                                     <  1:445090/68‑1 (MQ=255)
                                                                                                                   tCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTAc                                                                                     >  2:445090/1‑68 (MQ=255)
                                                                                                                   tCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgcg                                    <  2:1081706/117‑1 (MQ=255)
                                                                                                                   tCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgcg                                    >  1:1081706/1‑117 (MQ=255)
                                                                                                                    ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACg                                                                                                               >  2:1254166/1‑41 (MQ=255)
                                                                                                                    ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACg                                                                                                               <  1:1254166/41‑1 (MQ=255)
                                                                                                                    ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttcaggtt                                                          >  2:1268604/1‑94 (MQ=255)
                                                                                                                    ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttcaggtt                                                          <  1:1268604/94‑1 (MQ=255)
                                                                                                                       gATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgtt                                                                >  2:1008126/1‑85 (MQ=255)
                                                                                                                       gATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgtt                                                                <  1:1008126/85‑1 (MQ=255)
                                                                                                                       gATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTaa                                                   >  1:102686/1‑98 (MQ=255)
                                                                                                                       gATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTaa                                                   <  2:102686/98‑1 (MQ=255)
                                                                                                                        aTGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTg         <  2:1308102/139‑1 (MQ=255)
                                                                                                                          gAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                             >  2:467723/1‑53 (MQ=255)
                                                                                                                          gAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                             <  1:467723/53‑1 (MQ=255)
                                                                                                                             gCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaa                                                                                                   <  1:937196/44‑1 (MQ=255)
                                                                                                                             gCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaa                                                                                                   >  2:937196/1‑44 (MQ=255)
                                                                                                                              cGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTGCGTGGa   >  2:1378355/1‑139 (MQ=255)
                                                                                                                               gTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTGCGTGGAt  <  1:747067/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
TGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAGTTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTGCGTGGAT  >  NC_000913/3346939‑3347204

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