Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,482,227 C→A 10.0% intergenic (+26/+28) ynbD → / ← azoR putative phosphatase inner membrane protein/NADH‑azoreductase, FMN‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,482,2270CA10.0% 194.9 / 11.9 60intergenic (+26/+28)ynbD/azoRputative phosphatase inner membrane protein/NADH‑azoreductase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base C (31/23);  new base A (3/3);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

GTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATAAA  >  NC_000913/1482081‑1482374
                                                                                                                                                  |                                                                                                                                                   
gttatgttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTa                                                                                                                                                   <  3:71328/149‑1 (MQ=255)
   atgttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGAGTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAg                                                                                                                                                >  3:359234/1‑149 (MQ=255)
      ttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTcc                                                                                                                                             >  4:61439/1‑149 (MQ=255)
      ttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTcc                                                                                                                                             >  2:372848/1‑149 (MQ=255)
      ttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTcc                                                                                                                                             >  4:127580/1‑149 (MQ=255)
         cGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCcaca                                                                                                                                          <  3:341353/149‑1 (MQ=255)
             cacTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCaa                                                                                                                                      <  3:245275/149‑1 (MQ=255)
                tGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCCACATctcaat                                                                                                                                   <  1:100554/149‑6 (MQ=39)
                 gTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGAcc                                                                                                                                  <  3:326048/149‑1 (MQ=255)
                         gTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCgg                                                                                                                          <  2:218226/149‑1 (MQ=255)
                               gAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGc                                                                                                                    >  3:73420/1‑149 (MQ=255)
                                aaGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCa                                                                                                                   >  2:1493/1‑149 (MQ=255)
                                aaGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCa                                                                                                                   >  2:50293/1‑149 (MQ=255)
                                aaGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCa                                                                                                                   <  1:129287/149‑1 (MQ=255)
                                        aGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAAt                                                                                                           <  1:420966/149‑1 (MQ=255)
                                           tataTTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCt                                                                                                        <  3:381286/149‑1 (MQ=255)
                                           tataTTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCt                                                                                                        <  1:191461/149‑1 (MQ=255)
                                           tataTTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTCAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCt                                                                                                        <  1:372848/149‑1 (MQ=255)
                                              atTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTc                                                                                                     >  1:73557/1‑149 (MQ=255)
                                               ttCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCCACATc                                                                                                                                        >  2:416311/1‑113 (MQ=39)
                                               ttCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCCACATc                                                                                                                                        <  1:416311/113‑1 (MQ=38)
                                                    gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgc                                                                                               >  4:456034/1‑149 (MQ=255)
                                                       aGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgctgc                                                                                            >  1:211963/1‑149 (MQ=255)
                                                        gACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgctgct                                                                                           >  4:268538/1‑149 (MQ=255)
                                                           gCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCtttt                                                                                        <  3:61439/149‑1 (MQ=255)
                                                             cctcAGATTGTGCTGACAGAAGAACACAAAGTGATGATGAGATTATGGGAAAACAGGTAAGTGGAATGAGATGTGGCCTGAATATCTACAGTCCACACCAAGATCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGTTGCTTTTGc                                                                                      <  2:166496/146‑1 (MQ=255)
                                                                          tGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCt                                                                         >  1:120555/1‑149 (MQ=255)
                                                                           gACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCtt                                                                        >  1:186064/1‑149 (MQ=255)
                                                                            aCAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCcac                                                                                                                                           <  2:428912/81‑1 (MQ=37)
                                                                            aCAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCcac                                                                                                                                           >  1:428912/1‑81 (MQ=37)
                                                                             cAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTc                                                                      <  1:50293/149‑1 (MQ=255)
                                                                              aGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACGGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCg                                                                     >  1:157519/1‑149 (MQ=255)
                                                                              aGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCg                                                                     >  3:232118/1‑149 (MQ=255)
                                                                                          aGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCg                                                         <  3:127580/149‑1 (MQ=255)
                                                                                                  tGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGt                                                                                                      >  2:395181/1‑96 (MQ=255)
                                                                                                  tGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGt                                                                                                      <  1:395181/96‑1 (MQ=255)
                                                                                                  tGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTAt                                                 >  1:265124/1‑149 (MQ=255)
                                                                                                  tGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTAt                                                 >  4:401029/1‑149 (MQ=255)
                                                                                                  tGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTAt                                                 <  2:211963/149‑1 (MQ=255)
                                                                                                       ttATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGAt                                            >  3:370851/1‑149 (MQ=255)
                                                                                                            ggAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCtt                                       >  1:254100/1‑149 (MQ=255)
                                                                                                               aaaCAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCgg                                    >  3:97534/1‑149 (MQ=255)
                                                                                                                aaCAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGc                                   >  2:200975/1‑149 (MQ=255)
                                                                                                                   aGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGaa                                <  2:426041/149‑1 (MQ=255)
                                                                                                                         gTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACaaa                          >  4:262511/1‑149 (MQ=255)
                                                                                                                             aTTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTc                      >  4:359177/1‑149 (MQ=255)
                                                                                                                                 agaTGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCg                                     <  2:42008/130‑1 (MQ=255)
                                                                                                                                 agaTGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCg                                     >  1:42008/1‑130 (MQ=255)
                                                                                                                                   aTGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCg                <  4:119028/149‑1 (MQ=255)
                                                                                                                                     gtgGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGt              <  3:268538/149‑1 (MQ=255)
                                                                                                                                     gtgGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGt              >  3:103390/1‑149 (MQ=255)
                                                                                                                                       ggACTGAATATCTACAGTCCAAATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTaa            >  2:125137/1‑149 (MQ=255)
                                                                                                                                        gACTGTAGATATTCAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAAt           >  1:199546/1‑149 (MQ=255)
                                                                                                                                        gACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAAt           >  3:402014/1‑149 (MQ=255)
                                                                                                                                             tatatCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGCAGACAAATTTCACATCGGTAAGGACGa      >  4:307866/2‑149 (MQ=255)
                                                                                                                                             aaTATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTc                      <  3:427783/133‑1 (MQ=255)
                                                                                                                                             aaTATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTc                      >  4:427783/1‑133 (MQ=255)
                                                                                                                                              atatCTACAGTCCACATCAACACCGTGTCCGGTTATGCAGAAACAATGCTGTCAATGGCTGCTTTTGCGTCATACTGTGCTTTCGCTGCCATTTCCGCACCGTATGCGACCCCTTCGGCGAAGCCATCTTTCACATCGGTTATGCCGAt     <  4:359234/149‑1 (MQ=255)
                                                                                                                                                atCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGCATTACCGAtgt   <  2:186064/149‑3 (MQ=255)
                                                                                                                                                 tCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACGGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATaaa  >  4:389524/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATAAA  >  NC_000913/1482081‑1482374

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: