Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,334,705 C→T 5.7% intergenic (+24/‑204) ispB → / → sfsB octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,7050CT 5.7% 634.2 / 16.4 193intergenic (+24/‑204)ispB/sfsBoctaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (82/98);  new base T (5/6);  total (88/105)
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.84e-01

Only 100 of 202 total aligned reads displayed.

TGCGATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGG  >  NC_000913/3334435‑3334969
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                        
ggacggggaccaaggCACCCCAGAACAGGCCCATCTGTTCCGTACCGCCATCGAACAGGGTAACTGTCGCCATCTTCCGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcg                                                                                                                                                                                                                                                                       <  4:342541/261‑1 (MQ=255)
     aataacgaacacgacCGGAAGAGACACAGAGGGTCCGTAGCGCAATCGGACAGGGTAACGGTCGCCATCGTCGGGANCCGGTGCGGGAAGCAATGAACGCTTGTGGTTCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGGGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGa                                                                                                                                                                                                                                                                  <  3:163307/261‑1 (MQ=255)
                    cAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACATGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGACTCGCGCACATCGATGTTCAACGCGATCGTTAATGCCCTCCCCTCACCCCTCGCCGCCCGCGGGATGGTTCCAGTAAGCTc                                                                                                                                                                                                                                                  >  4:384669/1‑275 (MQ=255)
                          gggCCCAGATGACACGTACGGCCATCGAACAGGGTAACGGCCGCCATCTTCTGGAACCGGTGCTGGAAGCAAGGAAGGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGATAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAGCCCGCGCAGCGCGCGGGATGGTTCCAGGAAGTTCCATaaa                                                                                                                                                                                                                                            <  4:483120/274‑1 (MQ=255)
                          aGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                                                                                                                                            <  2:622946/275‑1 (MQ=255)
                            gCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGACGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcgc                                                                                                                                                                                                                                                                      >  1:163308/1‑247 (MQ=255)
                                      tCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                                                                                                                                                >  2:483122/1‑259 (MQ=255)
                                       ccggaacGCAATCGAAAAGGGTAACGGGCGCCAACTTATGGAAGCGGTTCTGGAAGCAAGGAACGCTTGGGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAg                                                                                                                                                                                                                               <  4:226500/269‑1 (MQ=255)
                                                          ggTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACCCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGaaa                                                                                                                                                                                                              >  2:167922/1‑273 (MQ=255)
                                                                                  acgggtagggAAGAAATTGACGTGTGTGGGTGTCTTGAATGGACGCGTCTGCGTGGCGAGGAATAAGTAGGCAAAGGCATCGCAGCGTTAGCGGTGACTCCTGGCACGCCATGGCGGGGAGCAGTCATCGGACTCGCGCAAATCGCTGGTCAACGCGATAGTTAATCCCCTCCCCTAATCCCGCGGAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTatatat                                                                                                                                                                                     <  4:112053/266‑1 (MQ=255)
                                                                                    ggTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                                                                                                                                                >  2:640080/1‑213 (MQ=255)
                                                                                        cTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                                                                                                                                            <  4:640078/213‑1 (MQ=255)
                                                                                                    aaCGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGt                                                                                                                                                                                                                                                        <  4:168009/189‑1 (MQ=255)
                                                                                                                   ctTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  2:647980/1‑224 (MQ=255)
                                                                                                                       aaTGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  4:647978/224‑1 (MQ=255)
                                                                                                                       aaTGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTagag                                                                                                                                                <  2:384671/274‑1 (MQ=255)
                                                                                                                               actCAGCTGGCCAAGGAAAGAGCGGCCAAAGGCCTAGCCGCGTTACAGGTGCGCCCGGAAACCCCTGGGCGAGACGCACTCATCGGCCTGGCGCACATGCCGGGTCAACGCGATCGTTAATCCACACACCTCATCACACGCTGCGCGCGGGATGGTTCCAGTAATTTCCATAAAAACAGTTATTCAGCTCTGTTAAACCTGAAAAATAGGGATGTGCATCTTTATAAATTCTGAATATTCACACTCGTTAAAGGAACTTTTTAGAGCAATAGGcc                                                                                                                                       <  4:538669/273‑1 (MQ=255)
                                                                                                                                 tCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCggg                                                                                                                                                                                                                                                                   >  8:626353/1‑149 (MQ=255)
                                                                                                                                  cAGCGTGCCGAGGAGGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCggg                                                                                                                                                                                                                                                                   <  6:148968/148‑1 (MQ=255)
                                                                                                                                  cAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                                                                                                                                              >  4:599981/1‑169 (MQ=255)
                                                                                                                                   aGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGAt                                                                                                                                                                                                                                                                 <  7:626353/149‑1 (MQ=255)
                                                                                                                                     cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGATga                                                                                                                                                                                                                                                               <  8:78118/149‑2 (MQ=255)
                                                                                                                                      gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGt                                                                                                                                                                                                                                                              >  7:87641/1‑149 (MQ=255)
                                                                                                                                           gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCa                                                                                                                                                                                                                                                          <  6:383705/148‑1 (MQ=255)
                                                                                                                                              gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                      <  8:170060/149‑1 (MQ=255)
                                                                                                                                               aagaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                      <  6:91249/148‑1 (MQ=255)
                                                                                                                                                      aGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                                                                                                                                              >  7:581278/1‑149 (MQ=255)
                                                                                                                                                       gACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                                                                                                                                                     >  2:588857/1‑189 (MQ=255)
                                                                                                                                                          aaaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                                                                                                                                                                           <  5:773637/148‑1 (MQ=255)
                                                                                                                                                           aaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTcc                                                                                                                                                                                                                                                           <  5:429463/131‑1 (MQ=255)
                                                                                                                                                           aaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaaa                                                                                                                                                                                                                                          >  5:34227/1‑148 (MQ=255)
                                                                                                                                                            aGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAaca                                                                                                                                                                                                                                        <  7:903165/149‑1 (MQ=255)
                                                                                                                                                                aTCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTa                                                                                                                                                                                                                                    >  8:490233/1‑149 (MQ=255)
                                                                                                                                                                 tCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                                   >  8:440617/1‑149 (MQ=255)
                                                                                                                                                                  cGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAtt                                                                                                                                                                                                                                  <  7:568193/149‑1 (MQ=255)
                                                                                                                                                                  cGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCta                                                                                                                                                                                                                          >  1:33377/1‑157 (MQ=255)
                                                                                                                                                                     aGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAg                                                                                                                                                                                                                               <  7:274025/149‑1 (MQ=255)
                                                                                                                                                                      gCGTTACAGGTGCTCCCGNACACCCNTTGNCGAGAAGNACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAac                                                                                                                                                                                                                      <  3:33376/157‑1 (MQ=255)
                                                                                                                                                                       cGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGc                                                                                                                                                                                                                              >  6:490932/1‑148 (MQ=255)
                                                                                                                                                                        gTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcgc                                                                                                                                                                                                                                                                      >  2:656931/1‑107 (MQ=255)
                                                                                                                                                                         ttACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctc                                                                                                                                                                                                                            <  6:64788/148‑1 (MQ=255)
                                                                                                                                                                          tACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAtt                                                                                                                                                                                                                                  <  2:84978/141‑1 (MQ=255)
                                                                                                                                                                            cAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGa                                                                                                                                                                                                                                                                  <  4:656929/107‑1 (MQ=255)
                                                                                                                                                                             aGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                                                                                                                                       >  7:195722/1‑149 (MQ=255)
                                                                                                                                                                              ggTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                                                                                                                                       <  6:34227/148‑1 (MQ=255)
                                                                                                                                                                                    tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCt                                                                                                                                                                                                                  >  5:557603/1‑148 (MQ=255)
                                                                                                                                                                                    tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                                                                                                                                                 >  8:42784/1‑149 (MQ=255)
                                                                                                                                                                                      ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATg                                                                                                                                                                                                                                                                <  6:676817/100‑1 (MQ=255)
                                                                                                                                                                                            acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTc                                                                                                                                                                                                                                                            <  8:582848/98‑1 (MQ=255)
                                                                                                                                                                                            acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                                                                                                                                         >  7:798121/1‑149 (MQ=255)
                                                                                                                                                                                             ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGATg                                                                                                                                                                                                                                                                <  5:690784/93‑1 (MQ=255)
                                                                                                                                                                                             ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                                                                                                                                         <  5:13961/148‑1 (MQ=255)
                                                                                                                                                                                               ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      <  8:195722/149‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      <  6:557603/148‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                                                                                                                                                     <  8:83828/149‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCt                                                                                                                                                                                             >  1:1501/1‑157 (MQ=255)
                                                                                                                                                                                                  tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgt                                                                                                                                                                                                   <  7:832679/149‑1 (MQ=255)
                                                                                                                                                                                                    gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcg                                                                                                                                                                                                                                                                       <  2:118825/79‑1 (MQ=255)
                                                                                                                                                                                                    gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  4:667252/1‑144 (MQ=255)
                                                                                                                                                                                                        gaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  2:667254/144‑1 (MQ=255)
                                                                                                                                                                                                          aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAACCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCtt                                                                                                                                                                                            <  6:651944/148‑1 (MQ=255)
                                                                                                                                                                                                              cTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  2:546492/1‑134 (MQ=255)
                                                                                                                                                                                                               tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  >  6:618439/1‑137 (MQ=255)
                                                                                                                                                                                                               tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                                                                                                                                                      >  7:403580/1‑149 (MQ=255)
                                                                                                                                                                                                                cATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                                                                                                                                                      >  5:295683/1‑148 (MQ=255)
                                                                                                                                                                                                                 aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                    >  8:531183/1‑149 (MQ=255)
                                                                                                                                                                                                                  tCGGCGTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                                                                                                                                                   <  7:149089/149‑1 (MQ=255)
                                                                                                                                                                                                                  tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                    >  6:168541/1‑148 (MQ=255)
                                                                                                                                                                                                                  tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                                                                                                                                                   <  8:798121/149‑1 (MQ=255)
                                                                                                                                                                                                                    ggCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTg                                                                                                                                                                                 <  8:87641/149‑1 (MQ=255)
                                                                                                                                                                                                                      ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGTCTGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTTCATCTTTTNATATTCTGAATATTCACAATCTTAACCGGAACTTTTTAGAGTAATAGGCCATCTGGAGAATAGTGATTCTCGAANGAAGCAGTGTCTTGAATGAAAGCAAATAAANNAAGAAGGAAGGAAAAAGGGtaaaacattt                                                 >  3:88229/1‑268 (MQ=255)
                                                                                                                                                                                                                      ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                                                                                                                                               <  7:490233/149‑1 (MQ=255)
                                                                                                                                                                                                                       cTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                                                                                                                                               <  5:889742/148‑1 (MQ=255)
                                                                                                                                                                                                                           cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  4:743242/1‑121 (MQ=255)
                                                                                                                                                                                                                               caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcg                                                                                                                                                                                                                                                                       <  4:29526/52‑1 (MQ=255)
                                                                                                                                                                                                                               caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  2:743244/121‑1 (MQ=255)
                                                                                                                                                                                                                                 tCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAACTTCCATAAAAACACTTATTGAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCGCtctc                                                                                                                                                                    <  8:573657/149‑1 (MQ=255)
                                                                                                                                                                                                                                  cGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGAGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCa                                     <  1:88230/275‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGct                                                                                                                                                                                                                             >  4:447259/1‑89 (MQ=255)
                                                                                                                                                                                                                                        tCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtat                                                                                                                                                                                                                         <  2:447261/89‑1 (MQ=255)
                                                                                                                                                                                                                                        tCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAgg                                                                                                                                         >  4:800805/1‑169 (MQ=255)
                                                                                                                                                                                                                                           aCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGaa                                                                                                                                                          >  7:730291/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCggg                                                                                                                                                                                                                                                                   <  6:284221/44‑1 (MQ=18)
                                                                                                                                                                                                                                            cgcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCa                                                                                                                                      >  2:270150/1‑168 (MQ=255)
                                                                                                                                                                                                                                             gcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTa                                                                                                                                                                >  5:739974/1‑141 (MQ=255)
                                                                                                                                                                                                                                                aTCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAg                                                                                                                                  <  4:270149/168‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                   <  8:635756/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                   <  8:403580/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        tcccctccccTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCaa                                                                                                                                             >  8:279176/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          ccctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAAt                                                                                                                                            >  5:616401/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           cctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                                                                                                                                           >  5:668523/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcac                                                                                                                                                                        <  5:649137/118‑1 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg                                                                                                                                          >  6:820800/1‑148 (MQ=255)
                                                                                                                                                                                                                                                              cccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGcc                                                                                                                                       <  7:42784/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaag      >  2:214760/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   cATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa                                                                                                                                   >  6:737285/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                    aTCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                                                                                                                                                    >  5:21542/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                     tCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCAATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAAATTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCCTCCCGCAGACATCATTGCGGGttttggacgtaagg  >  4:582359/1‑264 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGAt                                                                                                                      >  2:787312/1‑155 (MQ=255)
                                                                                                                                                                                                                                                                          cgcAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGtat                                                                                                                            <  5:490932/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                        
TGCGATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGG  >  NC_000913/3334435‑3334969

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