Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,547,609 G→A 6.9% intergenic (‑264/‑80) pydP ← / → pydB NCS1 family transporter PydP/bifunctional D‑hydantoinase/dihydropyrimidinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,547,6090GA 6.9% 160.6 / 2.0 59intergenic (‑264/‑80)pydP/pydBNCS1 family transporter PydP/bifunctional D‑hydantoinase/dihydropyrimidinase
Reads supporting (aligned to +/- strand):  ref base G (29/24);  new base A (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

AAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTACCC  >  NC_002947/4547465‑4547741
                                                                                                                                                 |                                                                                                                                      
aaaTACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGAtt                                                                                                                                    >  3:347456/1‑149 (MQ=255)
aaaTACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAtt                                                                                                                                    >  5:383667/1‑149 (MQ=255)
   tACCCGGACAACCCCGGCGCCCACTCGCCCGCTAATTTATTCCTATTGATATCAACACCTTAAAAACCACAAATTTAATCGGAAAATTTTCTTGCTGACTCCAGAATCCGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcga                                                                                                                                 <  8:112826/149‑4 (MQ=255)
   tACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcga                                                                                                                                 <  6:383667/149‑4 (MQ=255)
      cAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgc                                                                                                                              >  2:168418/1‑149 (MQ=255)
      cAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgc                                                                                                                              >  7:310375/1‑149 (MQ=255)
        ggAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccc                                                                                                                            >  2:273093/1‑149 (MQ=255)
        gcaaacCCCCGCCCCCCATCCCCCGGTAAATTTTTCCCAATTTATACAAAAACCTTAAAAACACAAAATTAAATGGCAAATTTTCCTGCCTTATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccct                                                                                                                           <  4:360733/147‑1 (MQ=255)
                        cATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATcc                                                                                                           <  1:347800/149‑1 (MQ=255)
                             gcgggTTAAATTTTTTCTTTTTATAATAACAACTTTAAAACCACAAATTTTATACATTAAATTAGATATCTTATTCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGa                                                                                                       <  1:320723/149‑1 (MQ=255)
                                cgTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAAt                                                                                                   <  2:172396/149‑1 (MQ=255)
                                 gTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATg                                                                                                  >  4:129524/1‑149 (MQ=255)
                                 gTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAACGCCCTGCCCGCCCTGCATCCATGACAATg                                                                                                  >  4:45556/1‑149 (MQ=255)
                                      ttttttCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCGCTGCATCCATGACAATGCa                                                                                                >  2:267906/1‑146 (MQ=255)
                                            ccTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCg                                                                                        <  1:273093/148‑1 (MQ=255)
                                             cTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCgg                                                                                       <  8:310375/148‑1 (MQ=255)
                                                  gATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTTTAAGGACTGACATGATTGACCGCCCTGCCCGCCCTGCATCAATGACAATGCAAAGAACCGGCCAAGa                                                                                 >  3:209975/1‑149 (MQ=255)
                                                   atatCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGAc                                                                                <  6:176439/149‑1 (MQ=255)
                                                    tatCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGAcc                                                                               >  7:89894/1‑149 (MQ=255)
                                                      tCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCg                                                                              <  2:513/148‑1 (MQ=255)
                                                        aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGt                                                                            >  1:196398/1‑148 (MQ=255)
                                                        aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGt                                                                            <  3:45556/148‑1 (MQ=255)
                                                        aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTc                                                                           <  2:196398/149‑1 (MQ=255)
                                                        aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATAAGCATCTACACACGAATCTTGTCAGGACTGACAGGATTGACCGCCCTTCACCACCTTCATCCATGACAAACCAAAGAACCGGACAAGACCgga                                                                            >  1:236673/1‑147 (MQ=255)
                                                         acaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCa                                                                          <  8:126548/149‑1 (MQ=255)
                                                           aacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAg                                                                         >  8:332272/1‑148 (MQ=255)
                                                              ttAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTTACAGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTg                                                                     >  7:143659/1‑149 (MQ=255)
                                                              ttAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTg                                                                     >  2:140672/1‑149 (MQ=255)
                                                                 aaaaCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGc                                                                    >  6:282399/1‑147 (MQ=255)
                                                                    aCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAgg                                                                >  7:317961/1‑148 (MQ=255)
                                                                      cacaAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCGAAGAACCGGCCAAGACCGGTCAGCCTGCGa                                                                  >  3:287775/1‑144 (MQ=255)
                                                                      cacaAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGaaa                                                            <  4:347456/149‑1 (MQ=255)
                                                                      cacaAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGaa                                                             >  7:98400/1‑148 (MQ=255)
                                                                       acaAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGaaaa                                                           >  3:232388/1‑149 (MQ=255)
                                                                           aTTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGc                                                       <  7:154616/149‑1 (MQ=255)
                                                                                 tGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTcc                                                 <  7:332272/149‑1 (MQ=255)
                                                                                    ggAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGGCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCt                                               >  2:244351/1‑149 (MQ=255)
                                                                                     gAAAATTTTCTTGCTGAATCCAGAATCATCATCTAGCCTCGAATCTTGTCAGGACTGACATGATTGACCGCCCTGCCCGCCCTGCTTCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCGCt                                               <  2:236673/148‑1 (MQ=255)
                                                                                     gAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGACCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTg                                              >  2:250902/1‑149 (MQ=255)
                                                                                            ttCTTGCTGAATGCAGAATCAGCATCTAGCCTCGAATGTTGTCAGGACTTACAGGATTGAACGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGACAAGACCGGTCAGCCTGCTAGGAAAAGGGCATGTCCCTGTTGATc                                        <  8:143659/148‑1 (MQ=255)
                                                                                            ttCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATc                                        <  1:168418/148‑1 (MQ=255)
                                                                                                 gCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGgc                                  >  5:129588/1‑149 (MQ=255)
                                                                                                  cTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGgcg                                 >  6:107889/1‑149 (MQ=255)
                                                                                                              tCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCGTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTa                     <  1:162048/149‑1 (MQ=255)
                                                                                                               cAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGTCCCGCCCTGCATCCATGACAATGCAAAGACCCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTa                     >  1:91539/1‑149 (MQ=255)
                                                                                                               cAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTCACCGCCCTGACCGCCCTGCATGAATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCACTGTTGATCCGTGCCGCCACCGTGGTTAc                    >  4:274717/1‑149 (MQ=255)
                                                                                                                      tAGGCTCGAATCTTGTCAGGAATGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGGAAAAGGAAAGAACCGGCCAAGACCGGTCAGCCGGGGAGGAAAACGGCATGGGCCGGTTGCTCCGTGGCGCCACCGTGGTTACCCACg               >  8:192287/1‑147 (MQ=255)
                                                                                                                        gCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga           <  4:232388/149‑1 (MQ=255)
                                                                                                                        gCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga           <  6:129588/149‑1 (MQ=255)
                                                                                                                        gCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga           <  3:129524/149‑1 (MQ=255)
                                                                                                                        gCCTCGAATCTTGTCACGACTGACAGGTATGACCGCCCAGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGTAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGaa             <  1:291339/148‑1 (MQ=255)
                                                                                                                         ccTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCACTGTTGATCCGTGGCGCCACCGTGGTTACCCACGaa             >  4:184285/1‑146 (MQ=255)
                                                                                                                          cTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagaga         <  3:84705/149‑1 (MQ=255)
                                                                                                                          cTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagaga         <  2:328851/149‑1 (MQ=255)
                                                                                                                          cTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga           >  3:264986/1‑147 (MQ=255)
                                                                                                                           tCGAATCTTGTCAGGACTGACAGGATTAACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagagag        >  2:44579/1‑149 (MQ=255)
                                                                                                                           gtcaaTCCTGTCAGTCCTGACAAGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagagag        >  4:344517/4‑149 (MQ=255)
                                                                                                                             gAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGt       >  5:299004/1‑148 (MQ=255)
                                                                                                                                    gTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTAccc  <  8:89894/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                      
AAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTACCC  >  NC_002947/4547465‑4547741

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