Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,192,754 C→A 9.4% intergenic (+135/‑67) PP_4572 → / → PP_4573 hypothetical protein/AAA family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,192,7540CA 9.4% 136.6 / 2.8 57intergenic (+135/‑67)PP_4572/PP_4573hypothetical protein/AAA family ATPase
Reads supporting (aligned to +/- strand):  ref base C (30/21);  new base A (3/3);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

CTGAGACAGCGTGACGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGC  >  NC_002947/5192616‑5192901
                                                                                                                                            |                                                                                                                                                    
ctgaGACAGCGTGACGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCttatta                                                                                                                                           <  7:160067/149‑1 (MQ=255)
     aCAGCGTGACGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGGCCAATTACCGGGGTAAAACCCCCCCAACAAACTTATTTGgc                                                                                                                                        >  2:311092/1‑149 (MQ=255)
       aGCGTGACGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGATACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTAAAACCTTATTAGGCGCgg                                                                                                                                    >  5:137247/1‑149 (MQ=255)
        gCGTGACGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCggg                                                                                                                                   <  3:247764/149‑1 (MQ=255)
             aCGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCGCCTAATAAGGTTGTAGGGGCGGGTTTAc                                                                                                                              <  1:313327/149‑1 (MQ=25)
                gggCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCcg                                                                                                                           >  8:251713/1‑149 (MQ=255)
                 ggCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCcgc                                                                                                                          <  1:16529/149‑1 (MQ=255)
                  gCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCGCCTAATAAGGTTGTAGGGGCGGGTTTACCcgcg                                                                                                                         <  2:18519/149‑1 (MQ=25)
                    aTTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCgg                                                                                                                                    >  6:2935/1‑136 (MQ=255)
                    aTTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCgg                                                                                                                                    <  5:2935/136‑1 (MQ=255)
                       cAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTTAAACCCCCCCGACAACCGTATTATGCGCGGGGGTACAcgcgcacac                                                                                                                    >  5:35490/1‑144 (MQ=255)
                           gACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTTGGCGCGGGTTTTCCCGGGAGGAGGccc                                                                                                                >  7:262694/1‑149 (MQ=255)
                            aCTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCg                                                                                                               >  4:204988/1‑149 (MQ=255)
                                aaTACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGGCCCGACACCCTTTTTAGGCGCGGGTTTTCCCGCGCAGCGGCCCggc                                                                                                             >  7:268257/1‑145 (MQ=255)
                                    cccGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTGTACCCGCGAAGAGGCCCGCAAAGccc                                                                                                       >  5:275033/1‑149 (MQ=255)
                                    cccGGGCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGTAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACACCCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGccc                                                                                                       <  3:217640/149‑1 (MQ=255)
                                      cGGNATGCGTTTCCAGCCGCCAGGCCCCGCGTTATAGTCGTAAACAAATCCAGATCATGTTCTGCCTTTTCGGGCCTCTTCGCGGTTAAACCCCCCCCTACAACCTTATTCGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTa                                                                                                     <  6:315763/149‑1 (MQ=38)
                                        gTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTTAACCCGCCCCTTCAACCGTATTTGGCGCGGGGTTACCCGCGAAGCGGCCCGCAAAGCCCGaaa                                                                                                   >  5:306685/1‑149 (MQ=255)
                                                   cAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACaaa                                                                                        >  2:246582/1‑149 (MQ=255)
                                                      ccgccAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAAccc                                                                                     <  7:251713/149‑1 (MQ=255)
                                                       cgccAGGCCCAGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGGTTAAACGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAAcccc                                                                                    >  3:169833/1‑149 (MQ=255)
                                                        gccAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAAcccc                                                                                    >  4:146735/1‑148 (MQ=255)
                                                             gCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGaa                                                                              <  4:117854/149‑1 (MQ=255)
                                                              ccccGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGCTAAACCCGCCCCTACCACAATATAAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCACCAGGaac                                                                             >  8:288464/1‑148 (MQ=38)
                                                                ccGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCTTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAAc                                                                                       >  5:78593/1‑137 (MQ=255)
                                                                ccGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCTTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAAc                                                                                       <  6:78593/137‑1 (MQ=255)
                                                                 cgcgTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAAACCGCCCCGACAACCTTATTATGCGCCGGGTTTCCCCCGAAGAGGCCCGCAAACCCCTAAACAGAGGACAAACCCCCATGAAGTTc                                                                          >  5:40688/1‑149 (MQ=255)
                                                                           tCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGaa                                                                       <  4:89320/142‑1 (MQ=255)
                                                                           tCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGaa                                                                       >  3:89320/1‑142 (MQ=255)
                                                                              gAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCcgc                                                              >  2:40027/1‑148 (MQ=255)
                                                                                 aCAAATCACTATCATGTGCTCCCTTTTCGGGCCTCTTCGGGGTTAACCCGCCCCTAAAACCATATTTCGGCGCGGTTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACt                                                          <  6:35490/149‑1 (MQ=255)
                                                                                     aTCCATATAATTTGCTGCCTTTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGt                                                                            <  5:14564/127‑1 (MQ=255)
                                                                                      tCCAGATCATATGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTTGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTc                                                     >  8:261290/1‑149 (MQ=255)
                                                                                         aGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTTAAACCGCCCCTACAACCGTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAAAGGACAAACCCCCAAGAAGTTAGAAGGCACCCGCCACTACGTCGc                                                   >  4:162107/1‑148 (MQ=38)
                                                                                         aGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAg                                                                      >  7:265715/1‑129 (MQ=255)
                                                                                         aGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACATTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAg                                                                      <  8:265715/129‑1 (MQ=255)
                                                                                                     gcCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCGACACCCTTATTAGGCGCGGGTGTACCCGCGAAGAGGCCCGCAACGCCCGAAACAGAGGCCCACCCCCCAGGCAGGTCGAAGGCCCCCGCGACGACGGCGCCACAGAGGAc                                         >  2:277402/1‑146 (MQ=25)
                                                                                                      ccTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCAGCCGCGACTACGTCGCCACAGATGACCTg                                      >  5:288975/1‑148 (MQ=255)
                                                                                                      ccTGTTCGGGCCTCTTCGCGGGTAAAACCGCCCCTAAAACCATATTAGGCGCGGGTTTTACCGCGAAGAGGCCCGCAAAACCCGAAAAAGAGGGCAAACCCCCATGAAGTTCGAAGGCACCCGCGACAAAGTCGCCACACATGCCCTg                                      >  4:122105/1‑148 (MQ=37)
                                                                                                        tGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGaaa                                   >  8:52042/1‑149 (MQ=255)
                                                                                                          tgcGGGCCTCTTCGCGGGTAAACCCGCGCCTAATAAGGTTGTAGGGGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGACACt                                 >  8:210513/3‑149 (MQ=21)
                                                                                                              gggCCTCTTCGCGGGTAAACCCTCCCCTACAACCTTATTCGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTgg                               <  4:148659/148‑1 (MQ=255)
                                                                                                              gggCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTTGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTgg                               <  7:261290/148‑1 (MQ=255)
                                                                                                                  ctctTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGTCAAACCCCCATGAAGTTCGAAGGCACCCGCGACTAAGTCGCCACAGATGCCCTGAAACTGGCGGt                           >  5:199091/1‑148 (MQ=255)
                                                                                                                  ctctTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGt                           >  5:80898/1‑148 (MQ=255)
                                                                                                                        cgcgGGTAAACCGGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTaaa                        <  3:122105/146‑1 (MQ=255)
                                                                                                                        cgcgGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAAc                       <  6:137247/147‑1 (MQ=255)
                                                                                                                          cgGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCgg                   >  7:192554/1‑149 (MQ=255)
                                                                                                                             gTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCCCCCGCGACTACGACGCCACAGCTGACCTGAAACTGGCGGTAAACGCGGCCa                >  6:278534/1‑149 (MQ=255)
                                                                                                                              tAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGTACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCAt               <  1:246582/149‑1 (MQ=255)
                                                                                                                                 ctccGCCCCTACAACCTTATTCGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTAc            <  1:277402/147‑1 (MQ=255)
                                                                                                                                 aCCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAATCAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGcc                 <  7:95135/144‑1 (MQ=255)
                                                                                                                                 aCCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAGCAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGcc                 >  8:95135/1‑144 (MQ=255)
                                                                                                                                   ccGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTAccc          <  2:22319/149‑1 (MQ=255)
                                                                                                                                    cGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCt         <  8:262694/149‑1 (MQ=255)
                                                                                                                                           aCAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGc  >  6:160110/1‑149 (MQ=255)
                                                                                                                                           aCAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGc  >  3:97227/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CTGAGACAGCGTGACGGGGCATTCAGGGACTGAATACCCGGTCTGCGTTTGCAGCCGCCAGGCCCCGCGTTATAGTCGGAAACAAATCCAGATCATGTGCTGCCTGTTCGGGCCTCTTCGCGGGTAAACCCGCCCCTACAACCTTATTAGGCGCGGGTTTACCCGCGAAGAGGCCCGCAAAGCCCTAAACAGAGGACAAACCCCCATGAAGTTCGAAGGCACCCGCGACTACGTCGCCACAGATGACCTGAAACTGGCGGTAAACGCGGCCATTACCCTCGAACGC  >  NC_002947/5192616‑5192901

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: