Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,874,886 C→G 7.6% G333R (GGC→CGC)  PP_2528 ← O‑acetylhomoserine (thiol)‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,874,8860CG 7.6% 157.0 / 4.1 53G333R (GGC→CGC) PP_2528O‑acetylhomoserine (thiol)‑lyase
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base G (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

GCGGGGTGACAGGCCAGTGACTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATGCGCAGGGCCAGGGTTTCCAG  >  NC_002947/2874750‑2875033
                                                                                                                                        |                                                                                                                                                    
gCGGGGTGACAGGCCAGTGACTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTg                                                                                                                                          <  4:125266/149‑1 (MQ=255)
 cGGGGTGACAGGCCAGTGACTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTgg                                                                                                                                         <  3:111746/149‑1 (MQ=255)
        aCAGGCCAGTGACTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGt                                                                                                                                   <  1:365945/148‑1 (MQ=255)
           ggCCAGTGACTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCAt                                                                                                                               <  6:161615/149‑1 (MQ=255)
              cAGTGACTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGtt                                                                                                                            >  8:337699/1‑149 (MQ=255)
                gTGACTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCg                                                                                                                          >  8:94999/1‑149 (MQ=255)
                 tGACTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCgg                                                                                                                         >  4:58139/1‑149 (MQ=255)
                   aCTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCgggg                                                                                                                       <  2:58150/149‑1 (MQ=255)
                         cGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCg                                                                                                                 <  6:376121/149‑1 (MQ=255)
                             gCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCa                                                                                                             <  3:365327/149‑1 (MQ=255)
                               cGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGc                                                                                                           >  6:312045/1‑149 (MQ=255)
                               cGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGc                                                                                                           <  3:179660/149‑1 (MQ=255)
                                 aTGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGccc                                                                                                         >  8:71965/1‑149 (MQ=255)
                                            gCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCa                                                                                             >  1:118611/1‑149 (MQ=255)
                                               tACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCAccc                                                                                          <  2:356025/149‑1 (MQ=255)
                                                 ccaccaACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAg                                                                                        >  2:164910/1‑149 (MQ=255)
                                                  caccaACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAgg                                                                                       >  8:211875/1‑149 (MQ=255)
                                                  caccaACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAgg                                                                                       >  3:186044/1‑149 (MQ=255)
                                                   accaACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGc                                                                                      >  8:320883/1‑149 (MQ=255)
                                                   accaACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGc                                                                                      <  8:330767/149‑1 (MQ=255)
                                                   accaACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGc                                                                                      >  6:337195/1‑149 (MQ=255)
                                                      aaCTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCAc                                                                                   <  7:225359/149‑1 (MQ=255)
                                                       aCTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCGCCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACt                                                                                  >  1:82779/1‑149 (MQ=255)
                                                          gCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGc                                                                               >  7:237809/1‑149 (MQ=255)
                                                            aaCGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTc                                                                             >  3:320994/1‑149 (MQ=255)
                                                              cgcgTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGt                                                                           <  2:118611/149‑1 (MQ=255)
                                                               gcgTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTg                                                                          >  1:194223/1‑149 (MQ=255)
                                                                cgTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTgg                                                                         <  5:337195/149‑1 (MQ=255)
                                                                  tCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAAGCCACTTGCTCGTGGGc                                                                       >  8:97888/1‑149 (MQ=255)
                                                                        aaaCGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGTCAGGATCGATGCCGGCTTGCCGCCGGTGTATCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGGCCGTGCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGc                                                                   <  1:263725/148‑1 (MQ=255)
                                                                        aaaCGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCa                                                                  <  4:186044/148‑1 (MQ=255)
                                                                        aaaCGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAg                                                                 <  1:19695/149‑1 (MQ=255)
                                                                            gcgcgcCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTaa                                                             <  7:320883/149‑1 (MQ=255)
                                                                               cgcCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAAGCCACTTGCTCGTGGGCTTGCAGGTAATgg                                                          <  7:97888/149‑1 (MQ=255)
                                                                                gcCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATggg                                                         >  7:74183/1‑149 (MQ=255)
                                                                                   aGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGa                                                      <  4:73550/149‑1 (MQ=255)
                                                                                     ccgccTGGCCACCCTTGATGCCGAACGCCAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGAcc                                                    <  2:84700/149‑1 (MQ=255)
                                                                                         cTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTa                                                <  6:323513/149‑1 (MQ=255)
                                                                                         cTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTa                                                <  7:337699/149‑1 (MQ=255)
                                                                                              cACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGcgcg                                           >  6:227297/1‑149 (MQ=255)
                                                                                                       tGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTa                                  <  7:211875/149‑1 (MQ=255)
                                                                                                                  aGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCAt                       >  4:192888/1‑149 (MQ=255)
                                                                                                                   ggATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATg                      <  7:94999/149‑1 (MQ=255)
                                                                                                                   ggATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATg                      >  1:286764/1‑149 (MQ=255)
                                                                                                                    gATCGATGCCGGCTTGCCGCGGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATgc                     >  2:326732/1‑149 (MQ=255)
                                                                                                                    gATCGATGCCGGCTTGCCGCGGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATgc                     >  4:331263/1‑149 (MQ=255)
                                                                                                                    gATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATgc                     >  1:287161/1‑149 (MQ=255)
                                                                                                                         aTGCCGGCTTGCCGCGGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATGCGCAgg                <  3:331263/149‑1 (MQ=255)
                                                                                                                         aTGCCGGCTTGCCGCGGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATGCGCAgg                <  1:326732/149‑1 (MQ=255)
                                                                                                                           gCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATg                      <  3:53951/141‑1 (MQ=255)
                                                                                                                           gCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGATCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCCCCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATg                      >  4:53951/1‑141 (MQ=255)
                                                                                                                              ggCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATGCGCAGGGCCAg           <  3:58139/149‑1 (MQ=255)
                                                                                                                               gCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATGCGCAGGGCCAgg          <  8:74183/149‑1 (MQ=255)
                                                                                                                                    ccgccgGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATGCGCGGGGCCAGGGTTTc     <  2:286764/149‑1 (MQ=255)
                                                                                                                                       ccgGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATGCGCAGGGCCAGGGTTTCCAg  <  3:175151/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
GCGGGGTGACAGGCCAGTGACTTGGCGTCGCCGATGTTCACCAGGCGTACCACCAACTGCAACGCGTCGATGAAACGCGCGCCAGCCGCCTGGCCACCCTTGATGCCGAACGACAGGATCGATGCCGGCTTGCCGCCGGTGTAGCGCTGGGCCAGTGCATGTTCGGGGTGGTCGGGCAGCCCG‑GCGAACTTCACCCAGGCCACTTGCTCGTGGGCTTGCAGGTAATGGGCGACCTTTAGCGCGTTCTCGGTATGGCGCTCCATGCGCAGGGCCAGGGTTTCCAG  >  NC_002947/2874750‑2875033

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