Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,365,290 T→A 7.0% H65L (CAC→CTC)  PP_1189 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,365,2900TA 7.0% 167.7 / 3.7 57H65L (CAC→CTC) PP_1189hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/29);  new base A (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 = 8.78e-01

GCTCCACGCCCTCCAGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGCGGCTGGCGCGGAATGGCGAACTTTTCCTTGAAACA  >  NC_002947/1365147‑1365438
                                                                                                                                                |                                                                                                                                                    
gCTCCACGCCCTCCAGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTggaaca                                                                                                                                                <  1:1084360/150‑1 (MQ=255)
  tCCACGCCCTCCAGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTggaacagg                                                                                                                                              <  1:1223436/150‑1 (MQ=255)
    cACGCCCTCCAGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTggaacaggaa                                                                                                                                            >  2:521070/1‑150 (MQ=255)
      cGCCCTCCAGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTggaacaggaaca                                                                                                                                          <  1:900224/150‑1 (MQ=255)
         ccTCCAGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAAcagca                                                                                                                                       <  1:659093/150‑1 (MQ=255)
           tCCAGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACagcagc                                                                                                                                     >  1:1165055/1‑150 (MQ=255)
           cacggtcGTACCGCCGACTGGCCTATGCCGGTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACagcagc                                                                                                                                     <  2:80735/144‑1 (MQ=255)
              aGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGCCCTGGAggaaa                                                                                                                                                 >  2:991676/1‑134 (MQ=255)
              aGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGCCCTGGAggaaa                                                                                                                                                 <  1:991676/135‑2 (MQ=255)
                    accaccGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATgg                                                                                                                            <  1:765493/150‑1 (MQ=255)
                      caccGACTGGCCGATGCCGTTGGGCCGGTGGCTGG‑CGCGGGTGGCGAACACACCCATGCTCTTGTTGGCGCCCGGACGCGGCGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCCGGTGGAACGAGAACAGCAGCCACACATGGcg                                                                                                                          >  2:1365563/1‑149 (MQ=255)
                         cGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGAc                                                                                                                       <  2:1318438/150‑1 (MQ=255)
                         cGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCCGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGAc                                                                                                                       <  1:241048/150‑1 (MQ=255)
                                  gATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAg                                                                                                              <  1:439612/150‑1 (MQ=255)
                                   aTGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAgg                                                                                                             <  1:521070/150‑1 (MQ=255)
                                        gTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGAc                                                                                                                       >  2:301565/1‑135 (MQ=255)
                                        gTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGAc                                                                                                                       <  1:301565/135‑1 (MQ=255)
                                        gTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTc                                                                                                        >  2:845928/1‑150 (MQ=255)
                                        gTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCACCAGGCCCTGGAggaa                                                                                                                                                  >  2:652746/1‑108 (MQ=255)
                                        gTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCACCAGGCCCTGGAggaa                                                                                                                                                  <  1:652746/108‑1 (MQ=255)
                                           gggCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTggaacgggaa                                                                                                                                            <  1:924920/111‑1 (MQ=255)
                                           gggCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTggaacgggaa                                                                                                                                            >  2:924920/1‑111 (MQ=255)
                                           gggCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGAc                                                                                                     >  2:465988/1‑150 (MQ=255)
                                               tcgttggtggcCGCGGGTGTGGAACCCCCCCATCCTCTTGTTGCCGCCCAGCCCCGGTGGCCGCACTTTCAGCCGCCGCTTTTCCTCCAGGGCATGTTGGAACAGGAACAGCAGCCACACATCGCTGACCTGTTCCAGGCCCTCGACcgc                                                                                                  <  2:655161/145‑1 (MQ=255)
                                                      tgg‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCtt                                                                                          <  1:582868/150‑1 (MQ=255)
                                                        g‑cgcGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACagcagc                                                                                                                                     <  2:460406/105‑1 (MQ=255)
                                                        g‑cgcGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACagcagc                                                                                                                                     >  1:460406/1‑105 (MQ=255)
                                                        g‑cgcGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGATCTGTTCCAGGCCCTCGACCGCGTCACCTTGa                                                                                        <  2:872728/150‑1 (MQ=255)
                                                                   cGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGATGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAAcggcg                                                                              >  2:822827/1‑150 (MQ=255)
                                                                    gAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAAcggcgg                                                                             >  2:83733/1‑150 (MQ=255)
                                                                      acacCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTggaaca                                                                                                                                                >  2:716048/1‑81 (MQ=255)
                                                                      acacCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTggaaca                                                                                                                                                <  1:716048/81‑1 (MQ=255)
                                                                        acCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCAAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAAc                                                                         >  2:324246/1‑150 (MQ=255)
                                                                          ccccATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAg                                                                       >  1:1190900/1‑150 (MQ=255)
                                                                          ccccATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAg                                                                       >  2:522795/1‑150 (MQ=255)
                                                                            ccATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCAcc                                                                                            <  2:1182710/127‑1 (MQ=255)
                                                                            ccATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCAcc                                                                                            >  1:1182710/1‑127 (MQ=255)
                                                                            ccATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGtt                                                                     <  2:1165055/150‑1 (MQ=255)
                                                                             cATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCAAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTc                                                                    <  1:324246/150‑1 (MQ=255)
                                                                                gCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAg                                                                 <  1:522795/150‑1 (MQ=255)
                                                                                  tctTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGa                                                                                                      >  2:221652/1‑111 (MQ=255)
                                                                                  tctTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGa                                                                                                      <  1:221652/111‑1 (MQ=255)
                                                                                                 gACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAg                                                <  2:1060437/150‑1 (MQ=255)
                                                                                                     cgGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGgcgc                                            <  1:419739/150‑1 (MQ=255)
                                                                                                       gatctcGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTCCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCa                                          <  1:635953/145‑1 (MQ=255)
                                                                                                            cGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGc                                     <  2:821211/150‑1 (MQ=255)
                                                                                                            cGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCCCTTCCAGTCCCTCGCCCCGGGCACCTTGAGCGTACGGTGGTACCGGTGCGGGGTGGCCACAGGCGTCACGCCCCAGTTGc                                     >  2:1365915/1‑150 (MQ=255)
                                                                                                             gCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAATTCGAGTACGCCACAGGCGGCGTGCACCAgagatc                                    >  2:635953/1‑145 (MQ=255)
                                                                                                              cACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGCgg                                   <  1:465988/150‑1 (MQ=255)
                                                                                                                 tttCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGAGCCAGTTGCGGCTg                                >  2:303235/1‑150 (MQ=255)
                                                                                                                  ttCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGCGGCTgg                               >  1:1074123/1‑150 (MQ=255)
                                                                                                                         gcgGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGAc                                                                                                                       <  1:970807/55‑1 (MQ=255)
                                                                                                                         gcgGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGAc                                                                                                                       >  2:970807/1‑55 (MQ=255)
                                                                                                                             cTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGCGGCTGGCGCGGAATggg                    >  2:470119/1‑149 (MQ=255)
                                                                                                                              ttttCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGCGGCTGGCGCGGAATGGCg                   >  1:553053/1‑150 (MQ=255)
                                                                                                                                   ctccAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGCGGCTGGCGCGGAATGGc                    >  2:241778/1‑144 (MQ=255)
                                                                                                                                   ctccAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGCGGCTGGCGCGGAATGGc                    <  1:241778/144‑1 (MQ=255)
                                                                                                                                             tggtggAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGCGGCTGGCGCGGAATGGCGAACTTTTCCTTGaaa    <  2:553053/150‑1 (MQ=255)
                                                                                                                                               gtggAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAATGCCTCCGAGAGGTTGAGTACGCCACAGGCGGCAGGCGCCAGTTGCGGCTGGCGCGGAATGGCGAACTTTTCCTTGAAACa  >  2:1118726/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCTCCACGCCCTCCAGGCGTACCACCGACTGGCCGATGCCGTTGGGCCGGTGGGTGG‑CGCGGGTGGCGAACACCCCCATGCTCTTGTTGCCGCCCAGACGCGGAGGCCGCACTTTCAGGCGCGGCTTTTCCTCCAGGGCCTGGTGGAACAGGAACAGCAGCCACACATGGCTGACCTGTTCCAGGCCCTCGACCGCGTCACCTTGATCGAACGGCGGCAACAGTTCGAGTACGCCACAGGCGGCAGGCGCCAGTTGCGGCTGGCGCGGAATGGCGAACTTTTCCTTGAAACA  >  NC_002947/1365147‑1365438

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