Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F41 I0 R1 212 10.4 395026 92.3% 364608 148.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,801 G→C 80.6% A188A (GCG→GCC talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,8010GC80.6% 166.4 / 35.4 67A188A (GCG→GCCtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base G (8/5);  new base C (32/22);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

TATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGG  >  CP009273/8654‑8946
                                                                                                                                                   |                                                                                                                                                 
tATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGaa                                                                                                                                                 <  2:15725/150‑1 (MQ=255)
    cGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATAcc                                                                                                                                             >  1:4973/1‑150 (MQ=255)
         cgctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGac                                                                                                                                          >  2:169550/5‑147 (MQ=255)
              aaCAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATaaga                                                                                                                                      >  2:5890/1‑147 (MQ=255)
                  gCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTa                                                                                                                               <  2:19050/150‑1 (MQ=255)
                      gAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGca                                                                                                                           >  2:80813/1‑149 (MQ=255)
                      gAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGca                                                                                                                           <  1:134773/150‑2 (MQ=255)
                       aaaaaGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGcac                                                                                                                          >  1:52724/1‑148 (MQ=255)
                        aaaaGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCAcc                                                                                                                         >  2:191553/1‑150 (MQ=255)
                         aaaGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCg                                                                                                                        >  1:67635/1‑150 (MQ=255)
                           agaagGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCgg                                                                                                                       >  2:49892/1‑149 (MQ=255)
                               ggCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCag                                                                                                                    >  1:107487/1‑148 (MQ=255)
                                 cATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCagaaga                                                                                                                <  1:191553/150‑1 (MQ=255)
                                   ttaaCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATc                                                                                                              <  2:194250/148‑1 (MQ=255)
                                   tCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATc                                                                                                              <  1:38708/150‑1 (MQ=255)
                                       cTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCggg                                                                                                          >  2:192567/1‑150 (MQ=255)
                                       cTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCggg                                                                                                          <  1:158248/150‑1 (MQ=255)
                                             ccTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGAGTGGTACAAAGCCAACTCCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGt                                                                                                    <  1:179167/150‑1 (MQ=255)
                                              cTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGtt                                                                                                   >  1:168965/1‑150 (MQ=255)
                                              cTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGtt                                                                                                   <  1:8436/150‑1 (MQ=255)
                                               tGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGttt                                                                                                  >  2:62865/1‑150 (MQ=255)
                                               tGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGttt                                                                                                  >  1:172228/1‑150 (MQ=255)
                                                   ccTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGt                                                                                              <  1:101798/150‑1 (MQ=255)
                                                    ctgctgTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGTGCGTGGTTTCCGtc                                                                                             <  2:145374/150‑2 (MQ=255)
                                                    ctgctgTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGtc                                                                                             <  1:49892/150‑2 (MQ=255)
                                                     tgctgTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGCTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGt                                                                                              >  1:149479/1‑148 (MQ=255)
                                                          ttCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGa                                                                                        >  2:160980/1‑149 (MQ=255)
                                                               cTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTa                                                                                  <  1:160980/150‑1 (MQ=255)
                                                                ttCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTAc                                                                                 >  1:62505/1‑150 (MQ=255)
                                                                    cTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGt                                                                             <  1:142050/150‑1 (MQ=255)
                                                                     tCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGta                                                                            >  1:142210/1‑150 (MQ=255)
                                                                          cTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGtactac                                                                        >  2:147552/1‑149 (MQ=255)
                                                                            cgcgCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACaag                                                                     >  1:124849/4‑149 (MQ=255)
                                                                              cgCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACaag                                                                     >  2:162050/2‑147 (MQ=255)
                                                                               gCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACaagcag                                                                  <  1:162850/150‑5 (MQ=255)
                                                                                    cgccgAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACgg                                                             <  2:140583/146‑1 (MQ=255)
                                                                                      ccgAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACggct                                                           <  1:62865/148‑3 (MQ=255)
                                                                                       cgAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTa                                                          >  1:148319/2‑150 (MQ=255)
                                                                                        gAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGCAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGCCAGCGAACTCTACGAGTATTACAAGCAGCACGGCTAt                                                         >  1:28535/1‑150 (MQ=255)
                                                                                             gggCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCTGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAAc                                                    <  2:168965/150‑1 (MQ=255)
                                                                                               gCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCg                                                  <  2:126244/150‑1 (MQ=255)
                                                                                                 gtgtTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTg                                                >  1:167824/1‑150 (MQ=255)
                                                                                                 gtgtTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTGCTGTATCTGAAATCTACCGGAAGAAGAAAGAGAACGGATAGGAAACCGTg                                                >  1:53826/1‑150 (MQ=255)
                                                                                                 gtcttCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGtc                                                >  1:19165/4‑149 (MQ=255)
                                                                                                   gtTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGAACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGt                                              >  1:35894/1‑150 (MQ=255)
                                                                                                          atttcGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATgg                                         <  2:38725/145‑1 (MQ=255)
                                                                                                            ttcGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGgcg                                      >  2:173326/2‑149 (MQ=255)
                                                                                                            ttcGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGgcg                                      >  1:156859/2‑149 (MQ=255)
                                                                                                              cGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAACGAGCACGGTTATGAAAACGTGGTTATGGgcgc                                     >  1:13499/1‑148 (MQ=255)
                                                                                                                 cgttcgttGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCt                                <  2:145644/150‑1 (MQ=255)
                                                                                                                 cgttcgttGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCt                                >  2:54252/1‑150 (MQ=255)
                                                                                                                 cgttcgttGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGAGCAAGCt                                >  1:56432/1‑150 (MQ=255)
                                                                                                                  gTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCtt                               <  1:192567/150‑1 (MQ=255)
                                                                                                                   ttcgttgGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTc                              <  1:122670/147‑1 (MQ=255)
                                                                                                                     cgttgGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTcc                             <  2:47141/148‑1 (MQ=255)
                                                                                                                         ggCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAc                        >  1:62989/1‑150 (MQ=255)
                                                                                                                          gCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcg                       >  1:7906/1‑149 (MQ=255)
                                                                                                                           ccGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcgt                      >  2:165241/1‑148 (MQ=255)
                                                                                                                           ccGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcgt                      >  1:48293/1‑148 (MQ=255)
                                                                                                                           ccGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcgt                      >  1:157449/1‑148 (MQ=255)
                                                                                                                            cGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcgtt                     >  1:181283/1‑147 (MQ=255)
                                                                                                                               atcctcgACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGc                  >  2:139394/7‑150 (MQ=255)
                                                                                                                                tcctcgACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCg                 <  1:80813/145‑1 (MQ=255)
                                                                                                                                          gTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGa       <  1:115272/150‑1 (MQ=255)
                                                                                                                                           tACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGaa      >  1:183717/1‑150 (MQ=255)
                                                                                                                                             cAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACt    >  1:3713/1‑150 (MQ=255)
                                                                                                                                               aaGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTgg  <  1:54252/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGG  >  CP009273/8654‑8946

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGG  >  CP009273/8653‑8946
                                                                                                                                                    |                                                                                                                                                 
gcattcgcgccgctGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCA                                                                                                                                                  >  SRR12117902.48021/15‑150 (MQ=60)
     CGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCG                                                                                                                                            >  SRR12117902.5041/1‑151 (MQ=60)
                      GGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCA                                                                                                                           <  SRR12117902.137230/151‑1 (MQ=60)
                        AAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACC                                                                                                                         >  SRR12117902.53683/1‑151 (MQ=60)
                          AAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGG                                                                                                                       >  SRR12117902.68853/1‑151 (MQ=60)
                                GGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGA                                                                                                                   >  SRR12117902.109447/1‑149 (MQ=60)
                                 GCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGA                                                                                                                <  SRR12117902.195034/151‑1 (MQ=60)
                                   ATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATC                                                                                                              <  SRR12117902.39427/151‑1 (MQ=60)
                                       ACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGG                                                                                                          <  SRR12117902.161135/151‑1 (MQ=60)
                                             ACCTGAGGCTGCTGTTCTCCTTCGCTCAGGCACGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGAGTGGTACAAAGCCAACTCCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGT                                                                                                    <  SRR12117902.182408/151‑1 (MQ=60)
                                              CCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTT                                                                                                   <  SRR12117902.8591/151‑1 (MQ=60)
                                               CTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTT                                                                                                  >  SRR12117902.172055/1‑151 (MQ=60)
                                                TGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTC                                                                                                 >  SRR12117902.175353/1‑151 (MQ=60)
                                                   CCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGT                                                                                              <  SRR12117902.103654/151‑1 (MQ=60)
                                                    GCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtc                                                                                             <  SRR12117902.50811/151‑5 (MQ=60)
                                                      TGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGCTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtc                                                                                             >  SRR12117902.152197/1‑145 (MQ=60)
                                                               CCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTA                                                                                  <  SRR12117902.163914/151‑1 (MQ=60)
                                                                 TTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacg                                                                                >  SRR12117902.63642/1‑134 (MQ=60)
                                                                    GCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagt                                                                             <  SRR12117902.144632/151‑21 (MQ=60)
                                                                      TCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtat                                                                           >  SRR12117902.144796/1‑129 (MQ=60)
                                                                             CGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtattacaagc                                                                    >  SRR12117902.127127/1‑122 (MQ=60)
                                                                               cgcttgcgccGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtattacaagcag                                                                  <  SRR12117902.165827/141‑32 (MQ=60)
                                                                                      GCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtattacaagcagcacggct                                                           <  SRR12117902.64004/151‑39 (MQ=60)
                                                                                        CGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtattacaagcagcacggctat                                                         >  SRR12117902.151021/1‑111 (MQ=60)
                                                                                         GAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGCAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgccagcgaactctacgagtattacaagcagcacggctatg                                                        >  SRR12117902.29042/1‑110 (MQ=60)
                                                                                                  GTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGG                                               >  SRR12117902.170903/1‑151 (MQ=60)
                                                                                                  GTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTGCTGTATCTGAAATCTACCggaagaagaaagagaacggataggaaaccgtgg                                               >  SRR12117902.54809/1‑118 (MQ=60)
                                                                                                  GTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtattacaagcagcacggctatgaaaccgtcg                                               >  SRR12117902.19490/1‑101 (MQ=60)
                                                                                                    GTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGAACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTT                                             >  SRR12117902.36543/1‑151 (MQ=60)
                                                                                                             TTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGC                                     >  SRR12117902.159718/1‑150 (MQ=44)
                                                                                                               CGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAACGAGCACGGCTATGAAAACGTGGTTATGGGCGCA                                    >  SRR12117902.13718/1‑149 (MQ=60)
                                                                                                                  CGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTT                               <  SRR12117902.196060/151‑1 (MQ=47)
                                                                                                                  CGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGAGCAAGCTT                               >  SRR12117902.57462/1‑151 (MQ=39)
                                                                                                                   GTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTC                              <  SRR12117902.124908/151‑1 (MQ=48)
                                                                                                                          GGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACG                       >  SRR12117902.64128/1‑151 (MQ=53)
                                                                                                                           GCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGT                      >  SRR12117902.8054/1‑151 (MQ=51)
                                                                                                                            CCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACgtt                     >  SRR12117902.160320/1‑148 (MQ=50)
                                                                                                                            CCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACgtt                     >  SRR12117902.49185/1‑148 (MQ=50)
                                                                                                                             CGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACgttg                    >  SRR12117902.184564/1‑147 (MQ=51)
                                                                                                                                ATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCG                 <  SRR12117902.82279/151‑1 (MQ=50)
                                                                                                                                          GGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGA       <  SRR12117902.117369/151‑1 (MQ=54)
                                                                                                                                            TACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAAC     >  SRR12117902.187068/1‑151 (MQ=54)
                                                                                                                                              CAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTG   >  SRR12117902.3764/1‑151 (MQ=54)
                                                                                                                                               AAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGG  <  SRR12117902.55250/151‑1 (MQ=56)
                                                                                                                                                    |                                                                                                                                                 
GTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGG  >  CP009273/8653‑8946

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: