Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,184,152 A→T 7.0% V78E (GTA→GAA)  mltF ← membrane‑bound lytic murein transglycosylase F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,184,1520AT 7.0% 173.9 / 2.9 57V78E (GTA→GAA) mltFmembrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base T (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

ACCTCCAGATAGGGATGGGAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGG  >  NC_002947/1184005‑1184293
                                                                                                                                                   |                                                                                                                                               
aCCTCCAGATAGGGATGGGAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTAc                                                                                                                                                >  6:262401/1‑149 (MQ=255)
    ccAGATAGGGATGGGAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTAGGCGGTCTCGATCTGCAGCTTTACGcc                                                                                                                                             >  8:267640/1‑148 (MQ=255)
         tAGGTATGGGAGTATCTCACCTGGTTGTTGCGCCGCTCACTGCTGACCAGGCGGGCTGCAGCCAGTACCGGTCCGGAGGGCTTTCCGAGTGCGTCTTACAGCTCGTCCTGGTTGTCGGCGGTCTCTATCTGCAGCTTTACGCCGAGAtc                                                                                                                                      <  6:169066/149‑1 (MQ=255)
                  gAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGa                                                                                                                              >  8:84580/1‑148 (MQ=255)
                   aGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGaaa                                                                                                                            >  8:31601/1‑149 (MQ=255)
                         tcACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGGTg                                                                                                                       >  2:301737/1‑148 (MQ=255)
                         tcACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCATTTCGCCGACGATCTCGGCGTAAAGCTgc                                                                                                                      >  4:47320/1‑148 (MQ=255)
                           aCCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                       <  7:84580/146‑1 (MQ=255)
                           aCCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCATTTCGCCGACGATCTCGGCGTAAAGCTgcag                                                                                                                    <  3:47320/149‑4 (MQ=38)
                           aCCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCATTTCGCCGACGATCTCGGCGTAAAGCTgc                                                                                                                      <  2:55713/147‑2 (MQ=255)
                             cTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGTCAGCCAGTACCGGGACGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGATTTACGCCGAGATCGTCGGCGAAATGCTGGAc                                                                                                                    >  2:203900/1‑147 (MQ=255)
                             cTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGaaa                                                                                                                            <  4:134000/139‑1 (MQ=255)
                             cTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGaaa                                                                                                                            >  3:134000/1‑139 (MQ=255)
                             cTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGAc                                                                                                                    >  1:110513/1‑147 (MQ=255)
                                ggTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGGTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGtt                                                                                                               <  8:257024/149‑1 (MQ=255)
                                               aCTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACTGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCgg                                                                                               >  1:89708/1‑149 (MQ=255)
                                               aCTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                       <  6:90200/126‑1 (MQ=255)
                                               aCTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                       >  5:90200/1‑126 (MQ=255)
                                               aCTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCg                                                                                                <  7:31601/148‑1 (MQ=255)
                                                 tgctgACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                       >  3:305543/1‑124 (MQ=255)
                                                 tgctgACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                       <  4:305543/124‑1 (MQ=255)
                                                      aCCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGcc                                                                                        >  3:252493/1‑149 (MQ=255)
                                                      aCCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGc                                                                                         <  5:286258/148‑1 (MQ=255)
                                                                  gcAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTg                                                                            <  8:161746/149‑1 (MQ=255)
                                                                    aGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCAGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                           <  8:150811/148‑1 (MQ=255)
                                                                    aGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCAGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                           >  7:150811/1‑148 (MQ=255)
                                                                     gCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGaa                                                                         <  7:28952/149‑1 (MQ=255)
                                                                          tACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCg                                                                                          >  4:188666/1‑127 (MQ=255)
                                                                          tACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCg                                                                                          <  3:188666/127‑1 (MQ=255)
                                                                                   ggAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGc                                                           <  5:262401/149‑1 (MQ=255)
                                                                                         cTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGc                                                     <  2:110513/149‑1 (MQ=255)
                                                                                          ttGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCg                                                    >  3:138760/1‑149 (MQ=255)
                                                                                             ccGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCgtgt                                                 <  4:252493/149‑1 (MQ=255)
                                                                                             ccGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCgtgt                                                 <  1:203900/149‑1 (MQ=255)
                                                                                               gAGGGCGTCGTACTGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGa                                               <  2:89708/149‑1 (MQ=255)
                                                                                                 gggCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATc                                             <  8:40531/149‑1 (MQ=255)
                                                                                                       cgtACAGTTCGTCTAGGTTGTAGGCTGTCTCGATCTGCAGCTTTAGGCCTAGATCGTCGGCGAAATGCTTTTCCTGTTCGTACTCAAAGCCTTTTTCGCCGTTACGTTCCTGTAACTATGTTGCCGGGCTGTTGCGTGTGATCActcg                                        <  4:174973/148‑1 (MQ=255)
                                                                                                                 gTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCAt                             >  3:237291/1‑149 (MQ=255)
                                                                                                                      ggTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCAtcctcc                        >  8:206642/1‑149 (MQ=255)
                                                                                                                       gTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATcctcct                       >  3:98096/1‑149 (MQ=255)
                                                                                                                          gTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCa                              >  1:235437/1‑139 (MQ=255)
                                                                                                                          gTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCa                              <  2:235437/139‑1 (MQ=255)
                                                                                                                              gcggTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTc                                                                                           <  8:215009/74‑1 (MQ=255)
                                                                                                                              gcggTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTc                                                                                           >  7:215009/1‑74 (MQ=255)
                                                                                                                                 gTCTCGATCTTCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGATCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                           <  1:145745/87‑1 (MQ=255)
                                                                                                                                 gTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGATCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                           >  2:145745/1‑87 (MQ=255)
                                                                                                                                  tctcGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTc             >  2:79899/1‑148 (MQ=255)
                                                                                                                                   ctcGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTc             >  3:105770/1‑147 (MQ=255)
                                                                                                                                    tcGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAg          >  7:10271/1‑149 (MQ=255)
                                                                                                                                      gATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCTTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAgg         <  7:267640/148‑1 (MQ=255)
                                                                                                                                      gATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAggg        >  4:263508/1‑149 (MQ=255)
                                                                                                                                      gATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAggg        >  4:259130/1‑149 (MQ=255)
                                                                                                                                      gATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAgg         <  7:206642/148‑1 (MQ=255)
                                                                                                                                            cAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAggg        <  8:10271/143‑1 (MQ=255)
                                                                                                                                            cAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAggg        <  4:98096/143‑1 (MQ=255)
                                                                                                                                            cAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTgg  >  3:157399/1‑149 (MQ=255)
                                                                                                                                             aGCATTTCGCCGACGATCTCGGCGTAAAGCTGGACCAGTTCGTACTCAAAGCCTAGTTTCGCCGTTACGGTCATGGAAGTAGGTTGCCGGCCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTgg  >  1:241875/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACCTCCAGATAGGGATGGGAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGG  >  NC_002947/1184005‑1184293

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: