Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 924,071 (A)5→6 intergenic (‑51/‑228) U712_04490 ← / → U712_04495 Vegetative catalase/Aliphatic sulfonates import ATP‑binding protein SsuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881924,0661.A98.2% 192.0 / ‑3.0 54intergenic (‑46/‑233)U712_04490/U712_04495Vegetative catalase/Aliphatic sulfonates import ATP‑binding protein SsuB
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (20/33);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAAT  >  CP006881/923932‑924196
                                                                                                                                       |                                                                                                                                   
aaaGTTGGTCCGCGAGAACCGGCAGTCATTGAGCTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGTAATTTATAAGAACATAATATCATCTCTAAGTaaa                                                                                                                                 <  6:335793/138‑1 (MQ=255)
 aaGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaa                                                                                                                                >  8:471598/1‑138 (MQ=255)
    ttGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaaaagt                                                                                                                            <  5:137279/139‑4 (MQ=255)
    ttGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaaaagt                                                                                                                            <  1:293051/139‑4 (MQ=255)
    ttGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaaaagt                                                                                                                            <  7:357135/139‑4 (MQ=255)
    ttGGTCCGCGAGAACCGGCAGTCATTGAGTTTAGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaaaagt                                                                                                                            <  3:373987/139‑4 (MQ=255)
      ggTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaaaagtc                                                                                                                           >  4:137930/1‑134 (MQ=255)
      ggTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaaaagtc                                                                                                                           >  6:80068/1‑134 (MQ=255)
      ggTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaaaagtc                                                                                                                           >  7:462648/1‑134 (MQ=255)
      ggTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaaaagtc                                                                                                                           >  7:456734/1‑134 (MQ=255)
         ccGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAata                                                                                                                        <  8:405719/139‑1 (MQ=255)
         ccGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAata                                                                                                                        <  4:374538/139‑1 (MQ=255)
          cgcgAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAata                                                                                                                        <  1:467830/138‑1 (MQ=255)
          cgcgAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAata                                                                                                                        <  1:187984/138‑1 (MQ=255)
                     gCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAat                                                                                                            >  8:273593/1‑139 (MQ=255)
                     gCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAat                                                                                                            <  8:114211/139‑1 (MQ=255)
                     gCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAat                                                                                                            <  4:408771/139‑1 (MQ=255)
                     gCAGTCAATGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTTGAAGTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAat                                                                                                            <  1:249535/139‑1 (MQ=255)
                       aGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAatta                                                                                                          <  5:223960/139‑1 (MQ=255)
                           aTTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAAttatta                                                                                                       >  6:270947/1‑138 (MQ=255)
                            ttGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATtata                                                                                                     <  5:270947/139‑1 (MQ=255)
                            ttGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATtata                                                                                                     <  4:145243/139‑1 (MQ=255)
                            ttGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATtata                                                                                                     <  8:361307/139‑1 (MQ=255)
                                   ttGATTATCTCCAACCGGAGCGCGCCAGCTAGTTGTCAGTTTATATGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAg                                                                                               <  3:81969/138‑1 (MQ=255)
                                      aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAt                                                                                           <  4:190827/139‑1 (MQ=255)
                                      aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAt                                                                                           <  7:471598/139‑1 (MQ=255)
                                      aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAt                                                                                           <  3:137930/139‑1 (MQ=255)
                                      aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAt                                                                                           <  6:437546/139‑1 (MQ=255)
                                       ttATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAtc                                                                                          <  7:273593/139‑1 (MQ=255)
                                           ctcCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTa                                                                                      >  3:471537/1‑139 (MQ=255)
                                              cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTAtt                                                                                    <  8:462648/138‑1 (MQ=255)
                                               aaCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGTACTCATGTTATCACCTCTTTGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa                                                                                  <  3:132772/139‑1 (MQ=255)
                                               aaCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa                                                                                  <  3:457302/139‑1 (MQ=255)
                                               aaCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa                                                                                  <  2:161686/139‑1 (MQ=255)
                                                aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt                                                                                 <  7:200035/139‑1 (MQ=255)
                                                aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt                                                                                 <  1:141533/139‑1 (MQ=255)
                                                aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt                                                                                 <  8:456734/139‑1 (MQ=255)
                                                aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt                                                                                 <  5:349625/139‑1 (MQ=255)
                                                  cGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCg                                                                               <  4:471537/139‑1 (MQ=255)
                                                     aGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCg                                                                               <  3:155236/136‑1 (MQ=255)
                                                     aGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCg                                                                               >  4:155236/1‑136 (MQ=255)
                                                      gcgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGata                                                                            >  7:334195/1‑138 (MQ=255)
                                                        gcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGatata                                                                          >  7:337301/1‑138 (MQ=255)
                                                           ccAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAt                                                                                           >  6:33399/1‑118 (MQ=255)
                                                           ccAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAt                                                                                           <  5:33399/118‑1 (MQ=255)
                                                                   ttGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAAt                                                              <  8:31837/139‑1 (MQ=255)
                                                                   ttGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAAt                                                              <  2:228694/139‑1 (MQ=255)
                                                                   ttGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAAt                                                              <  8:337301/139‑1 (MQ=255)
                                                                   ttGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAAt                                                              <  5:80068/139‑1 (MQ=255)
                                                                   ttGTCAGTTTATTTGAACTCATGTTATCACCGCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGatat                                                              >  3:8813/1‑139 (MQ=255)
                                                                                            aTCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAg                                      <  2:284361/138‑1 (MQ=255)
                                                                                               aCCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTATAAAAGTCTATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAAGCTTGTTGTAACAGCttt                                  >  8:143770/1‑139 (MQ=255)
                                                                                                 ctctTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCttttt                                >  5:75200/1‑139 (MQ=255)
                                                                                                 ctctTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCtttt                                 >  3:242301/1‑138 (MQ=255)
                                                                                                   ctTGGAATTTATAAGAACATAATATCATCTCTAATTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTAt                              >  2:112907/1‑139 (MQ=255)
                                                                                                   ctTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTAt                              >  8:62677/1‑139 (MQ=255)
                                                                                                      ggAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGac                           >  8:481377/1‑139 (MQ=255)
                                                                                                                       aataTCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTg            >  1:101309/1‑137 (MQ=255)
                                                                                                                            catcTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGaaa      >  2:111825/1‑138 (MQ=255)
                                                                                                                            catcTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGaaa      >  1:5229/1‑138 (MQ=255)
                                                                                                                             atcTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaat     >  4:413736/1‑138 (MQ=255)
                                                                                                                              tctcTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataa   >  3:252922/1‑139 (MQ=255)
                                                                                                                              tctcTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataa   >  1:81622/1‑139 (MQ=255)
                                                                                                                               ctctAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat  >  6:159474/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                   
AAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAAT  >  CP006881/923932‑924196

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