Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,129,226 C→T 13.9% intergenic (+32/+35) PP_5376 → / ← PP_5377 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,129,2260CT13.9% 138.5 / 11.6 57intergenic (+32/+35)PP_5376/PP_5377hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (3/5);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGC  >  NC_002947/6129077‑6129364
                                                                                                                                                     |                                                                                                                                          
cTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGcc                                                                                                                                            >  2:776901/1‑150 (MQ=255)
   gACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCAGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAAAGAAAACGGCGCCTTa                                                                                                                                         >  1:207744/1‑150 (MQ=255)
               cTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGttt                                                                                                                             <  1:224187/150‑1 (MQ=255)
                    aCGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGtt                                                                                                                        <  2:149772/150‑1 (MQ=255)
                    aCGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGtt                                                                                                                        >  1:1048800/1‑150 (MQ=255)
                     cGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCCAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGttt                                                                                                                       <  1:1066139/150‑2 (MQ=38)
                      gTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCCAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGtttg                                                                                                                      >  2:1225724/1‑148 (MQ=38)
                      gTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCCAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGtttg                                                                                                                      <  1:1225724/150‑3 (MQ=38)
                         aTGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGaa                                                                                                                   >  2:657489/1‑150 (MQ=255)
                             cTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTg                                                                                                                            <  1:64788/137‑1 (MQ=255)
                             cTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTg                                                                                                                            >  2:64788/1‑137 (MQ=255)
                              tGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGtt                                                                                                                        <  1:842622/140‑1 (MQ=39)
                              tGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGtt                                                                                                                        >  2:842622/1‑140 (MQ=39)
                                   aaaTGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGtt                                                                                                                        <  2:321406/135‑1 (MQ=255)
                                   aaaTGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGtt                                                                                                                        >  1:321406/1‑135 (MQ=255)
                                   aaaTGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTggg                                                                                                                     <  1:80419/138‑1 (MQ=255)
                                   aaaTGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTggg                                                                                                                     >  2:80419/1‑138 (MQ=255)
                                       gACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGt                                                                                                     <  1:1148942/150‑1 (MQ=255)
                                        aCCTATGAAGACTCGACAGGCAAGCTCNACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTg                                                                                                    >  1:1201651/1‑150 (MQ=255)
                                            aTGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCt                                                                                                <  2:1297801/150‑1 (MQ=255)
                                               aaGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGa                                                                                             <  1:422391/150‑1 (MQ=255)
                                                   cTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATAct                                                                                         >  2:937788/1‑150 (MQ=255)
                                                          ggCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAAc                                                                                  <  1:930252/150‑1 (MQ=255)
                                                                 tCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCa                                                                           >  1:1195952/1‑150 (MQ=255)
                                                                        tctcGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTc                                                                    <  2:1201651/150‑1 (MQ=255)
                                                                         ctcGAATACCGCGCAGCTGGCATCAATTGCCGCGTCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTcc                                                                   >  2:576422/1‑150 (MQ=255)
                                                                          tcGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCAAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCg                                                                  >  1:450094/1‑150 (MQ=255)
                                                                           cGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGc                                                                 >  1:1233020/1‑150 (MQ=255)
                                                                                  cgcgCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgtt                                                          <  1:776901/150‑1 (MQ=255)
                                                                                  cgcgCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgtt                                                          >  2:1233925/1‑150 (MQ=255)
                                                                                  cgcgCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgtt                                                          <  2:1048800/150‑1 (MQ=255)
                                                                                   gcgcAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgttg                                                         >  1:98562/1‑150 (MQ=255)
                                                                                    cgcAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCAAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTga                                                        <  2:450094/150‑1 (MQ=255)
                                                                                     gcAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCtt                                                             <  1:1406453/144‑1 (MQ=255)
                                                                                     gcAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCtt                                                             >  2:1406453/1‑144 (MQ=255)
                                                                                          tGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATTCCGTCCGCTCTTGTTGAGACCGa                                                  <  1:550255/150‑1 (MQ=255)
                                                                                               tCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGAcc                                             <  1:657489/150‑1 (MQ=255)
                                                                                                                     aGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCAt                       <  2:1233020/150‑1 (MQ=255)
                                                                                                                      gCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATc                      <  1:1233925/150‑1 (MQ=255)
                                                                                                                       ctctGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATcc                     >  1:882591/1‑150 (MQ=255)
                                                                                                                        tctGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCa                    >  1:483298/1‑150 (MQ=255)
                                                                                                                          tGAATGCTGTCCAAACGAAACCGGCGCCATAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCCCTTGGCGAGACCGATGACCTGGTATGCATTAGAGACACCATCCACt                  >  2:665854/1‑150 (MQ=255)
                                                                                                                          tGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGAcc                                             <  1:409382/123‑1 (MQ=255)
                                                                                                                          tGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGAcc                                             >  2:409382/1‑123 (MQ=255)
                                                                                                                              tGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgtt                                                          <  2:952797/106‑1 (MQ=255)
                                                                                                                              tGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgtt                                                          >  1:952797/1‑106 (MQ=255)
                                                                                                                                 tGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATAc           <  2:483298/150‑1 (MQ=255)
                                                                                                                                   tCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCt         >  1:637330/1‑150 (MQ=255)
                                                                                                                                     cAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGTCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTgg       >  2:418841/1‑150 (MQ=255)
                                                                                                                                      caaCGCAAACGGCGCTTTTAAGGCGCCGTTTTCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGa                                               >  2:1158745/2‑109 (MQ=37)
                                                                                                                                      caaCGCAAACGGCGCTTTTAAGGCGCCGTTTTCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGa                                               <  1:1158745/108‑1 (MQ=37)
                                                                                                                                      aaaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGTCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGga      <  1:418841/150‑1 (MQ=255)
                                                                                                                                      aaaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGga      >  1:1299554/1‑150 (MQ=255)
                                                                                                                                          gAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGc  >  1:424324/1‑150 (MQ=255)
                                                                                                                                          gAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGc  <  2:1195952/150‑1 (MQ=255)
                                                                                                                                          gAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGc  <  1:937788/150‑1 (MQ=255)
                                                                                                                                                   gcCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTg        >  2:828580/1‑135 (MQ=255)
                                                                                                                                                   gcCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTg        <  1:828580/135‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
CTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGC  >  NC_002947/6129077‑6129364

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