Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,239,889 T→C 8.6% P142P (CCA→CCG ldh ← leucine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,239,8890TC 8.6% 158.0 / 5.3 58P142P (CCA→CCGldhleucine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (25/28);  new base C (2/3);  total (27/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.22e-01

GCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAG  >  NC_002947/5239745‑5240027
                                                                                                                                                |                                                                                                                                          
gcCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGc                                                                                                                                         >  8:7007/1‑148 (MQ=255)
  ccACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGAt                                                                                                                                      <  6:173462/149‑1 (MQ=255)
     ccGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATggg                                                                                                                                   <  6:61646/149‑1 (MQ=255)
      cGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCACATGCGGCGATgggg                                                                                                                                  <  7:283041/149‑2 (MQ=255)
      cGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCACATGCGGCGATgggg                                                                                                                                  <  3:320719/149‑2 (MQ=255)
       gCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACAGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTc                                                                                                                                 <  6:352763/149‑1 (MQ=255)
         aTGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGc                                                                                                                               >  6:287072/1‑149 (MQ=255)
                  cTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCa                                                                                                                       <  8:152751/148‑1 (MQ=255)
                      gCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGAGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgtggt                                                                                                                  >  1:307243/1‑149 (MQ=255)
                        cAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGc                                                                                                                <  5:136867/149‑1 (MQ=255)
                                  ccAACATTGTCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgtg                                                                                                                    >  4:296157/1‑135 (MQ=255)
                                  ccAACATTGTCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgtg                                                                                                                    <  3:296157/135‑1 (MQ=255)
                                  ccAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACg                                                                                                        >  4:279212/1‑147 (MQ=255)
                                  ccAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTa                                                                                                      <  7:7007/149‑1 (MQ=255)
                                    aaCATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAgg                                                                                                    <  3:146203/149‑1 (MQ=255)
                                           cccAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCt                                                                                             >  2:37350/1‑149 (MQ=255)
                                               aCCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                         <  3:279212/149‑1 (MQ=255)
                                                   cTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTggggga                                                                                     >  7:181952/1‑146 (MQ=255)
                                                   cTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                         <  8:181952/145‑1 (MQ=255)
                                                   cTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCaa                                                                                     >  5:288408/1‑149 (MQ=255)
                                                   cTGGACCGCTACCCCCAAGCCTCGCCGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCaa                                                                                     <  4:64043/149‑1 (MQ=255)
                                                        ccGCTACCCGCAAGCCTCTCATGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAg                                                                                <  6:174031/149‑1 (MQ=255)
                                                        ccGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAg                                                                                <  1:37350/149‑1 (MQ=255)
                                                              cccGCGAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCANGTAGGGCGTGCTGTGGGCAATGCAGTCCAt                                                                           <  1:34971/148‑1 (MQ=255)
                                                                       ctctCAGGTCTTCGCTGCCCACCCGGAACGTGGTTGTGGCACGTAGGCCCGCTACCACCCCCATCGCCGCTTGTGGGGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTc                                                                 <  2:178159/149‑1 (MQ=255)
                                                                       ctctCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATg                                                                                   >  3:173019/1‑131 (MQ=255)
                                                                       ctctCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATg                                                                                   <  4:173019/131‑1 (MQ=255)
                                                                        tctcAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTc                                                                 >  1:73728/1‑148 (MQ=255)
                                                                            aGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGt                                                            <  4:235528/149‑1 (MQ=255)
                                                                              gTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGc                                                          >  3:176075/1‑149 (MQ=255)
                                                                              gTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGc                                                          >  8:264435/1‑149 (MQ=255)
                                                                               tCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                              >  5:329889/1‑144 (MQ=255)
                                                                               tCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                              <  6:329889/144‑1 (MQ=255)
                                                                                 aTCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCAc                                                       >  6:130248/1‑149 (MQ=255)
                                                                                     cggCCCAGCCGGATCGTGGTTGGGGCACGTTTGCCCGCGAACCCCCCCATCGGCGCAGGTGGCTATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTc                                                   <  2:307243/147‑1 (MQ=255)
                                                                                        cccAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCa                                                                                 >  1:157461/1‑116 (MQ=255)
                                                                                        cccAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCa                                                                                 <  2:157461/116‑1 (MQ=255)
                                                                                         ccAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACg                                               >  8:66708/1‑149 (MQ=255)
                                                                                               ggAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgatt                                         <  7:264435/149‑2 (MQ=255)
                                                                                                       gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgcg                                 >  3:160578/1‑149 (MQ=255)
                                                                                                            ggCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCAt                            >  2:17910/1‑149 (MQ=255)
                                                                                                             gCACGTATGCCCGAGAACACCCACATCGCCGATTGTGGCTATGGGTCGCCCGCTGAAGTGGTTCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCAt                            <  3:243525/148‑1 (MQ=255)
                                                                                                                 gTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAg                       <  6:288408/149‑1 (MQ=255)
                                                                                                                 gTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAg                       <  5:287072/149‑1 (MQ=255)
                                                                                                                      cccGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCACGGTATc                  >  3:206493/1‑149 (MQ=255)
                                                                                                                       ccGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAgg                      >  7:60549/1‑144 (MQ=255)
                                                                                                                       ccGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAgg                      <  8:60549/144‑1 (MQ=255)
                                                                                                                               acCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCg         <  1:17910/149‑1 (MQ=255)
                                                                                                                               acCCCCATCGCCGCATGTGGCAATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCg         >  2:347405/1‑149 (MQ=255)
                                                                                                                                 ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGc       >  2:364365/1‑149 (MQ=255)
                                                                                                                                  cccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGc       >  3:309171/1‑148 (MQ=255)
                                                                                                                                   ccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGcc      >  8:206552/1‑148 (MQ=255)
                                                                                                                                   ccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTAt                   <  2:114530/135‑1 (MQ=255)
                                                                                                                                   ccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTAt                   >  1:114530/1‑135 (MQ=255)
                                                                                                                                   ccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGcc      >  4:117606/1‑148 (MQ=255)
                                                                                                                                    cATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACgg                                              >  5:119043/1‑107 (MQ=255)
                                                                                                                                    cATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACgg                                              <  6:119043/107‑1 (MQ=255)
                                                                                                                                      tCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAg  <  2:73728/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                          
GCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAG  >  NC_002947/5239745‑5240027

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