Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,508,572 A→T 6.6% L70Q (CTG→CAG)  cysG ← multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,508,5720AT 6.6% 178.1 / 2.2 60L70Q (CTG→CAG) cysGmultifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base T (2/2);  total (28/32)
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.18e-01

CGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCT  >  NC_002947/4508425‑4508716
                                                                                                                                                   |                                                                                                                                                
cGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGAGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAg                                                                                                                                                 <  5:366067/149‑1 (MQ=255)
         aTCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGc                                                                                                                                        <  4:253557/149‑1 (MQ=255)
          tCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGcc                                                                                                                                       >  6:118123/1‑149 (MQ=255)
            gCCGGGAATATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAGCCGACCTGGTCGGTTGccg                                                                                                                                      <  6:260195/148‑2 (MQ=255)
               gggAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGa                                                                                                                                    >  5:119924/1‑147 (MQ=255)
               gggAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATa                                                                                                                                  <  5:196672/149‑1 (MQ=255)
                              cacaGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCa                                                                                                                    <  4:47557/148‑1 (MQ=255)
                              cacaGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCAc                                                                                                                   <  6:119924/149‑1 (MQ=255)
                               acaGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCAcc                                                                                                                  >  6:175606/1‑149 (MQ=255)
                                 aGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGa                                                                                                                                    >  3:316661/1‑129 (MQ=255)
                                 aGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGa                                                                                                                                    <  4:316661/129‑1 (MQ=255)
                                    gCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCAcc                                                                                                             <  5:118123/149‑1 (MQ=255)
                                             tCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTc                                                                                                    >  2:118545/1‑149 (MQ=255)
                                             tCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTc                                                                                                    >  4:215667/1‑149 (MQ=255)
                                                accacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCAcc                                                                                                                        <  3:1403/126‑1 (MQ=255)
                                                accacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCAcc                                                                                                                        >  4:1403/1‑126 (MQ=255)
                                                accacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCg                                                                                                 <  1:197523/149‑1 (MQ=255)
                                                 ccacATTGACCGGCAGGCTAAGTGCCTGCGCATCGGCCGACACCTGGGTATTCAGCGCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCATCCGACGAGGTCGGCTGCCTGATAGGCAGGCTGCAGGTGCTCACCGGCTCCTTCACgg                                                                                                <  8:281811/149‑1 (MQ=255)
                                                 ccacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCgg                                                                                                >  6:196537/1‑149 (MQ=255)
                                                     aTTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCt                                                                                                                                      >  4:274048/1‑107 (MQ=255)
                                                     aTTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCt                                                                                                                                      <  3:274048/107‑1 (MQ=255)
                                                                 gCTAAGGGCCTGCGCATCGGCCGACATCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACt                                                                                >  2:266995/1‑149 (MQ=255)
                                                                 gCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATa                                                                                                                                  >  6:313587/1‑99 (MQ=255)
                                                                 gCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATa                                                                                                                                  <  5:313587/99‑1 (MQ=255)
                                                                 gCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACt                                                                                >  6:212718/1‑149 (MQ=255)
                                                                   tAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTgg                                                                              >  8:116540/1‑149 (MQ=255)
                                                                            gcgcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCa                                                                                                                    >  6:100710/1‑102 (MQ=255)
                                                                            gcgcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCa                                                                                                                    <  5:100710/102‑1 (MQ=255)
                                                                               cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCTTCAGTCGCCCCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTcc                                                                  <  3:215667/149‑1 (MQ=255)
                                                                               cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTcc                                                                  <  1:118545/149‑1 (MQ=255)
                                                                                      cGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGa                                                                                                                                    >  8:129303/1‑76 (MQ=255)
                                                                                      cGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGa                                                                                                                                    <  7:129303/76‑1 (MQ=255)
                                                                                      cGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGcca                                                           <  2:251137/149‑1 (MQ=255)
                                                                                            ccgggCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCg                                                                 <  6:189058/135‑1 (MQ=255)
                                                                                            cTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCAcgc                                                     >  5:243039/1‑149 (MQ=255)
                                                                                             tGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCAcgcg                                                    >  4:289921/1‑149 (MQ=255)
                                                                                              gggCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTc                                                                         >  4:233402/1‑127 (MQ=255)
                                                                                              gggCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTc                                                                         <  3:233402/127‑1 (MQ=255)
                                                                                                   ttCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCAc                                              >  7:191551/1‑149 (MQ=255)
                                                                                                         ccAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACgccgcc                                        >  4:95787/1‑149 (MQ=255)
                                                                                                          cAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGccgccg                                       >  4:324188/1‑149 (MQ=255)
                                                                                                             ggTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCa                                    >  3:317395/1‑149 (MQ=255)
                                                                                                              gtcgtcAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCTGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCAt                                   <  6:243039/149‑1 (MQ=255)
                                                                                                                 gtcAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCgg                                >  2:35866/1‑149 (MQ=255)
                                                                                                                  tcAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGc                               <  5:175606/149‑1 (MQ=255)
                                                                                                                            cgATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGgcgc                     <  8:191551/149‑1 (MQ=255)
                                                                                                                            cgATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGgcgc                     <  7:116540/149‑1 (MQ=255)
                                                                                                                            cgATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGgcgc                     <  5:262061/149‑1 (MQ=255)
                                                                                                                                 aCCAGCCTGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCgcagc                                                >  5:337068/1‑117 (MQ=255)
                                                                                                                                 aCCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCgcagc                                                <  6:337068/117‑1 (MQ=255)
                                                                                                                                  ccAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTg               <  3:289921/149‑1 (MQ=255)
                                                                                                                                   cAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTgc              <  4:317395/149‑1 (MQ=255)
                                                                                                                                   cAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTgc              >  5:129425/1‑149 (MQ=255)
                                                                                                                                   cAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTgc              <  5:212718/149‑1 (MQ=255)
                                                                                                                                         gCAGCCGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGcc                                                            <  6:36895/97‑1 (MQ=255)
                                                                                                                                         gCAGCCGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGcc                                                            >  5:36895/1‑97 (MQ=255)
                                                                                                                                         gCAGCCGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAgg        >  1:194239/1‑149 (MQ=255)
                                                                                                                                           aaCCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCg      <  5:196537/149‑1 (MQ=255)
                                                                                                                                               gACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGgc                       <  8:269575/128‑1 (MQ=255)
                                                                                                                                                aCCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGAtct  >  5:67303/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCT  >  NC_002947/4508425‑4508716

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: