Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,945,198 C→G 9.2% intergenic (+558/+240) PP_5673 → / ← flhB D‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,945,1980CG 9.2% 152.7 / 4.4 53intergenic (+558/+240)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (21/27);  new base G (2/3);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGC  >  NC_002947/4945053‑4945327
                                                                                                                                                  |                                                                                                                                   
tCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGccgcc‑‑cgc                                                                                                                                <  1:127263/149‑1 (MQ=255)
tCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGccgcc‑‑cgc                                                                                                                                >  5:34436/1‑149 (MQ=255)
 ccACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑cgcg                                                                                                                               >  8:151037/1‑149 (MQ=255)
  cACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑cgcgc                                                                                                                              >  8:153267/1‑149 (MQ=255)
    ccGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCgcgg                                                                                                                            <  7:153267/149‑1 (MQ=255)
     cgcgAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCgcggc                                                                                                                           >  7:290233/1‑149 (MQ=255)
        gAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCa                                                                                                                        <  3:1969/149‑1 (MQ=255)
         aGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGATGCGCCGCG‑‑CGGGCGgcgct                                                                                                                        <  4:11078/149‑2 (MQ=255)
             ctctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATc                                                                                                                      >  3:287692/1‑146 (MQ=255)
             ctctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATc                                                                                                                      <  4:287692/146‑1 (MQ=255)
              tctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATcgc                                                                                                                    >  3:82815/1‑147 (MQ=255)
                     gCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAgct                                                                                                           >  7:219709/1‑149 (MQ=255)
                       aaaaCTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCG‑‑CGGGCGgcgctcga                                                                                                                     >  1:288126/1‑133 (MQ=255)
                       aaaaCTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCG‑‑CGGGCGgcgctcga                                                                                                                     <  2:288126/137‑5 (MQ=255)
                       aaaaCTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGa                                                                                                                  <  1:36836/140‑1 (MQ=255)
                       aaaaCTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGa                                                                                                                  >  2:36836/1‑140 (MQ=255)
                            tttGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTcc                                                                                                     >  4:107144/1‑148 (MQ=255)
                              tGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTa                                                                                                   <  3:274871/148‑1 (MQ=255)
                                 gtcccAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTATTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCTCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACAtt                                                                                               <  4:178646/147‑1 (MQ=255)
                                      aGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATc                                                                                                                      >  6:152362/1‑121 (MQ=255)
                                      aGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATc                                                                                                                      <  5:152362/121‑1 (MQ=255)
                                        ggggCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCg                                                                                        <  4:170730/149‑1 (MQ=255)
                                        ggggCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCg                                                                                        <  4:82815/149‑1 (MQ=255)
                                         gggCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCgg                                                                                       <  3:107144/149‑1 (MQ=255)
                                         gggCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCgg                                                                                       <  7:164515/149‑1 (MQ=255)
                                                    gCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTc                                                                            <  6:34436/149‑1 (MQ=255)
                                                         cATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGGCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCg                                                                       <  8:290233/149‑1 (MQ=255)
                                                         cATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGGGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCg                                                                       <  6:258606/149‑1 (MQ=255)
                                                             gCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGAt                                                                   >  8:178816/1‑149 (MQ=255)
                                                               tAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCTGTTATTCCTGCGGGAtt                                                                  >  3:64957/1‑148 (MQ=255)
                                                                 gAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGCGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTAc                                                               >  5:71036/1‑149 (MQ=255)
                                                                 gAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACAttt                                                                                              <  6:266218/118‑1 (MQ=255)
                                                                 gAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACAttt                                                                                              >  5:266218/1‑118 (MQ=255)
                                                                 gAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTAc                                                               >  4:282710/1‑149 (MQ=255)
                                                                 gAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTAc                                                               <  7:151037/149‑1 (MQ=255)
                                                                 gAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTAc                                                               >  1:145555/1‑149 (MQ=255)
                                                                   aGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTAcgc                                                             <  3:282710/149‑1 (MQ=255)
                                                                   aGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTAcgc                                                             <  2:145555/149‑1 (MQ=255)
                                                                        tCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCTTCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTAcgc                                                             <  1:94933/146‑1 (MQ=255)
                                                                           aaCTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGc                                                     >  5:219197/1‑149 (MQ=255)
                                                                                 gtgAGCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTcg                                                <  7:1055/148‑1 (MQ=255)
                                                                                  tgATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTcgc                                               <  8:219709/148‑1 (MQ=255)
                                                                                    aTCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATTGAGCGCCGCC‑‑CGAGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTa                                            >  1:245282/1‑149 (MQ=255)
                                                                                            tgAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCg                                     <  7:178816/148‑1 (MQ=255)
                                                                                                                 aGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAgg               <  6:219197/149‑1 (MQ=255)
                                                                                                                   gTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCg             >  6:89321/1‑149 (MQ=255)
                                                                                                                   gTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACAGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCg             >  5:128307/1‑149 (MQ=255)
                                                                                                                     ttttGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGt           >  3:125268/1‑149 (MQ=255)
                                                                                                                       ttGTCGCCTGTGAGATCGAGCGCCGGC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCTGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTc          <  4:64957/148‑1 (MQ=255)
                                                                                                                          tCGCCTGTGAGATTGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGGGCCGTACAAACAAGGCGGTcgcg       <  2:245282/148‑1 (MQ=255)
                                                                                                                              cTGTTAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCTCGCATGCTCGCTCCTACATTTGTTTCTGGCCAGTTATTCCTGCGGGATTTACGCGAGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCgcgcgcgc  <  6:71036/149‑1 (MQ=25)
                                                                                                                                        gatgcgcCGCG‑‑CGGGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAg                <  3:185091/123‑1 (MQ=255)
                                                                                                                                        gatgcgcCGCG‑‑CGGGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAg                >  4:185091/4‑126 (MQ=255)
                                                                                                                                           gcgcCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATAtt                                >  3:242283/1‑107 (MQ=255)
                                                                                                                                           gcgcCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATAtt                                <  4:242283/107‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                   
TCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCC‑‑CGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGC  >  NC_002947/4945053‑4945327

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