Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,545,358 C→A 8.5% K164N (AAG→AAT hyuC ← bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,545,3580CA 8.5% 187.1 / 5.7 59K164N (AAG→AAThyuCbifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base A (2/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

TTGGCCTGGTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGCCCAGGTCGTTGAGGGTGCGG  >  NC_002947/4545218‑4545505
                                                                                                                                             |                                                                                                                                                   
tgggCCTGATCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGGAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGAATCACGCTTGGCCAgg                                                                                                                                            >  2:161543/3‑150 (MQ=255)
 tGGCCTGGTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAggg                                                                                                                                           <  2:1367579/150‑1 (MQ=255)
       ggTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCtt                                                                                                                                     >  1:722121/1‑150 (MQ=255)
        gTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTc                                                                                                                                    >  2:1003347/1‑150 (MQ=255)
          ccTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCa                                                                                                                                  <  1:1003347/150‑1 (MQ=255)
            tCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGTACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                                >  2:525526/1‑150 (MQ=255)
            tCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                                <  1:820601/150‑1 (MQ=255)
            tCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                                <  1:601078/150‑1 (MQ=255)
             cAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAggg                                                                                                                               >  1:985510/1‑150 (MQ=255)
                ggATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGa                                                                                                                            >  1:480782/1‑150 (MQ=255)
                      gCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGc                                                                                                                      >  1:51099/1‑150 (MQ=255)
                      gCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGc                                                                                                                      >  1:498341/1‑150 (MQ=255)
                           tGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCgg                                                                                                                   <  2:1264950/148‑1 (MQ=255)
                           tGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCgg                                                                                                                   >  1:1264950/1‑148 (MQ=255)
                              tcgaagtggtggtgtggaccaacgaagaaggttcgcgttttgcgccatgcatgatggggtctggggtattcgccggcaagttcaccctggaaGACACCCTGGGCATCGCGCTTGGCCAGGgt                                                                                                                                          >  2:69006‑M1/93‑122 (MQ=255)
                              tcgaagtggtggtgtggaccaacgaagaaggttcgcgttttgcgccatgcatgatggggtctggggtattcgccggcaagttcaccctggaaGACACCCTGGGCATCGCGCTTGGCCAGGgt                                                                                                                                          <  1:69006‑M1/30‑1 (MQ=255)
                                gATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAgg                                                                                                                                            <  1:317159/118‑1 (MQ=255)
                                gATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAgg                                                                                                                                            >  2:317159/1‑118 (MQ=255)
                                  tatGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATAccc                                                                                                          <  2:51099/150‑1 (MQ=255)
                                    tGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCa                                                                                                        <  1:406405/150‑1 (MQ=255)
                                    tGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCa                                                                                                        <  2:722121/150‑1 (MQ=255)
                                         ttCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATAccc                                                                                                          >  2:1090749/1‑143 (MQ=255)
                                         ttCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATAccc                                                                                                          <  1:1090749/143‑1 (MQ=255)
                                         ttCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGAccc                                                                                                   >  2:734286/1‑150 (MQ=255)
                                                tACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatg                                                                                            >  1:1248640/1‑150 (MQ=255)
                                                   gCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcat                                                                                         >  1:1370800/1‑150 (MQ=255)
                                                     ccccACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCTCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcatg                                                                                        >  1:754694/1‑149 (MQ=255)
                                                     ccccACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCTCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcatg                                                                                        <  2:754694/149‑1 (MQ=255)
                                                     ccccACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATgg                                                                                       >  1:129580/1‑150 (MQ=255)
                                                                   ccGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                                <  1:231306/95‑1 (MQ=255)
                                                                   ccGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                                >  2:231306/1‑95 (MQ=255)
                                                                    cGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCa                                                                                                                                  >  2:451866/1‑92 (MQ=255)
                                                                    cGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCa                                                                                                                                  <  1:451866/92‑1 (MQ=255)
                                                                             aCGAGGACCGGCGTAGCCAATGGAGTGCAACGCTTCGCCAACGCTGACACCCTGAGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGCCCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGtt                                                               <  1:655116/150‑1 (MQ=255)
                                                                              cgcgGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTg                                                              <  2:480782/150‑1 (MQ=255)
                                                                               gcgGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAggg                                                                                                                                           <  1:1147700/72‑1 (MQ=255)
                                                                               gcgGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAggg                                                                                                                                           >  2:1147700/1‑72 (MQ=255)
                                                                                     cGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTc                                                           <  1:1032423/146‑1 (MQ=255)
                                                                                     cGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTc                                                           >  2:1032423/1‑146 (MQ=255)
                                                                                       gCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCAcacc                                                     >  2:1264651/1‑150 (MQ=255)
                                                                                        cGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCAcacca                                                    <  2:129580/150‑1 (MQ=255)
                                                                                         gTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGTCGCAAAACGCGAACCTTCTTCGTTGGTCCAcaccac                                                   <  2:985510/150‑1 (MQ=255)
                                                                                            gCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGaa                                                                                                               >  2:1195735/1‑87 (MQ=255)
                                                                                            gCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGaa                                                                                                               <  1:1195735/87‑1 (MQ=255)
                                                                                                atgtggTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCGTCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCg                                            >  2:1066355/6‑150 (MQ=255)
                                                                                                 tGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGa                                           <  1:91336/150‑1 (MQ=255)
                                                                                                 tGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGa                                           <  1:620236/150‑1 (MQ=255)
                                                                                                             gttcaccctggaaGACACCCTGGCCAAGCGCGATGCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcat                                                                                         <  2:535795‑M1/79‑1 (MQ=255)
                                                                                                             gttcaccctggaaGACACCCTGGCCAAGCGCGATGCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcat                                                                                         >  1:535795‑M1/14‑92 (MQ=255)
                                                                                                              ttCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCacac                                                      <  2:1032840/126‑1 (MQ=255)
                                                                                                              ttCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCacac                                                      >  1:1032840/1‑126 (MQ=255)
                                                                                                                     tggaaGACACCCTGGCCAAGCGCGATGCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTg                                                              >  2:476964‑M1/6‑111 (MQ=255)
                                                                                                                     tggaaGACACCCTGGCCAAGCGCGATGCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTg                                                              <  1:476964‑M1/106‑1 (MQ=255)
                                                                                                                     tggaaGACACCCTGGCCAAGCGCGATGCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACg                       <  2:182926‑M1/145‑1 (MQ=255)
                                                                                                                      cGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGc                      >  1:1103206/1‑150 (MQ=255)
                                                                                                                             ACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGCCCGACGCTTCGGTTTCTACCCCCAGGTc               <  2:139348‑M1/150‑1 (MQ=255)
                                                                                                                              CCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCTACGCTTCGGTTTCCACGCCCAGGTCg              <  2:1008308‑M1/150‑1 (MQ=255)
                                                                                                                                   GGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGCCCAGGTCGTTGAg         >  1:734286‑M1/1‑150 (MQ=255)
                                                                                                                                         cgcgCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGCCCAGGTCGTTGAGGGTGc    <  2:498341/150‑1 (MQ=255)
                                                                                                                                          gcgcTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTc                                                           >  1:1358330/1‑94 (MQ=255)
                                                                                                                                          gcgcTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTc                                                           <  2:1358330/94‑1 (MQ=255)
                                                                                                                                           cgcTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGGATACCGCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGCCCAGGTCGTTGAGGGTGCgg  >  2:984743/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TTGGCCTGGTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCATCGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGCCCAGGTCGTTGAGGGTGCGG  >  NC_002947/4545218‑4545505

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