Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,204,071 C→T 7.0% R250W (CGG→TGG)  PP_1948 → benzaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,204,0710CT 7.0% 169.8 / 2.3 57R250W (CGG→TGG) PP_1948benzaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base T (2/2);  total (26/31)
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.26e-01

CCAGTTCGCGGAACTGGTTGTGAAAGAGGCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCACCGTCAACACTGGTC  >  NC_002947/2203923‑2204215
                                                                                                                                                     |                                                                                                                                                
ccAGTTCGCGGAACTGGTTGTGAAAGAGGCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTACGTGCCGGAc                                                                                                                                                  >  3:71020/1‑149 (MQ=255)
   gTTCGCGGAACTGGTTGTGAAAGAGGCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACgg                                                                                                                                                >  2:66621/1‑148 (MQ=255)
         ggAACTGGTTGTGAAAGAGGCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGTCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATc                                                                                                                                          <  1:54178/149‑1 (MQ=255)
             cTGGTTGTGAAAGAGGCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCaa                                                                                                                                     <  6:62692/149‑1 (MQ=255)
             cTGGTTGTGAAAGAGGCAGGGTTGCCACCAGGCGTATTAAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCaa                                                                                                                                     <  7:270062/149‑1 (MQ=255)
                    tgAAAGAGGCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTcgcc                                                                                                                              <  2:134552/149‑1 (MQ=255)
                         gagGCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCg                                                                                                                                         <  7:199698/133‑1 (MQ=255)
                         gagGCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCg                                                                                                                                         >  8:199698/1‑133 (MQ=255)
                            gCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCg                                                                                                                      <  3:318542/149‑1 (MQ=255)
                                  ttGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGcatcat                                                                                                                >  1:256124/1‑149 (MQ=255)
                                    gCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTc                                                                                                              >  5:178626/1‑149 (MQ=255)
                                     ccaccaGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTcc                                                                                                             >  4:258943/1‑149 (MQ=255)
                                         caGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTg                                                                                                         >  2:234090/1‑149 (MQ=255)
                                          aGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGt                                                                                                        <  3:234820/149‑1 (MQ=255)
                                            gCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCg                                                                                                      <  5:355008/149‑1 (MQ=255)
                                            gCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCg                                                                                                      <  3:352190/149‑1 (MQ=255)
                                                    aaCGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGATGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAAtt                                                                                              >  3:175343/1‑149 (MQ=255)
                                                    aaCGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGATGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAAtt                                                                                              >  3:180943/1‑149 (MQ=255)
                                                      cGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTgg                                                                                            <  6:178626/149‑1 (MQ=255)
                                                        tGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGggc                                                                                          <  6:70906/149‑1 (MQ=255)
                                                              ccGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCaaa                                                                                    <  3:258943/149‑1 (MQ=255)
                                                               cGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAAt                                                                                   >  7:238079/1‑149 (MQ=255)
                                                                       gTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGaa                                                                            >  8:121211/1‑148 (MQ=255)
                                                                       gTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGaa                                                                            >  3:175613/1‑148 (MQ=255)
                                                                           ggAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTg                                                                       <  1:234090/149‑1 (MQ=255)
                                                                                 cGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTTGGCGGCAAATCCCCGAACATTGTCTTCg                                                                 <  2:256124/149‑1 (MQ=255)
                                                                                 cGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTc                                                                  >  7:25915/1‑148 (MQ=255)
                                                                                   gTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGaa                                                               <  8:238079/149‑1 (MQ=255)
                                                                                          tctGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTc                                                                                                   >  8:128084/1‑106 (MQ=255)
                                                                                          tctGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTc                                                                                                   <  7:128084/106‑1 (MQ=255)
                                                                                           ctGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTc                                                                  <  8:125173/138‑1 (MQ=255)
                                                                                           ctGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTc                                                                  >  7:125173/1‑138 (MQ=255)
                                                                                              gTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCtt                                                                   <  2:17055/134‑1 (MQ=255)
                                                                                              gTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCtt                                                                   >  1:17055/1‑134 (MQ=255)
                                                                                                  aaCACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCTAATCCCCAAACATTGTCTTCGAAGACGCAGATCTGGa                                                 <  1:66621/148‑1 (MQ=255)
                                                                                                  aaCACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGAt                                                >  1:186805/1‑149 (MQ=255)
                                                                                                             aTCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTc                                       <  4:23077/147‑1 (MQ=255)
                                                                                                              tCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTcgcg                                    >  5:175812/1‑149 (MQ=255)
                                                                                                                  aaaGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTcgc                                     <  4:295462/144‑1 (MQ=255)
                                                                                                                  aaaGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTcgc                                     >  3:295462/1‑144 (MQ=255)
                                                                                                                  aaaGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTa                                <  6:328415/149‑1 (MQ=255)
                                                                                                                    aGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGc                              <  7:121211/149‑1 (MQ=255)
                                                                                                                            tttCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTg                       >  3:185793/1‑149 (MQ=255)
                                                                                                                                aCCGGTTCTGTCCGTGCCGGATGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCtt                   <  4:175343/149‑1 (MQ=255)
                                                                                                                                aCCGGTTCTGTCCGTGCCGGATGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCtt                   <  4:180943/149‑1 (MQ=255)
                                                                                                                                aCCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCtt                   >  7:238110/1‑149 (MQ=255)
                                                                                                                                  cGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTCTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCa                 >  3:206247/1‑149 (MQ=255)
                                                                                                                                     ttCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCACCg              >  7:305115/1‑149 (MQ=255)
                                                                                                                                       cTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCACCgc             <  4:71020/148‑2 (MQ=255)
                                                                                                                                         gTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGaa                                                               >  4:233423/1‑96 (MQ=255)
                                                                                                                                         gTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGaa                                                               <  3:233423/96‑1 (MQ=255)
                                                                                                                                          tCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCACCGTCaa          >  4:207746/1‑148 (MQ=255)
                                                                                                                                            cGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCACCGTCAaca        <  4:185793/148‑1 (MQ=255)
                                                                                                                                             gTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGa                                             >  2:316340/1‑110 (MQ=255)
                                                                                                                                             gTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGa                                             <  1:316340/110‑1 (MQ=255)
                                                                                                                                             gTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCACCGTCAACACt      <  7:212660/149‑1 (MQ=255)
                                                                                                                                                ccGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCACCGTCAACACTGGt   >  4:47215/1‑149 (MQ=255)
                                                                                                                                                 cGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTCTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCACCGTCAACACTGGTc  <  4:206247/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CCAGTTCGCGGAACTGGTTGTGAAAGAGGCAGGGTTGCCACCAGGCGTATTGAACGTGGTTACCGGCACCGGTAAGGAAACCGGTGAGCCTCTGGTTAAACACCCTCTGATCCGAAAGGTTGCTTTCACCGGTTCTGTCCGTGCCGGACGGGAGATCGGCAAGCTCGCCGCAGATCGCATCATTCCGCTGTCGCTCGAATTGGGCGGCAAATCCCCGAACATTGTCTTCGAAGACGCAGATCTGGATCGAGCTGTCGCGGGTAGCGTCTTTGCCTTCACCGTCAACACTGGTC  >  NC_002947/2203923‑2204215

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