Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 6,013,885 T→C 100% intergenic (+90/+20) xpt → / ← PP_5266 Xanthine phosphoribosyltransferase/putative acetyl‑CoA hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154516,013,8850TC100.0% 157.6 / NA 43intergenic (+90/+20)xpt/PP_5266Xanthine phosphoribosyltransferase/putative acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/27);  total (16/27)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCT  >  AE015451/6013656‑6014110
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
ccACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAt                                                                                                                                                                                                                             >  1:850437/1‑236 (MQ=255)
    cGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGcc                                                                                                                                                                                                                         <  1:1080100/236‑1 (MQ=255)
                           tCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGcc                                                                                                                                                                                                                         <  1:615399/213‑1 (MQ=255)
                           tCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgccg                                                                                                                                                                                                  >  1:546176/1‑236 (MQ=255)
                              agaAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATcctgtctcttat                                                                                                                                                                                                >  1:927551/1‑224 (MQ=255)
                                      ttCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGcc                                                                                                                                                                                                                         <  2:927546/202‑1 (MQ=255)
                                                   gTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAg                                                                                                                                                                              >  1:1209829/1‑232 (MQ=255)
                                                              gaaaGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgccgcc                                                                                                                                                                                                <  2:576213/200‑1 (MQ=255)
                                                                           aCCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGcccc                                                                                                                                                   >  1:609621/1‑235 (MQ=255)
                                                                                                                  tgggtgtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAgg                                                                                                                     <  1:75334/224‑1 (MQ=255)
                                                                                                                  ggggtgtggtCAGCTTCATCGAATAAACGCGACACCTGCGTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaac                                                                                                                                                                                 <  2:8645/164‑1 (MQ=255)
                                                                                                                             aGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGccc                                                                                                                         <  2:609616/211‑1 (MQ=255)
                                                                                                                             aGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGCACCTGTGCGAGCGCGCCTGCCCGCGAATGGcc                                                                                                                                                                                                                         >  1:318924/1‑115 (MQ=37)
                                                                                                                                                     cTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTg                                                                        <  1:523397/236‑1 (MQ=255)
                                                                                                                                                              gggCACGCCCGCTCCCACAGGTGTAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCt                                                                 <  1:508821/234‑1 (MQ=255)
                                                                                                                                                              gggCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGTAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCTCGCCAACAACAAGCGCTGGTaacgct                                                                                                                                                                 >  2:271801/1‑134 (MQ=255)
                                                                                                                                                                             cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCt                                                    >  2:908316/1‑232 (MQ=255)
                                                                                                                                                                                aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTc                                                                                                             <  1:400622/172‑1 (MQ=255)
                                                                                                                                                                                       gcaCAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCCGCACCTCGCTCAGCTTgca                                           >  2:984451/1‑229 (MQ=255)
                                                                                                                                                                                         acaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGct                                     >  1:668018/1‑235 (MQ=255)
                                                                                                                                                                                               ctcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGc                                                            >  1:66048/1‑206 (MQ=255)
                                                                                                                                                                                                tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaacaa                                                                                                                                                                               <  1:867556/90‑1 (MQ=255)
                                                                                                                                                                                                tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaacaa                                                                                                                                                                               >  2:867551/1‑90 (MQ=255)
                                                                                                                                                                                                tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTg                                                                        <  1:15315/193‑1 (MQ=255)
                                                                                                                                                                                                tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTg                                                                        >  2:15314/1‑193 (MQ=255)
                                                                                                                                                                                                     cGCTGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTc                                                       <  1:918377/205‑1 (MQ=255)
                                                                                                                                                                                                        tgtgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAg                      <  1:1146980/235‑1 (MQ=255)
                                                                                                                                                                                                           gCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGcc                                                                                               >  1:1189800/1‑159 (MQ=255)
                                                                                                                                                                                                           gCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGcc                                                                                               <  2:1189795/159‑1 (MQ=255)
                                                                                                                                                                                                            ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTg                                                                        <  2:850432/181‑1 (MQ=255)
                                                                                                                                                                                                            ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTc                                                       <  2:546171/198‑1 (MQ=255)
                                                                                                                                                                                                                 ccTGTGGGGGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGccc                                                                                              <  2:109662/154‑1 (MQ=255)
                                                                                                                                                                                                                 ccTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGc                                                                                                <  2:623913/152‑1 (MQ=255)
                                                                                                                                                                                                                  cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaacaa                                                                                                                                                                               >  2:1279092/1‑72 (MQ=255)
                                                                                                                                                                                                                  cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaacaa                                                                                                                                                                               <  1:1279097/72‑1 (MQ=255)
                                                                                                                                                                                                                   tgtgGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTg                                                                        <  1:78278/174‑1 (MQ=255)
                                                                                                                                                                                                                     gggggGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgccgcc                                                                                                                                                                                                <  2:251076/51‑1 (MQ=25)
                                                                                                                                                                                                                       ggAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCa                                                                                                                       <  1:133240/123‑1 (MQ=255)
                                                                                                                                                                                                                       ggAGCGGGCGGGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCt                                                                 <  1:153436/177‑1 (MQ=255)
                                                                                                                                                                                                                           cGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGct  <  1:456172/236‑1 (MQ=255)
                                                                                                                                                                                                                            gggCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAg                                                                                                                                                                   <  1:1154260/74‑1 (MQ=255)
                                                                                                                                                                                                                              gCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCt                                                        <  2:550970/179‑1 (MQ=255)
                                                                                                                                                                                                                               cGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCAt                                                                                      >  2:1275377/1‑148 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
CCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCT  >  AE015451/6013656‑6014110

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