Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 1,932,640 +C intergenic (‑107/+43) rluA ← / ← minE 23S rRNA (746) and tRNA (32) pseudouridine synthase/Cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,932,6381.C96.4% 184.3 / ‑2.9 57intergenic (‑105/+45)rluA/minE23S rRNA (746) and tRNA (32) pseudouridine synthase/Cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (28/26);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

AGGATGATGCCGGAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAGAAGTCAAAAAG  >  AE015451/1932340‑1932931
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                       
gtgtgctcttccgatctttCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcc                                                                                                                                                                                                                                                                                                       <  1:1181817/285‑1 (MQ=255)
   atgatgCCGGAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACACAGGCTATGCGTCCCCTGGCGGAGCGGGCTTGcccgc                                                                                                                                                                                                                                                                                                    >  1:2267273/1‑298 (MQ=255)
        gCCGGAGGTTTCCCAGTCCCGCCGGTGCACGATCAGCTCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCGTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAgg                                                                                                                                                                                                                                                                                               <  2:472437/301‑1 (MQ=255)
          cGGAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGTATCGGGCTTGCCCGCGAAGGGc                                                                                                                                                                                                                                                                                             >  2:835665/1‑301 (MQ=255)
                 ttCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCTGGCAAGCCCGCGCCCACAAAGGATATGCGTCCCCTGTAGGGGCGGGCTTGcgcgatcggaagagcggcg                                                                                                                                                                                                                                                                                      >  2:1181369/1‑283 (MQ=255)
                        cccAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGcccc                                                                                                                                                                                                                                                                               <  1:2239722/300‑1 (MQ=255)
                                              tctcggtgaaGCCGGTTTTCTGCAGACGGGGACACAGGCATTACTTGTTTCACTCGGGACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTAGATCACCAGGTTGTCGGCGTATTCGAAGGGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTACTACAAAGGATATGCGTCCCCTGGAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCGCCAAATTGTCGACTACGCGGATCa                                                                                                                                                                                                                                                         <  1:61453/292‑1 (MQ=255)
                                              tatcgtggtggccggtagccTGCAGACGGGTGTTCAGGCAGTCTGTGTTGTCCTCGGCACGGCCAGGTACCGCCAGCAGCAGGGTCCGTTTGTTGATCACCAGGATGGGCGCGTCTTCAAAGAGGATCTGGACATTCGACAGCAGCATGGCAGGTCTCAAGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCa                                                                                                                                                                                                                                                         <  1:463068/283‑1 (MQ=255)
                                               aTCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCATGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAACATAGCAGTAAACGCTGGGGCCGCTTTGCGGCTCTTCGCGGGCACGCCCGCTCCTACAAAGGCTCCGCGCCCCCTGTAGGAGCGGGCTTCCCCGCGAAGGGCCGCCCTGAAGCCCCACATTGTCGACTACGCGGATccc                                                                                                                                                                                                                                                        >  2:154853/1‑299 (MQ=255)
                                                 cGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCg                                                                                                                                                                                                                                                      >  1:1063017/1‑301 (MQ=255)
                                                  ggggTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAgatcggaagagcacacgtctgaactccagtc                                                                                                                                                                                                                                                     >  1:557313/1‑271 (MQ=255)
                                                         ccGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCGCGCTCCTACCAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCACGGCAGCCCCAAATTGTCGACGACGCGGATCAGCCATCCgtcc                                                                                                                                                                                                                                              >  2:1738015/1‑298 (MQ=255)
                                                           gTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAACCCCGCTCCTACCAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCCGGCCAGCCCCACATTGTGGACTACGCGGATcgggcatccgcgacc                                                                                                                                                                                                                                            >  2:263845/1‑287 (MQ=255)
                                                                 cTGCAGCCGGGTGTTCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                                                                                                                                                                                      <  2:1062627/301‑1 (MQ=255)
                                                                     aGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCCCCTACAAAGGATATGCGTCCCCCGTAGGAGCGGGCTTGCCCGCGAAGGGCCGCAAGGCAGCCCCCAATTGGCGACTACGCGGATCAGCGATCCGGCAAGGCGATgctcc                                                                                                                                                                                                                                  >  1:1260698/1‑297 (MQ=255)
                                                                                          tttgtggTCATCGGCACGGCCAGTTACGGAGAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTCCGAAGAGCATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTg                                                                                                                                                                                                              <  2:1331639/300‑1 (MQ=255)
                                                                                          gttgttgTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGAATCTGGACATTGGAGAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAGCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGAAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGAAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTg                                                                                                                                                                                                              <  1:154897/300‑1 (MQ=255)
                                                                                            tgttgtCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCACCGAGCCGGCAGGGTGATGGTCAGCCTCAGAATCGAGCAGCTGcc                                                                                                                                                                                                            >  2:1739958/1‑301 (MQ=255)
                                                                                               tgtCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTg                                                                                                                                                                                                         >  2:1802936/1‑301 (MQ=255)
                                                                                                     gcGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTGTTCGAAGAGGATCTGGACATTCGACAGCGGCATGGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtttt                                                                                                                                                                                                    <  2:1514250/300‑1 (MQ=255)
                                                                                                                aGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGAGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGata                                                                                                                                                                                         <  2:475007/301‑1 (MQ=255)
                                                                                                                  gtagcggcagtggcagGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTGGAAGAGGATCGGGACATTCGACAGCGGCATTGCAGGGCTCACGAAAAGAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGAGATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                                                                                                                                                                                       <  1:263946/289‑1 (MQ=255)
                                                                                                                       ggcagcagcagGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcat                                                                                                                                                                                  <  1:835952/299‑1 (MQ=255)
                                                                                                                           gcagcagGGTGGGTTTGTTGATCACCAGGATGGAGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAACCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACatcatc                                                                                                                                                                              <  2:1260227/301‑1 (MQ=255)
                                                                                                                                 ggggcgggttgttggacACCAGGATGGCGGGGTCTTCGAAGAGGATCTGGACATTCGACGGGGGCATGGCAGGTCTCACGAAAATAGCAGTAAGCGCTGGGGCCGTTTTGCGGCCCTTCGCAGGCAAGCCCGCTCCTACAAAGGATAGGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGcc                                                                                                                                                                        <  1:662198/285‑1 (MQ=255)
                                                                                                                                     cGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATg                                                                                                                                                                    <  1:1803662/301‑1 (MQ=255)
                                                                                                                                         ttgttgATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGTGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTCGCCGATGTTca                                                                                                                                                                >  1:622226/1‑301 (MQ=255)
                                                                                                                                                     aGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTCTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCCCCAGAATCCAGCAACTGCTCTGGTTTTCCCGCTCGATCTGCACatcatctctgtcgatggtcaacttttttctgt                                                                                                                                                    >  2:138357/1‑275 (MQ=255)
                                                                                                                                                            gggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTc                                                                                                                                             <  1:1738714/300‑1 (MQ=255)
                                                                                                                                                            cggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGCGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGTCTGCAAGGCATCCCCAAATTGTCGACTACTCGGATCAGCGATCCGGCACGGGGATGGTCAGCTCCATAATCGAGCAGCCGCCCTGGTGCTCCAGCTCCATATGCAAATCATGGCTGCCGATGTTCACATGTTTGCCGATCACTTc                                                                                                                                             >  2:2100655/1‑301 (MQ=255)
                                                                                                                                                             ggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCCAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGGGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACCTTTTTGCGGATCCCTTcc                                                                                                                                            >  1:1554702/1‑301 (MQ=255)
                                                                                                                                                                ctggagttcagacgtgtgctcttccgatctcGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCcagc                                                                                                                                         <  1:2325629/271‑1 (MQ=255)
                                                                                                                                                                     cGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCCGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCCCGATATGCACATCATCCTTGCCTATGTTCACATATTTGCGGATCACTTCCAGCAGttc                                                                                                                                    >  1:1491542/1‑301 (MQ=255)
                                                                                                                                                                      gAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCTGCTCGATATGCACATCATCGTTGCCGATGTTCCCATATTTGCGGATCACTTCCAGCAGttct                                                                                                                                   >  1:743067/1‑301 (MQ=255)
                                                                                                                                                                         gagGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGttctttt                                                                                                                                >  2:1197118/1‑300 (MQ=255)
                                                                                                                                                                          agGACCTGGCCATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGttctttct                                                                                                                               <  1:1740657/301‑1 (MQ=255)
                                                                                                                                                                               cagacgtgtgcgcttccgatctGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAgc                                                                                                                          <  1:1750786/280‑1 (MQ=255)
                                                                                                                                                                               cTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGACGTTCACATATTTGCGGATCACTTCCAGCCGTTCTTTCTGCacc                                                                                                                          >  1:2078921/1‑299 (MQ=255)
                                                                                                                                                                                              cGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTTCCGATGTTCACATATTTGCGGATCACTTCcagcagatcggaagagcgtcgtgtagggaaagag                                                                                                           >  2:2324695/1‑273 (MQ=255)
                                                                                                                                                                                                     gCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCAATTCCAGCAGTTCTTTCTGCaccagatcggaagatcgatgggtag                                                                                                    >  2:1750087/1‑277 (MQ=255)
                                                                                                                                                                                                           ctcACGAAAATAGCAGTAACCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGtt                                                                                              <  1:1197566/301‑1 (MQ=255)
                                                                                                                                                                                                                    aTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGTTAGTCGGGTTCGCTGCGTTTGCCgcgcc                                                                                     >  2:632942/1‑300 (MQ=255)
                                                                                                                                                                                                                              aCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCGCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCAc                                                                           <  1:138400/301‑1 (MQ=255)
                                                                                                                                                                                                                                        cGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAACGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATCCACATCATCGTTGCCGATGTTCCCATATTTGCGGATCACTTCCACCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCTCGCTCGCTCGCCACgatgttgtgt                                                                 >  2:2254411/1‑295 (MQ=255)
                                                                                                                                                                                                                                               cGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACCTATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCAGCATGCTCTGTAGACGcc                                                          >  1:757401/1‑300 (MQ=255)
                                                                                                                                                                                                                                                   ccTTAGCGGGCGAGCCCGCTCCTACAAAGGATATGCGGCCCCTGTAGGAGCTGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAATTTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGCTGATCTGTAGACGCTCttt                                                      <  1:2101495/301‑1 (MQ=255)
                                                                                                                                                                                                                                                     acgtgtgctcttccgatctcTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCg                                                   <  1:901799/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ctttccctacacgacgctcttccgatcTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCt                                      <  2:276173/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccggtacgaaaaaggatattgttCCCCTGTAGGAGCGGGCTTGTCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCt                                      <  2:757131/279‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCCGGTAGTCCGGTTCGCTGCGTTTGCCGCGCTCGTTCGCCACGATCATCTGTAGCCGCTCttccgagatcggaagagcgtcgtg                                 >  2:901489/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                              aGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCACGCAGTTCTTTCTGCAGCGGGGGCAGGTAGTCCGTTTCCCTGCGTTGGCCGCGCTCGTGCGCCAGGATGGTCTGTAGCCGCTCTTTCGCTACCGAccccgtggttggtttctt                          >  2:1896106/1‑284 (MQ=255)
                                                                                                                                                                                                                                                                                  ttggCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTtctgtctg                      <  2:742797/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   atGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCCTTTCCGCTACCGCCGCGCTGGTCTGTTTCTGTCTGc                     >  2:1674362/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                      ggtcggcgggaGGAGCGGACGTGCCAGCGCAGGGCTGCAAGGCAGCCCCAAACTGTCGACTACGCGGATGAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGGTGCCGATGTTCACATATTTGCGGATCGCTTCCAGCAGTGCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGCTCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCAc                  <  1:398217/291‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       gTCCCCTGTAGGGGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACg                 <  2:1490953/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          cccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGaaa              <  1:633177/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            cTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCGCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAaaga            >  2:305546/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                              gTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCtagatcggaagagcacacgtctgaactcc          >  1:276274/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                                                aGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTTCGCCACGATGATCTTTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAGAAGTc        >  1:2219120/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                      gggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAGAAGTCAAAAAg  >  1:127713/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                       
AGGATGATGCCGGAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAGAAGTCAAAAAG  >  AE015451/1932340‑1932931

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