Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,000,345 (A)5→4 coding (722/1188 nt) ygeW → knotted carbamoyltransferase YgeW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,000,3410A.100.0% 269.9 / NA 61coding (718/1188 nt)ygeWknotted carbamoyltransferase YgeW
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (32/29);  total (32/29)

CCTACTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCGTACT  >  NZ_CP009273/3000207‑3000466
                                                                                                                                      |                                                                                                                             
ccTACTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGaaaac                                                                                                                           <  2:487971/139‑2 (MQ=255)
 cTACTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGaaaaca                                                                                                                          <  2:591462/139‑3 (MQ=255)
   aCTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGaaaacaat                                                                                                                        >  2:5989/1‑135 (MQ=255)
   aCTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGaaaacaat                                                                                                                        >  2:104128/1‑135 (MQ=255)
   aCTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCACCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGaaaacaat                                                                                                                        >  2:337196/1‑135 (MQ=255)
                 ggCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCgg                                                                                                         >  2:79025/1‑139 (MQ=255)
                 ggCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCgg                                                                                                         >  1:198251/1‑139 (MQ=255)
                                  aCCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAgg                                                                                        >  2:199667/1‑139 (MQ=255)
                                       aaGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTc                                                                                                            >  1:118856/1‑114 (MQ=255)
                                       aaGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTc                                                                                                            <  2:118856/114‑1 (MQ=255)
                                         ggCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCt                                                                                                                    >  2:149573/1‑104 (MQ=255)
                                         ggCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCt                                                                                                                    <  1:149573/104‑1 (MQ=255)
                                         ggCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAg                                                                                 >  1:331294/1‑139 (MQ=255)
                                         ggCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAg                                                                                 >  1:483594/1‑139 (MQ=255)
                                           catcatCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCa                                                                               <  1:432733/139‑1 (MQ=255)
                                              catcGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATg                                                                             >  2:393893/1‑138 (MQ=255)
                                                 cGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAg                                                                                         <  2:349347/123‑1 (MQ=255)
                                                 cGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAg                                                                                         >  1:349347/1‑123 (MQ=255)
                                                        atgaCTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCtt                                                                  >  1:360426/1‑139 (MQ=255)
                                                           aCTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCa                                                                >  1:447418/1‑138 (MQ=255)
                                                            cTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTACGTCAGGTCACCAGCATGGAAGAAGCCTTCaaa                                                              <  1:312394/139‑1 (MQ=255)
                                                                   cGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCa                                                        <  2:60463/138‑1 (MQ=255)
                                                                    ggTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGa                                                      <  2:442916/139‑1 (MQ=255)
                                                                           aTGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGtt                                               >  2:110520/1‑139 (MQ=255)
                                                                           aTGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGtt                                               <  1:363910/139‑1 (MQ=255)
                                                                                aCCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGgaag                                                                         <  2:119257/108‑1 (MQ=255)
                                                                                aCCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGgaag                                                                         >  1:119257/1‑108 (MQ=255)
                                                                                aCCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATcc                                          >  2:587004/1‑139 (MQ=255)
                                                                                aCCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATcc                                          >  1:193041/1‑139 (MQ=255)
                                                                                aCCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATcc                                          >  1:466101/1‑139 (MQ=255)
                                                                                    tGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAg                                      >  1:138836/1‑139 (MQ=255)
                                                                                      gCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCa                                                        <  1:25026/119‑1 (MQ=255)
                                                                                      gCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCa                                                        >  2:25026/1‑119 (MQ=255)
                                                                                      gCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGt                                     >  1:25512/1‑138 (MQ=255)
                                                                                       cccATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCa                                   >  1:22453/1‑139 (MQ=255)
                                                                                           tCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACg                                                          <  1:501462/112‑1 (MQ=255)
                                                                                           tCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACg                                                          >  2:501462/1‑112 (MQ=255)
                                                                                             cGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGc                              <  2:483594/138‑1 (MQ=255)
                                                                                             cGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGc                              <  2:360426/138‑1 (MQ=255)
                                                                                             cGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGc                              <  1:367391/138‑1 (MQ=255)
                                                                                             cGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGc                              <  1:5989/138‑1 (MQ=255)
                                                                                             cGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCa                             <  2:73786/139‑1 (MQ=255)
                                                                                               gAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCAcc                           <  1:99495/139‑1 (MQ=255)
                                                                                                aaGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCt                          >  1:489008/1‑139 (MQ=255)
                                                                                                 aGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCt                          >  2:485651/1‑138 (MQ=255)
                                                                                                     tACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACa                      >  1:216527/1‑138 (MQ=255)
                                                                                                     tACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACa                      >  1:548536/1‑138 (MQ=255)
                                                                                                       cGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACaa                     <  2:193041/137‑1 (MQ=255)
                                                                                                       cGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAg                   <  2:192125/139‑1 (MQ=255)
                                                                                                       cGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAg                   <  1:337196/139‑1 (MQ=255)
                                                                                                          ccTGATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGa                <  1:393893/139‑1 (MQ=255)
                                                                                                             gATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCAc                                                                                    <  1:597162/68‑1 (MQ=255)
                                                                                                             gATCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCAc                                                                                    >  2:597162/1‑68 (MQ=255)
                                                                                                               tCCCGGATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGaa           <  2:198251/139‑1 (MQ=255)
                                                                                                                    gATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCg      <  1:199667/139‑1 (MQ=255)
                                                                                                                    gATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCg      >  2:341280/1‑139 (MQ=255)
                                                                                                                    gATGTGGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCg      <  1:79025/139‑1 (MQ=255)
                                                                                                                        tgGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGa                                                      <  2:97470/87‑1 (MQ=255)
                                                                                                                        tgGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGa                                                      >  1:97470/1‑87 (MQ=255)
                                                                                                                        tgGTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCGTACt  >  1:452852/1‑139 (MQ=255)
                                                                                                                         ggTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCg                                         <  2:511838/99‑1 (MQ=255)
                                                                                                                         ggTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCg                                         >  1:511838/1‑99 (MQ=255)
                                                                                                                          gTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGa                                                      <  2:314451/85‑1 (MQ=255)
                                                                                                                          gTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGa                                                      >  1:314451/1‑85 (MQ=255)
                                                                                                                          gTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGt                                                >  1:618240/1‑91 (MQ=255)
                                                                                                                          gTTGAAGTGGCGAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGt                                                <  2:618240/91‑1 (MQ=255)
                                                                                                                                      |                                                                                                                             
CCTACTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCGTACT  >  NZ_CP009273/3000207‑3000466

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