Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,359,141 T→C 7.5% intergenic (‑9/+73) PP_1183 ← / ← PP_1184 siderophore biosynthesis protein/dienelactone hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,359,1410TC 7.5% 150.4 / 2.9 53intergenic (‑9/+73)PP_1183/PP_1184siderophore biosynthesis protein/dienelactone hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (27/22);  new base C (2/2);  total (29/24)
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.52e-01

ATACCGGCGCGCTGGAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAAT  >  NC_002947/1358996‑1359288
                                                                                                                                                     |                                                                                                                                                   
aTACCGGCGCGCTGGAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATtgt                                                                                                                                                  <  4:48159/149‑1 (MQ=255)
 tACCGGCGCGCTGGAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATtgtg                                                                                                                                                 >  4:289587/1‑149 (MQ=255)
    cGGCGCGCTGGAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGc                                                                                                                                                >  4:287663/1‑147 (MQ=255)
     ggCGCGCTGGAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTACTGTGCCAg                                                                                                                                             >  5:119958/1‑149 (MQ=255)
     ggCGCGCTGGAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTACTGTGCCAg                                                                                                                                             >  5:114744/1‑149 (MQ=255)
     gccgcgcTGGAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAg                                                                                                                                             <  8:187867/147‑1 (MQ=255)
      gcgcgcTGGAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGc                                                                                                                                            >  5:364846/1‑149 (MQ=255)
             ggAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCt                                                                                                                                           >  8:121000/1‑143 (MQ=255)
             ggAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCt                                                                                                                                           <  7:121000/143‑1 (MQ=255)
                 gtcgtcGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGTGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGgcgc                                                                                                                                 >  8:140476/1‑149 (MQ=255)
                  tcgtcgGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGTGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGgcgc                                                                                                                                 <  7:140476/148‑1 (MQ=255)
                   cgtcgGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCt                                                                                                                                >  8:326115/1‑148 (MQ=255)
                           gCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGt                                                                                                                       >  1:293039/1‑149 (MQ=255)
                            ccAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCCAAGGGAGGt                                                                                                                       <  4:158749/148‑1 (MQ=255)
                                     cAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGTTGGCGTGTTac                                                                                                             >  7:294992/1‑148 (MQ=255)
                                      aGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCTGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAAc                                                                                                            <  5:202529/149‑1 (MQ=255)
                                             aaGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGcctcct                                                                                                     >  7:244576/1‑149 (MQ=255)
                                                   cAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTACTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCca                                                                                               <  6:119958/149‑1 (MQ=255)
                                                   cAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTACTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCca                                                                                               <  6:114744/149‑1 (MQ=255)
                                                    aaCTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCcac                                                                                              <  1:63761/149‑1 (MQ=255)
                                                         aCCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGa                                                                                         >  5:263625/1‑149 (MQ=255)
                                                          ccAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGAt                                                                                        >  8:131954/1‑149 (MQ=255)
                                                          ccAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGAt                                                                                        <  1:314740/149‑1 (MQ=255)
                                                           cctgcccGCGGCGGGGAATGGCCGGGGGAGGGGCCCTGTG‑‑GGGCGGTTCGGGGGGGCTGCTGGGGGTACTGTGTTCATGGGCGGCTTTTGTGCATCTGGGGATGGGGTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGAtt                                                                                       <  8:294992/144‑1 (MQ=255)
                                                           cAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGAtt                                                                                       >  3:78588/1‑149 (MQ=255)
                                                           cAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGAtt                                                                                       <  6:175764/149‑1 (MQ=255)
                                                             gCAACGCGGCGGGGAATGGACGGGGCAGGGGCCAGTGCACTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGAtt                                                                                       <  1:124059/149‑1 (MQ=255)
                                                              cACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGca                                                                                    <  3:289587/149‑1 (MQ=255)
                                                                ccGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGGTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGAGTGGTAACGCCTCCTTTCCCACAACGATTGcaca                                                                                  >  8:281124/1‑149 (MQ=255)
                                                                      cGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCttt                                                                            <  7:326115/149‑1 (MQ=255)
                                                                         gCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTAt                                                                         <  8:244576/149‑1 (MQ=255)
                                                                             tGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTg                                                                      >  3:206992/1‑148 (MQ=255)
                                                                               gACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGaa                                                                    <  7:64568/148‑1 (MQ=255)
                                                                                aCGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAATTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAaca                                                                  >  4:45879/1‑149 (MQ=255)
                                                                                     gCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGAGTGCACAATCTTTTATTTGACCACcttc                                                             >  2:76769/1‑149 (MQ=255)
                                                                                      cAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACcttct                                                            <  6:364846/149‑1 (MQ=255)
                                                                                       aGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACcttctt                                                           >  3:98545/1‑149 (MQ=255)
                                                                                        ggggCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACcttctt                                                           >  3:114328/1‑148 (MQ=255)
                                                                                                     gtgCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCt                                                <  4:114328/149‑1 (MQ=255)
                                                                                                          ggAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAac                                                                   <  3:93072/125‑1 (MQ=255)
                                                                                                          ggAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAac                                                                   >  4:93072/1‑125 (MQ=255)
                                                                                                           gAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCatat                                          >  3:345689/1‑149 (MQ=255)
                                                                                                           gAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCata                                           >  8:299994/1‑148 (MQ=255)
                                                                                                                gggCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGc                                      >  3:216778/1‑148 (MQ=255)
                                                                                                                              ggAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGc                                              <  1:42232/126‑1 (MQ=255)
                                                                                                                              ggAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGc                                              >  2:42232/1‑126 (MQ=255)
                                                                                                                                aaGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCg                     <  4:216778/149‑1 (MQ=255)
                                                                                                                                       ttCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCtt              >  3:10728/1‑149 (MQ=255)
                                                                                                                                              gcggcTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTtata        <  4:10728/149‑1 (MQ=255)
                                                                                                                                              gcggcTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTtata        <  5:152914/149‑1 (MQ=255)
                                                                                                                                                ggcTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGc      >  7:22676/1‑149 (MQ=255)
                                                                                                                                                    aTTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAAt  >  1:293186/1‑149 (MQ=255)
                                                                                                                                                    aTTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAAt  >  8:312933/1‑149 (MQ=255)
                                                                                                                                                     ttGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCaa   <  7:131954/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATACCGGCGCGCTGGAAGTCGTCGGTTGCCAGGTGGGCAGGGTCGAAGGCGCAACTGACCAGCACCGCGCCGGGCAATGGACGGGGCAGGGGCCAGTG‑‑GTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAAT  >  NC_002947/1358996‑1359288

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: