Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,623,585 A→G E63G (GAA→GGA)  yhhH → putative NTF2 fold immunity protein for polymorphic toxin RhsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,623,5850AG98.3% 194.1 / NA 59E63G (GAA→GGA) yhhHputative NTF2 fold immunity protein for polymorphic toxin RhsB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (30/28);  minor base T (0/1);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTA  >  NC_000913/3623439‑3623728
                                                                                                                                                  |                                                                                                                                               
tAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGAt                                                                                                                                               <  2:1048336/149‑1 (MQ=255)
    aaCTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGc                                                                                                                                           <  2:1222740/149‑1 (MQ=255)
             aTTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGatata                                                                                                                                  >  2:873198/1‑149 (MQ=255)
                     tttgtttgTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGa                                                                                                                          >  1:768600/1‑149 (MQ=255)
                      ttgtttgtTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGaa                                                                                                                         <  1:591421/149‑1 (MQ=255)
                         tttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTCGCTGAGATATATGTTAAGAATCg                                                                                                                      >  2:1282757/1‑149 (MQ=255)
                         tttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCg                                                                                                                      >  1:749284/1‑149 (MQ=255)
                              ttttGCATACAACAGATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCACGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATTTGTTAAGAATCGTTATg                                                                                                                 <  2:1031570/149‑1 (MQ=255)
                               tttACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAACAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATgg                                                                                                                >  1:1124607/1‑149 (MQ=255)
                                ttACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGa                                                                                                               <  2:402053/149‑1 (MQ=255)
                                ttACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGa                                                                                                               >  1:381985/1‑149 (MQ=255)
                                   catacaAAACATGGGCGAAGTGGAATAATAATATTAAAATAATGCGCACGTATGAAAGTGAATGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGACTCGTTATGGACAg                                                                                                            <  2:849100/149‑1 (MQ=255)
                                    atacaACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTGAGAATCGTTATGGACAgg                                                                                                           >  1:494520/1‑149 (MQ=255)
                                        aacaCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGc                                                                                                       <  2:1200729/149‑1 (MQ=255)
                                        aacaCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGAGTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGc                                                                                                       >  2:170520/1‑149 (MQ=255)
                                          cacaTGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATCGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCg                                                                                                     <  1:806801/149‑1 (MQ=255)
                                            caTGGGCGCAGTGTAATAATAATATTAAAATACTGCGCCAGTATGAAAGGGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCa                                                                                                   <  1:947637/149‑1 (MQ=255)
                                             aTGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCAGTGGGATGAGCTGAGAGATATGTTAAGAATCGTTATGGACAGGATGCCGCag                                                                                                  >  1:368801/1‑149 (MQ=255)
                                                  cgcAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCagaagaa                                                                                             <  2:1150975/149‑1 (MQ=255)
                                                  ccgaGTGTAAGAATAATATTAAAATAATGCGCAAGTAGGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCagaagaa                                                                                             <  1:968532/146‑1 (MQ=255)
                                                          aataataatATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAAcc                                                                                     >  2:1120009/1‑149 (MQ=255)
                                                           ctaaaactatTAAAATAATGCGCAAGTATGAAAGTGAACGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCa                                                                                    <  2:975430/142‑1 (MQ=255)
                                                           ataataataTTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATCCCGGTAGAAATTTGGTTAAAAATAAAGCTAGAGCATTGGTATTAGCTGAAATATATGTTTAGAATCGTTATGGACAGGATTCCGCAGAAGAAGAAAACCCa                                                                                    <  2:326775/149‑1 (MQ=255)
                                                             aataataTTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGCAGAAGAAAAACCATa                                                                                  >  1:126083/1‑149 (MQ=255)
                                                             aataataTTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATa                                                                                  >  1:869569/1‑149 (MQ=255)
                                                                  tatTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCGAAAGAAGAAAAACCATACGaaa                                                                             >  1:1012244/1‑149 (MQ=255)
                                                                   atTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGCGATACATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAAt                                                                            >  1:816751/1‑149 (MQ=255)
                                                                   atTAAAATAATGCGCAAATATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAACAATAAAGCTATAGCATTGGGATTAGCTGAGATATACGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAAt                                                                            >  2:1031272/1‑149 (MQ=255)
                                                                    ttAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATc                                                                           >  2:328474/1‑149 (MQ=255)
                                                                     tAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCa                                                                          >  2:69088/1‑149 (MQ=255)
                                                                      aaaaTAATGCGCAAGTATGAAAGGGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTCTAGCATTGGGAATAGCTGAGAAACATGTTAAGAATCGTTATGGACAGGATGCAGCAGAAGAAGAAAAACCATCAGAAAtgaa                                                                         >  1:455177/1‑146 (MQ=255)
                                                                       aaaTAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACt                                                                        >  1:522609/1‑149 (MQ=255)
                                                                        aataatGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCAccg                                                                       >  2:3582/1‑147 (MQ=255)
                                                                         ataatGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGa                                                                      <  2:1115761/149‑1 (MQ=255)
                                                                          taatGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAg                                                                     <  1:743627/149‑1 (MQ=255)
                                                                               cgcAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTaac                                                                <  2:1006035/149‑1 (MQ=255)
                                                                               cgcAAGTATGAAAGTGAAGCTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTaac                                                                <  1:89838/149‑1 (MQ=255)
                                                                                       tGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATAGGAAATCACTGAGTTAACAACAAGtt                                                        >  2:103149/1‑149 (MQ=255)
                                                                                        gAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTg                                                       <  1:873198/149‑1 (MQ=255)
                                                                                        gAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTg                                                       <  2:762089/149‑1 (MQ=255)
                                                                                         aaaGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACGGAGTTAACAACAAGTTgg                                                      >  2:583307/1‑149 (MQ=255)
                                                                                         aaaGTGAAGGGAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTgg                                                      <  1:1044830/149‑1 (MQ=255)
                                                                                          aaGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTAGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGCTGCCGCAGAAGAAGAAAAACCCTACGAAATCACTGAGTTAACAACAAGTTggg                                                     >  2:736414/1‑149 (MQ=255)
                                                                                                  gTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACGGAGATAACAACAGGTTGGGTTGTTGaa                                             >  1:1106000/1‑149 (MQ=255)
                                                                                                    gaaTATCCTGTTAGAAATTTTGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATTGACAGGATGCCGCAGAATAAGACAAACCATACGAAAGCACTGAGTTAACAACAAGTTGGGTTGTTGAAgg                                           <  2:582247/148‑1 (MQ=255)
                                                                                                          aCCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGATGGGTGGTTGAAGGGACCAt                                     >  1:590008/1‑149 (MQ=255)
                                                                                                                 gAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCa                              <  2:749284/149‑1 (MQ=255)
                                                                                                                 gAAATTTGGTTAAAAATAAAGCTATAGCATTCGGATTAGCTGGTATTCATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCa                              <  2:292875/149‑1 (MQ=255)
                                                                                                                  aaaTTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAg                             <  1:1226575/149‑1 (MQ=255)
                                                                                                                    aTTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGAc                           <  2:1042328/149‑1 (MQ=255)
                                                                                                                          ttAAAAATAAAGCTATAGCATTGAGATTAGTTGAGATATAAGTAAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGCCCAAAtt                     >  2:650831/1‑149 (MQ=255)
                                                                                                                           tAAAAATAAAGCTATAGCATTGGGGTTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTAAGTCAGAGCAAAtgc                    >  1:552991/1‑147 (MQ=255)
                                                                                                                            taaaaTAAAACGATTGCATTGGGATTAGCTGAGATATGTGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAATCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCTCTCAGACCAAATTGc                   <  2:285511/148‑1 (MQ=255)
                                                                                                                                aTAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGAGGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtggt               >  1:1284631/1‑149 (MQ=255)
                                                                                                                                  aaaGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtggtgg             <  2:760338/149‑1 (MQ=255)
                                                                                                                                       tataGCAGTGGGATTAGCTGAGATATATGTTAAGAGTCGTGATGGACAGGAGGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGttt        <  2:494520/149‑1 (MQ=255)
                                                                                                                                          aGCATTGGGATTAGCTGAGATCTATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTtta     <  2:525108/149‑1 (MQ=255)
                                                                                                                                           gCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTttat    >  2:1087976/1‑149 (MQ=255)
                                                                                                                                             aTTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTAACTCAGACAGAAATGCTGGTGGGGTCTttatta  >  1:651226/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTA  >  NC_000913/3623439‑3623728

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