Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,045,126 C→A 7.0% intergenic (+25/+32) bglA → / ← ygfF 6‑phospho‑beta‑glucosidase A/putative NAD(P)‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,045,1260CA 7.0% 427.5 / 15.2 122intergenic (+25/+32)bglA/ygfF6‑phospho‑beta‑glucosidase A/putative NAD(P)‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (60/53);  new base A (5/4);  total (65/57)
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.79e-01

Only 100 of 123 total aligned reads displayed.

CGCCCATATCGAAGAGATGAAAAAAGCGGTGACTTACGATGGCGTGGATCTGATGGGCTACACACCGTGGGGCTGCATCGACTGCGTGTCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGGCTCGCCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTC  >  NC_000913/3044867‑3045396
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                               
cgcCCATATCGAAGAGATGAAAAAAGCGGTGACTTACGATGGCGTGGATCTGATGGGCTACACACCGTGGGGCTGCATCGACTGCGTGTCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTGAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGGGAAGGCTTACGTAGATGAATTCACGGATGCAATACTTGTTGAATcc                                                                                                                                                                                                                                                                      >  2:58456/1‑271 (MQ=255)
            aggaataaaaaagGCGGGGACTTACGATGGGGTGGGTCTGAGGGGGTAAACACCGTGGGGCTGCATCGACTGCGTGTCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGGCGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTcc                                                                                                                                                                                                                                                       <  4:521926/267‑1 (MQ=255)
             gagaTGAAAAAAGCGGTTACTTACGATGGCGTGGATCTGATGGGCTACACACCGTGGGGCGGCATCGACTGCGTGTCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCt                                                                                                                                                                                                                                                     <  2:406196/275‑1 (MQ=255)
              gaaggaagaaAGGGGTGGCTTCGGAGGNNGGGGATGTGTAGGGCTACCACCCGTGGGGCNGCCTCGCCTGCGTGTCTGTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACNACGGTACTGGAGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTccc                                                                                                                                                                                                                                                      <  3:63623/269‑1 (MQ=255)
                               aCTTACGATGGCGTGGATCTGATGGGCTACACACCGTGGGGCTGCATCGACTGCGTGTCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGCATAAACATGACGACGGTACTGGCGATATGTCTCGTTCACGTAAGAAGAGCTTTAACTGGTACAACGAGGTGATTGCCAGCAAAGTCGAGAANGTTTAACTTGATGAATTACCGGATGCAAATATTTTTGTAACCGGTCATAAACGTTCACTGttttttttatcccccc                                                                                                                                                                                                                                    >  3:125265/1‑266 (MQ=255)
                                                         cTACACACCGTGGGGCTGCATCGACTGCGTGTCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGAAAAAACTTCCCGTGACg                                                                                                                                                                                                          >  4:284970/1‑274 (MQ=255)
                                                         cTACACACCGTGGGGCTGCATCGACTGCGTGTCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACAGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACg                                                                                                                                                                                                          >  2:156114/1‑274 (MQ=255)
                                                             cgacgggggggGTGCTTCGAGTGGGTGTTGTTCACCACGCGTGAGTACAGCAAACGCTACGGCTGTATATATGTTAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAaga                                                                                                                                                                                                     <  1:12117/268‑1 (MQ=255)
                                                                   tGGGGCTGCATCGACTGCGTGTCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAAATTCCCGTGACGTAAGAGGCttt                                                                                                                                                                                               >  1:92613/1‑275 (MQ=255)
                                                                   tGGGGCTGCATCGACTGCGTGTCCTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACCGGCGATATGTCGCGTTCACGCAAGAAGACCTGTAACTGGTACAAAGAGCTGATTGCCAGCAACGCCGAGAAGCTTTAAGGCGATGAAGTAACGAATGCAATACCTGTTGCATCGCGACAGACAACTGCCCTGTTTTATTTCCCGCCCGCCAAACGGAAAACACTACACGAGACCTaacggactt                                                                                                                                                                                                >  4:27812/1‑267 (MQ=255)
                                                                          gCATCGACTGCGGGTCGTTCACGACCGGGCAGTACAGCAAACGCAACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATATCGCGTTCACGTAAGAAGTGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACtt                                                                                                                                                                                        <  2:376565/275‑1 (MQ=255)
                                                                                  ggcgttccgtgaaccaccGGGCAGTAAGGCAAACGATACGGCTTTATCTATGTGAATAAAAATGACGACGGTACTGGCGATATGTAGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTAGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTAGTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGcc                                                                                                                                                                                 <  4:410441/262‑1 (MQ=255)
                                                                                     gtgtCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTAtt                                                                                                                                                                                                                                             >  2:387105/1‑211 (MQ=255)
                                                                                      tgtCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTcc                                                                                                                                                                                                                 >  2:53631/1‑238 (MQ=255)
                                                                                         cGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTccc                                                                                                                                                                                                                                         <  4:387104/211‑1 (MQ=255)
                                                                                          gTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTg                                                                                                                                                                                                             <  4:53630/238‑1 (MQ=255)
                                                                                                               cAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCt                                                                                                                                                   <  2:334378/275‑1 (MQ=255)
                                                                                                                    gCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgt                                                                                                                                                                                                                                                                             <  7:778586/148‑1 (MQ=255)
                                                                                                                    gCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgt                                                                                                                                                                                                                                                                             <  7:163987/148‑1 (MQ=255)
                                                                                                                         ggCTTTATCTATGTGAATATACATGACGACGGTACTGGAGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACg                                                                                                                                        <  1:86074/275‑1 (MQ=255)
                                                                                                                         ggCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACg                                                                                                                                        <  1:125265/275‑1 (MQ=255)
                                                                                                                              tATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAAc                                                                                                                                                                                                                                                           >  4:429060/1‑156 (MQ=255)
                                                                                                                               aTCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTc                                                                                                                                                                                                                                                                  <  8:407173/148‑1 (MQ=255)
                                                                                                                                  tATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTcc                                                                                                                                                                                                                                                       <  2:429061/156‑1 (MQ=255)
                                                                                                                                   aTGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATcc                                                                                                                                                                                                                                                                      >  6:1119753/1‑140 (MQ=255)
                                                                                                                                     gtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGGAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCa                                                                                                                                                                                                                                                                         <  7:674360/135‑1 (MQ=255)
                                                                                                                                     gtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGTATGCAATACTTGTTGCa                                                                                                                                                                                                                                                                         >  6:577129/1‑135 (MQ=255)
                                                                                                                                     gtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCa                                                                                                                                                                                                                                                                         <  5:577129/135‑1 (MQ=255)
                                                                                                                                     gtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCCCCCTCCATGGGGa                                                                                                                             >  4:426466/1‑274 (MQ=255)
                                                                                                                                     gtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGCGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCTCGCTGCAGGGGGAt                                                                                                                            >  2:194864/1‑275 (MQ=255)
                                                                                                                                              cATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgttg                                                                                                                                                                                                                                                                           <  6:787340/124‑1 (MQ=255)
                                                                                                                                                ttgcgtcgGTACTGGTGGTGTGGCTCGTTCACGGAAAAAGAGCTTTATCGGGTCAAAAGTGGTGATTGCCAGCAACGGCTAGAAGCTTTAAGTCGATGAATTACCGGATGCAATACTTGTTGCTTCCGGTCAGACAACTTCCCTTTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTaa                                                                                                                 <  4:177113/269‑1 (MQ=255)
                                                                                                                                                tGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttg                                                                                                                                                                                                                                                                              >  7:1088414/1‑119 (MQ=255)
                                                                                                                                                  acgacgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAAcg                                                                                                               >  2:498866/1‑275 (MQ=255)
                                                                                                                                                          aCTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTcc                                                                                                                                                                                                                                                       >  7:565219/1‑132 (MQ=255)
                                                                                                                                                          aCTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTcc                                                                                                                                                                                                                                                       <  8:565219/132‑1 (MQ=255)
                                                                                                                                                            tGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACTGCGAGAAGCTTTAAGTCGATGAAGTACCAGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTCTATTTCCCGCCCGCAAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATTACTTATTAGCCTGATAATCGCCGTCGCGACCTCTTCTTCCTTTCCACCGCTCCTGCTGGGGGTGTTCGTCTTAACCCGAGCGatg                                                                                                      >  4:197755/1‑272 (MQ=255)
                                                                                                                                                            tGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgt                                                                                                                                                                                                                                                                             >  8:834942/1‑108 (MQ=255)
                                                                                                                                                            tGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGa                                                                                                                                                                                                                                                               >  6:725004/1‑122 (MQ=255)
                                                                                                                                                            tGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGa                                                                                                                                                                                                                                                               <  5:725004/122‑1 (MQ=255)
                                                                                                                                                            tGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAAAACttgt                                                                                                                                                                                                                                                                             <  7:834942/108‑1 (MQ=255)
                                                                                                                                                                   atGTCGCGTTCACGTAACAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAATTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTTTTTCCCGCCCGCCAAATCGATAAAACTTCCCCTGACGTACTAGGCTTTATCACTTTGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTCTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGCTCGACGCGTCCAg                                                                                               >  1:36151/1‑274 (MQ=255)
                                                                                                                                                                       cgcgTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGAt                                                                                                                                                                                                                          <  5:1046687/148‑1 (MQ=255)
                                                                                                                                                                       cgcgTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGAt                                                                                                                                                                                                                          >  8:135274/1‑148 (MQ=255)
                                                                                                                                                                           ttCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCa                                                                                                                                                                                                                                                                         >  5:732912/1‑97 (MQ=255)
                                                                                                                                                                           ttCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCa                                                                                                                                                                                                                                                                         <  6:732912/97‑1 (MQ=255)
                                                                                                                                                                                  aagaagAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgt                                                                                                                                                                                                                                                                             >  6:953989/1‑86 (MQ=255)
                                                                                                                                                                                  aagaagAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgt                                                                                                                                                                                                                                                                             <  5:953989/86‑1 (MQ=255)
                                                                                                                                                                                    gaagaGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGGCTCGCCGCCGCTGGCGt                                                                             >  4:205360/1‑275 (MQ=255)
                                                                                                                                                                                      agagCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCaaa                                                                                                                                                                                                                               <  8:1116835/128‑1 (MQ=255)
                                                                                                                                                                                      agagCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCaaa                                                                                                                                                                                                                               >  7:1116835/1‑128 (MQ=255)
                                                                                                                                                                                        agCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGc                                                                                                                                                                                                                                                                          >  5:296427/1‑83 (MQ=255)
                                                                                                                                                                                        agCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGc                                                                                                                                                                                                                                                                          <  6:296427/83‑1 (MQ=255)
                                                                                                                                                                                          cTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGAt                                                                                                           >  2:51806/1‑239 (MQ=255)
                                                                                                                                                                                          cTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCACCCCCGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGGCCCCCCCCCGCTGGCGTTCAtt                                                                        >  4:220469/1‑274 (MQ=255)
                                                                                                                                                                                            ttAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGt                                                                                                                                                                                                                                                                   >  5:264155/1‑86 (MQ=255)
                                                                                                                                                                                            ttAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGt                                                                                                                                                                                                                                                                   <  6:264155/86‑1 (MQ=255)
                                                                                                                                                                                            ttAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCg                                                                                                                                                                                                                            >  8:496777/1‑125 (MQ=255)
                                                                                                                                                                                              aaCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGTCTTTATCCCTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCCCGCTGCATGGGGATGTTCGACTTAACGCGATCGCCCCGTCCAGCCTCGCCCCCGCTGGCTTGCATTTCGGt                                                                   >  4:165375/1‑275 (MQ=255)
                                                                                                                                                                                              aaCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGGCTCGCCGCCGCTGGCGTGCATTTCGGt                                                                   <  2:216993/275‑1 (MQ=255)
                                                                                                                                                                                               aCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGGCTCGCCGCCGCTGGCGTGCATTTCGGt                                                                   <  2:205361/274‑1 (MQ=255)
                                                                                                                                                                                                  ggTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgt                                                                                                                                                                                                                                                                             <  7:684392/70‑1 (MQ=255)
                                                                                                                                                                                                  ggTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgt                                                                                                                                                                                                                                                                             >  8:684392/1‑70 (MQ=255)
                                                                                                                                                                                                    tACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGtt                                                                                                                                                                                                                                                  >  7:330257/1‑95 (MQ=255)
                                                                                                                                                                                                    tACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAAt                                                                                                                                                                                                                              <  6:130526/115‑1 (MQ=255)
                                                                                                                                                                                                    tACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAAt                                                                                                                                                                                                                              >  5:130526/1‑115 (MQ=255)
                                                                                                                                                                                                    tACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGGCTCGCCGCCGCTGGCGTGCATTTCGGTataaa                                                              >  1:118540/1‑274 (MQ=255)
                                                                                                                                                                                                    tACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTATTCTGCCTGTCCACCAAGCTGCATGGGGATGTTCGACTTAACGCGAGCGACGCGTCCAGGCCCNCCGCCNNCGGCGGGCCGTTCGGtggaa                                                               >  3:1194/1‑269 (MQ=255)
                                                                                                                                                                                                        aaGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGGCTCGCCGCCGCTGGCGTGCATTTCGGta                                                                  <  1:1194/266‑1 (MQ=255)
                                                                                                                                                                                                          gagGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATCCGGt                                                                                                                                                                                                                                                                   >  8:119821/1‑72 (MQ=25)
                                                                                                                                                                                                          gagGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATCCGGt                                                                                                                                                                                                                                                                   <  7:119821/72‑1 (MQ=37)
                                                                                                                                                                                                            ggTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCAt                                                                                                                                                                                                                                                                        <  6:829809/65‑1 (MQ=37)
                                                                                                                                                                                                            ggTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCAt                                                                                                                                                                                                                                                                        >  5:829809/1‑65 (MQ=37)
                                                                                                                                                                                                             gTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgt                                                                                                                                                                                                                                                                             <  6:428492/59‑1 (MQ=255)
                                                                                                                                                                                                             gTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACtt                                                                                                                                                                                        >  6:946401/1‑144 (MQ=255)
                                                                                                                                                                                                             gTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACtt                                                                                                                                                                                        <  5:946401/144‑1 (MQ=255)
                                                                                                                                                                                                                  atcactcgtcGGCGAGAATATTTCATTCGNNGAAGTACCGGCTTCCACACATGTTACATCCNGTAAACCAAATTCCCTATTTTATCTCCCGCCCGTGATAGCCATACAACTACCCGTGACGTAAGAGGCTTTGTCANTTAGTAGCCAGACAATGACCTGCGCGGCCTCTTGTGCCTGTCCACCACGGTGCATGGGGATGTTCGACTTAACGCGATCACCGCGTCCAGTCTCGCCGCCGCTGGGCTGCATTGCGGTATAAATAGACCGTGGcggaa                                               <  3:60591/266‑5 (MQ=255)
                                                                                                                                                                                                                      aGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATc                                                                                                                                                                                                                                                                       >  8:356963/1‑56 (MQ=25)
                                                                                                                                                                                                                      aGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATc                                                                                                                                                                                                                                                                       <  8:691873/56‑1 (MQ=25)
                                                                                                                                                                                                                      aGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATc                                                                                                                                                                                                                                                                       <  7:356963/56‑1 (MQ=25)
                                                                                                                                                                                                                      aGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATc                                                                                                                                                                                                                                                                       >  7:691873/1‑56 (MQ=25)
                                                                                                                                                                                                                                          ggaagttgtctgaCCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGTCTCGCCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTgcg                       >  2:252917/13‑275 (MQ=255)
                                                                                                                                                                                                                                           gTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTcc                                                                                                                                                                                                                 >  7:980463/1‑89 (MQ=255)
                                                                                                                                                                                                                                               aTGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTc                                                                                                                                                                                                                                                        >  7:160907/1‑46 (MQ=255)
                                                                                                                                                                                                                                               aTGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTc                                                                                                                                                                                                                                                        <  8:160907/46‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACtt                                                                                                                                                                                                                   >  6:487846/1‑78 (MQ=255)
                                                                                                                                                                                                                                                    gTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACtt                                                                                                                                                                                                                   <  5:487846/78‑1 (MQ=255)
                                                                                                                                                                                                                                                        cGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAagag                                                                                                                                                                                                    <  6:1007255/89‑1 (MQ=255)
                                                                                                                                                                                                                                                        cGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAagag                                                                                                                                                                                                    >  5:1007255/1‑89 (MQ=255)
                                                                                                                                                                                                                                                           aTGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTc                                                                                                                                                                                                                  >  7:213358/1‑72 (MQ=255)
                                                                                                                                                                                                                                                           aTGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTc                                                                                                                                                                                                                  <  8:213358/72‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaTACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTgcg                                                                                                                                                                  >  5:785522/1‑116 (MQ=255)
                                                                                                                                                                                                                                                               aaTACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTgc                                                                                                                                                                   >  2:10105/1‑115 (MQ=255)
                                                                                                                                                                                                                                                               aaTACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGGCTCGCCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCGGCGAGGCCACTTCCGCCGATGGTc  >  4:388568/1‑275 (MQ=255)
                                                                                                                                                                                                                                                               aaTACTTGTTGCATCAGGTCAGACAACTTCCCTGTTTTATTTCCCGCTCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCCCTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCTTCCAGCCGCGCCGCCGCCGGCGTGACTTTCGGTATAATTAACCCCTGGCCCCAACCACTCAACGCGGCACCCCTTGGCCGCCGCTTCCAACGaacgcc  >  4:279281/1‑270 (MQ=255)
                                                                                                                                                                                                                                                                 tACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCAc                                                                                                                                                                                          >  8:104576/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                 tACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATggg                                                                                                                               >  7:1142678/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                  aCTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTg                                                                                                                                                      <  8:706195/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aCTTGTTGCATCCGGACAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTgcgc                                                                                                                                                                 >  5:416259/1‑114 (MQ=255)
                                                                                                                                                                                                                                                                   cTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGa                                                                                                                                                               <  4:10105/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                               
CGCCCATATCGAAGAGATGAAAAAAGCGGTGACTTACGATGGCGTGGATCTGATGGGCTACACACCGTGGGGCTGCATCGACTGCGTGTCGTTCACCACCGGGCAGTACAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGTTCGACTTAACGCGATCGACGCGTCCAGGCTCGCCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTC  >  NC_000913/3044867‑3045396

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