Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 826,689 C→G 8.2% intergenic (‑88/‑124) PP_0710 ← / → ycaC‑I hypothetical protein/putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947826,6890CG 8.2% 166.7 / 5.7 57intergenic (‑88/‑124)PP_0710/ycaC‑Ihypothetical protein/putative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base G (3/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

TTGCTGGCGAACGCCGCGACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGTCTGG  >  NC_002947/826540‑826837
                                                                                                                                                     |                                                                                                                                                    
ttGCTGGCGAACGCCGCGACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATc                                                                                                                                                      <  2:618607/150‑1 (MQ=255)
        gAACGCCGCGACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGa                                                                                                                                              <  2:1540160/150‑1 (MQ=255)
        gAACGCCGCGACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGGTTTTTTGTGATTGTGGTGAAGCTATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGa                                                                                                                                              >  1:1013611/1‑150 (MQ=255)
             ccgcGACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCGGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAAttt                                                                                                                                         <  2:475742/150‑1 (MQ=255)
               gcgACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTtgt                                                                                                                                       >  2:1067681/1‑150 (MQ=255)
                  acgaGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGaa                                                                                                                                    <  2:1337455/150‑1 (MQ=255)
                  acgaGGCCGGCGGCGGGTTGCTGCACTTGAAGGTTCTGGTCCACGGGAGGTGCCTCACGGTTTTTTTGTGAGTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGaa                                                                                                                                    <  2:178842/150‑1 (MQ=255)
                       gTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGAt                                                                                                                               >  1:752826/1‑150 (MQ=255)
                                       cGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTa                                                                                                               <  1:1067681/150‑1 (MQ=255)
                                         cACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACa                                                                                                             <  1:580012/150‑1 (MQ=255)
                                            ttGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATGAGTTCAGAAATTTGTGAACTGATGGGGCACAggg                                                                                                                    >  1:1189035/1‑140 (MQ=255)
                                            ttGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATGAGTTCAGAAATTTGTGAACTGATGGGGCACAggg                                                                                                                    <  2:1189035/140‑1 (MQ=255)
                                            ttGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGc                                                                                                          >  2:939618/1‑150 (MQ=255)
                                              gAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGTTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCaa                                                                                                        >  2:211739/1‑150 (MQ=255)
                                                   gtCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGCTTTTGGTGGCGGTGGGCGAGGCGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGa                                                                                                   <  1:1601466/150‑1 (MQ=255)
                                                   gtCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGa                                                                                                   <  2:752826/150‑1 (MQ=255)
                                                        gTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGAc                                                                                              >  2:280512/1‑150 (MQ=255)
                                                            aCGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGtc                                                                                          >  1:1195560/1‑150 (MQ=255)
                                                             cGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGtct                                                                                         >  2:1320993/1‑150 (MQ=255)
                                                              gggAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGtctc                                                                                        >  1:1601634/1‑150 (MQ=255)
                                                                gAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGtctctc                                                                                      <  1:939618/150‑1 (MQ=255)
                                                                  gTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTccc                                                                                    <  1:63908/150‑1 (MQ=255)
                                                                        tCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTc                                                                              >  1:737618/1‑150 (MQ=255)
                                                                            ggTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGTCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGATATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATc                                                                          >  2:25117/1‑150 (MQ=255)
                                                                             gTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTtgt                                                                                                                                       >  2:451406/1‑88 (MQ=255)
                                                                             gTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTtgt                                                                                                                                       <  1:451406/88‑1 (MQ=255)
                                                                             gTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATcc                                                                         <  2:1601634/150‑1 (MQ=255)
                                                                             gTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATcc                                                                         >  2:564267/1‑150 (MQ=255)
                                                                               ttttttGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCtt                                                                       <  1:280512/150‑1 (MQ=255)
                                                                                      tgATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATc                                                                          <  1:915733/140‑1 (MQ=255)
                                                                                      tgATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATc                                                                          >  2:915733/1‑140 (MQ=255)
                                                                                        aTTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGATGGTCTCTCCCACCCTCCATCTTTGCCTACGCt                                                              >  1:90409/1‑150 (MQ=255)
                                                                                              gtgAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGtctct                                                                                       >  2:894853/1‑119 (MQ=255)
                                                                                              gtgAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGtctct                                                                                       <  1:894853/119‑1 (MQ=255)
                                                                                              gtgAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACt                                                        >  1:1110479/1‑150 (MQ=255)
                                                                                              gtgAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACt                                                        >  2:997290/1‑150 (MQ=255)
                                                                                                   gCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATGAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCa                                                                                                                  >  2:1378775/1‑87 (MQ=255)
                                                                                                   gCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATGAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCa                                                                                                                  <  1:1378775/87‑1 (MQ=255)
                                                                                                   gCGAATCTGGTGGCGGTGAGCGAGGGGAAAAATCGCTGGCGGCAAATATTGAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACACTTATACTCCCaa                                                   >  2:432689/1‑150 (MQ=255)
                                                                                                     gATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAAcc                                                 >  1:436343/1‑150 (MQ=255)
                                                                                                        tCTGGTGGCTGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGc                                              <  1:564267/150‑1 (MQ=255)
                                                                                                            gtggCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGaa                                          <  2:436343/150‑1 (MQ=255)
                                                                                                                 ggTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAg                                     <  2:737618/150‑1 (MQ=255)
                                                                                                                  gTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAgg                                    <  1:1320993/150‑1 (MQ=255)
                                                                                                                           gggAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATcc                                                                         <  1:417874/104‑1 (MQ=255)
                                                                                                                           gggAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATcc                                                                         >  2:417874/1‑104 (MQ=255)
                                                                                                                           gggAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACt                           <  1:211739/150‑1 (MQ=255)
                                                                                                                            ggAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACtg                          >  1:7949/1‑150 (MQ=255)
                                                                                                                             gAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGc                                                               >  1:651928/1‑112 (MQ=255)
                                                                                                                             gAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGc                                                               <  2:651928/112‑1 (MQ=255)
                                                                                                                                  aTCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCt                    <  1:25117/150‑1 (MQ=255)
                                                                                                                                       tGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCaa               >  2:1246958/1‑150 (MQ=255)
                                                                                                                                          cggcAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATAc            <  2:7949/150‑1 (MQ=255)
                                                                                                                                           ggcAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACa           >  1:1057286/1‑150 (MQ=255)
                                                                                                                                            gcAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACaa          <  1:997290/150‑1 (MQ=255)
                                                                                                                                             cAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAg         >  2:433854/1‑150 (MQ=255)
                                                                                                                                                   atCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGTCTg   >  1:542080/1‑150 (MQ=255)
                                                                                                                                                    tCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGTCTgg  >  1:659003/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGCTGGCGAACGCCGCGACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGTCTGG  >  NC_002947/826540‑826837

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