Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,830,641 C→T Q277* (CAA→TAA)  gutQ → D‑arabinose 5‑phosphate isomerase GutQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,830,6410CT97.1% 96.3 / NA 35Q277* (CAA→TAA) gutQD‑arabinose 5‑phosphate isomerase GutQ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (13/20);  minor base G (0/1);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GCTGAATAAAGTGCATCATCTGATGCGCCGTGACGATGCCATCCCACAGGTGGCGTTAACCGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGCAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGAC  >  NC_000913/2830397‑2830884
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
gctgaaaaaaggccttcaacTGATGCGCCGGGACGATGCCGCCCCACAGGTGGCTTTAACCGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCg                                                                                                                                                                                                                                                <  1:244143/231‑1 (MQ=255)
  ggAATAAAGGGCATCATCTGGTGCGCCGGGACCATGCCAGCCCACAGGTGGCGTTAACCGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCaa                                                                                                                                                                                                                                             <  1:8537/250‑1 (MQ=255)
              atcatcTGATGCGCCGTGCCGATGCCACCCCACAGGTGGCGTTAACCGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATc                                                                                                                                                                                                                                 <  1:264621/251‑1 (MQ=255)
                                    ggccttCCCAGAGGGGGCGTTAACCGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCtgatga                                                                                                                                                                                                           <  1:8773/246‑1 (MQ=255)
                                                    gCGTTAACCGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAAGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGCCAATGCAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCACAATCACTGc                                                                                                                                                                                           >  1:278134/1‑251 (MQ=255)
                                                          aCCGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTACTCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCAcc                                                                                                                                                                                     >  2:26097/1‑251 (MQ=255)
                                                          aCCGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCAcc                                                                                                                                                                                     >  1:13705/1‑251 (MQ=255)
                                                           ccGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCg                                                                                                                                                                                    <  1:68775/251‑1 (MQ=255)
                                                                     gatggatgCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAAt                                                                                                                                                                                                                                                  <  1:137931/179‑1 (MQ=255)
                                                                     gatggatgCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAAt                                                                                                                                                                                                                                                  >  2:137931/1‑179 (MQ=255)
                                                                     ganggatgCGATGCTGGAACTCAGCCGCGCCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGAt                                                                                                                                                                           <  2:259218/250‑1 (MQ=255)
                                                                       tggatgCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGa                                                                                                                                                                         >  1:124104/1‑250 (MQ=255)
                                                                             cgatgcTGGAACTCAGCCGCACCGGTCTGCGGCTGGTGGCGGTATGTGACGCTCACCACCAGGAACAAGGCGTATTTACCGCCGGCGATTTACGTACCTGGATGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGGCCGCGGCCCCATGTTGTAATCGCAAAGTCGCGCCCTCGGCGCCACAGAGATCCTGCTGAAGCGCAAAATAACTGCCGCACCGGTGGCGGCTGAAAACGGc                                                                                                                                                                  >  2:215525/1‑251 (MQ=255)
                                                                                       aCTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGGTGAAAACGGCAAACTCAccc                                                                                                                                                        >  2:169808/1‑250 (MQ=255)
                                                                                                         ggggCTGGTGGCGGGATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGGCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCGCCGGCGCAATAAACCTGc                                                                                                                                         <  2:267326/248‑1 (MQ=255)
                                                                                                           ggCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGAt                                                                                                                                    <  2:267544/251‑1 (MQ=255)
                                                                                                            gggggtggCGGTATGTGACGCTCAACAACAGGGACAAGGCGGCGTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGCAATCGCAAAGTCGCGCCAGCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGGTGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGAtt                                                                                                                                   <  2:112579/248‑1 (MQ=255)
                                                                                                                         tgtgACGCTCAACAACAGGTACAAGGAGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCCGGATTTCTATCAGGCCgg                                                                                                                      >  2:74544/1‑251 (MQ=255)
                                                                                                                               gCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGattat                                                                                                                >  1:215957/1‑251 (MQ=255)
                                                                                                                                            tACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAACCTCACCGGCGCAATAAACCTGCCGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAAt                                                                                                   >  2:19183/1‑251 (MQ=255)
                                                                                                                                                         ttccGACGGCGATTTACGTCGCTGGCCGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTc                                                                                       <  1:265045/248‑1 (MQ=255)
                                                                                                                                                          aCCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGc                                                                                     >  1:121582/1‑251 (MQ=255)
                                                                                                                                                                ggCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCg                                                                               >  1:42961/1‑251 (MQ=255)
                                                                                                                                                                        ccGTCGCTGGTTGGTTGGCGGCGGCGCACGCACCACGACAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCGGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGt                                                                       <  2:246073/250‑1 (MQ=255)
                                                                                                                                                                         cGTCGCTGGCTGGTTGGCGGCGGCGCACTCCCCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAGTCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCGCTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAAACCCAGACGTTTCGCCAGCCGATGCAGGtt                                                                      <  2:82572/251‑1 (MQ=255)
                                                                                                                                                                                              ggcgCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCAt                                                 >  2:263808/1‑251 (MQ=255)
                                                                                                                                                                                               gcgcCCTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCTCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATc                                                <  1:128372/251‑1 (MQ=255)
                                                                                                                                                                                                cgcACTCACCCCGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATcc                                               <  2:220622/251‑1 (MQ=255)
                                                                                                                                                                                                   ccccaccacGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCcgcg                                            <  2:199578/248‑1 (MQ=255)
                                                                                                                                                                                                      caccacGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCcgcgc                                           <  2:281169/249‑1 (MQ=255)
                                                                                                                                                                                                        cccngCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCgcag                                       <  2:261122/247‑1 (MQ=255)
                                                                                                                                                                                                          ccGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGATGTCGGTTTCCAGCATCCGCGCgcaggc                                     <  2:278134/250‑1 (MQ=255)
                                                                                                                                                                                                                  cAATGAAGCGATGACGGTCGGCGGCACCACGTTGTAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAg                              <  1:169808/250‑1 (MQ=255)
                                                                                                                                                                                                                              gACGGTCGGCGGCACCACGTTGTACTCGCAAAGTCGCGCCATCGACGCAAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAGCTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAAGCCCAGACTTTTCGCCAGCCGATGCAGGTTGGCGACGTCTGTGGCCAGCATCCGCGCGCAGGCAGACCAGTTGTGATGAtttt                 <  2:167945/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cacGTTGGAATCGCAAAGCCGCGCCATCGACGCCAAAGAGATCCTGATGAATCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGAc  <  2:215678/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
GCTGAATAAAGTGCATCATCTGATGCGCCGTGACGATGCCATCCCACAGGTGGCGTTAACCGCCAGCGTGATGGATGCGATGCTGGAACTCAGCCGCACCGGTCTGGGGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGCAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGAC  >  NC_000913/2830397‑2830884

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