Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 211,090 T→C 6.8% A5532A (GCT→GCC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451211,0900TC 6.8% 264.7 / 4.6 90A5532A (GCT→GCCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (45/34);  new base C (2/4);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGTTGGC  >  AE015451/210860‑211317
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                   
aaCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                                                                                                                 <  2:661265/235‑1 (MQ=34)
 aCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGtg                                                                                                                                                                                                                               <  1:1319682/236‑2 (MQ=25)
    tGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGTGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACt                                                                                                                                                                                                                            <  1:1017546/236‑1 (MQ=25)
     gCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTg                                                                                                                                                                                                                           >  1:688635/1‑236 (MQ=34)
              cggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCgg                                                                                                                                                                                                                     <  2:903241/233‑1 (MQ=35)
                          ggTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGACAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTa                                                                                                                                                                                                           <  1:567669/231‑1 (MQ=18)
                                               cgacgaAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTg                                                                                                                                                                                  <  2:722426/235‑1 (MQ=38)
                                                   gaAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGc                                                                                                                                                                                                    <  2:1319533/213‑1 (MQ=35)
                                                    aaGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCg                                                                                                                                                                                                                      <  1:1008/194‑1 (MQ=34)
                                                         gtggCCAAGCTGACTGCGACCCAGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGa                                                                                                                                                                           <  1:1381531/232‑1 (MQ=39)
                                                               aaGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAAc                                                                                                                                                                    >  2:1348842/1‑233 (MQ=255)
                                                                    gACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAg                                                                                                                                                            >  1:1288360/1‑236 (MQ=40)
                                                                    gACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAg                                                                                                                                                            >  1:634453/1‑236 (MQ=255)
                                                                      cTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACTCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATc                                                                                                                                                                         <  2:1246334/221‑1 (MQ=38)
                                                                      cTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTg                                                                                                                                                                                                                           <  2:289344/171‑1 (MQ=34)
                                                                         cGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGGTGGAg                                                                                                                                                            >  1:343486/1‑231 (MQ=37)
                                                                          gACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                                                                                                                 <  2:792103/161‑1 (MQ=33)
                                                                           aCCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGccc                                                                                                                                                                                                                        >  2:723016/1‑169 (MQ=33)
                                                                              ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTc                                                                                                                                                                                                             >  1:34037/1‑177 (MQ=33)
                                                                              ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGt                                                                                                                                                  >  1:1153009/1‑236 (MQ=38)
                                                                               cGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATc                                                                                                                                                                         >  2:1110313/1‑212 (MQ=39)
                                                                               cGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATc                                                                                                                                                                         <  1:1110316/212‑1 (MQ=38)
                                                                                  cGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGa                                                                                                                                                                           >  1:1518218/1‑207 (MQ=40)
                                                                                  cGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGa                                                                                                                                                                           <  2:1518215/207‑1 (MQ=40)
                                                                                  cGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACCTCCCCGGACAACgcc                                                                                                                                                                                                             >  2:370431/1‑171 (MQ=18)
                                                                                  cGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCAGGACAACGt                                                                                                                                                                                                              >  2:1063356/1‑172 (MQ=18)
                                                                                               cggcgAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGc                                                                                                                                    <  2:684431/233‑1 (MQ=37)
                                                                                                  cgAGATCACCTACACCATCACGCTGACCACCAAAGACGGTCTGCCAATCAACAACCACAGTGCATCGACCTTCACGCTGAGCGACGGCAAAACCGTCACCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCAc                                                                                                                                                                                                                             <  2:96056/141‑1 (MQ=255)
                                                                                                    cgaTCACCTACACCACCACGCTGACCAACAAAGCCGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTa                                                                                                                                             <  2:1471309/218‑1 (MQ=11)
                                                                                                    agaTCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAACACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGGCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTg                                                                                                                                                                                  >  1:1434523/1‑182 (MQ=255)
                                                                                                    agaTCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCCACGTCGGCACCCACGACCCCGTCgcg                                                                                                                                                                                  >  2:304223/1‑180 (MQ=255)
                                                                                                                 cATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGAcg                                                                                                                <  1:1348845/235‑2 (MQ=39)
                                                                                                                 cATCACCCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGAcc                                                                                                                <  1:950233/235‑1 (MQ=37)
                                                                                                                  aTCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGAcccc                                                                                                              <  1:1087311/236‑1 (MQ=37)
                                                                                                                  aTCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGAcc                                                                                                                <  2:891373/234‑1 (MQ=38)
                                                                                                                   tCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGAcccc                                                                                                              >  2:1055443/1‑235 (MQ=38)
                                                                                                                      cGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTg                                                                                                                         >  2:967794/1‑221 (MQ=38)
                                                                                                                      cGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                                                                                                                                      >  2:689792/1‑128 (MQ=17)
                                                                                                                       gCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTc                                                                                                                                                                             <  1:1332966/168‑1 (MQ=37)
                                                                                                                       gCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTc                                                                                                                                                                             >  2:1332963/1‑168 (MQ=39)
                                                                                                                         tGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGc                                                                                                       <  2:634451/236‑1 (MQ=38)
                                                                                                                           aCCAACAAAGACGGTCTGCCTATCAACAACCATGCTGCATTGACCTTTACGCTGATCGACGGCAAAACCGACATCACCGTCCCGGCTAACGGCACCACCGGCACTGCCACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTc                                                                                                                                                                             <  1:1008075/164‑1 (MQ=11)
                                                                                                                           aCCAACAAAGACGGTCTGCAAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGGCATCACCGTGGAGGCCAACGGCACCGTGGGCACTGCCACTGTCACTGACCCg                                                                                                                                                                                                                      >  2:1447712/1‑123 (MQ=11)
                                                                                                                            ccAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCCGACAACGTCTACgg                                                                                                                                                                                                        >  2:1029366/1‑135 (MQ=17)
                                                                                                                              aacaaAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTAGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTccctt                                                                                                                                                                                 >  2:1000472/1‑153 (MQ=2)
                                                                                                                              aacaaAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCGACCGGCACCGTGGGCACTGCGAcc                                                                                                                                                                                                                                  >  2:171706/1‑108 (MQ=255)
                                                                                                                              aacaaAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTAGATCGCAACCgat                                                                                                                                                                >  2:308979/1‑172 (MQ=25)
                                                                                                                                       ggccTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGtt                                                                                                                                                                <  1:1374781/162‑1 (MQ=17)
                                                                                                                                        gTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCAcc                                                                                           >  2:856168/1‑233 (MQ=38)
                                                                                                                                        gTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCAcc                                                                                           >  2:978098/1‑233 (MQ=37)
                                                                                                                                        gTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACcgac                                                                                        <  1:856170/236‑1 (MQ=38)
                                                                                                                                        gTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACcga                                                                                         >  1:1099634/1‑235 (MQ=38)
                                                                                                                                             ccgaTCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGcc                                                                                   <  1:40397/233‑1 (MQ=11)
                                                                                                                                                       acCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGGCGGCACCAACGACCCTGTCGTGATGTCGATc                                                                                                                                                                         <  2:1356773/140‑1 (MQ=17)
                                                                                                                                                          acaGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGt                                                                                                                                                          >  2:695126/1‑152 (MQ=32)
                                                                                                                                                          acaGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGt                                                                                                                                                          <  1:695128/152‑1 (MQ=31)
                                                                                                                                                                cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGt                                                                                                                                                                 >  2:590170/1‑139 (MQ=32)
                                                                                                                                                                cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCg                                                                                                                                                >  1:51426/1‑156 (MQ=31)
                                                                                                                                                                cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAAccc                                                                 >  1:1016994/1‑235 (MQ=38)
                                                                                                                                                                cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCg                                                                >  1:576901/1‑236 (MQ=37)
                                                                                                                                                                   tGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCgg                                                               >  2:1044325/1‑234 (MQ=37)
                                                                                                                                                                   tGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTCGAGGTGCGGACGTCGGTAAGTTCGAg                                                                                                                                     >  2:532124/1‑164 (MQ=14)
                                                                                                                                                                      ccTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTc                                                                                                          >  1:101320/1‑188 (MQ=37)
                                                                                                                                                                      ccTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGggcgg                                                           >  1:1299325/1‑235 (MQ=38)
                                                                                                                                                                        ttCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACa                                                                                                                     <  2:1402066/175‑1 (MQ=37)
                                                                                                                                                                        ttCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACa                                                                                                                     >  1:1402069/1‑175 (MQ=21)
                                                                                                                                                                           aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTCCCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCg                                                                                                                                                                  >  2:1018995/1‑127 (MQ=17)
                                                                                                                                                                             gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACg                                                   >  1:628939/1‑236 (MQ=38)
                                                                                                                                                                             ccTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCCCTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCAcccc                                                                          <  1:330133/212‑1 (MQ=25)
                                                                                                                                                                                 aGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAgg                                               <  1:1404294/236‑1 (MQ=39)
                                                                                                                                                                                  gCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACg                                                   >  2:1372727/1‑231 (MQ=38)
                                                                                                                                                                                    gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGCCGCTg                                                                                                                         >  2:521447/1‑159 (MQ=14)
                                                                                                                                                                                    gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGt                                                                                               >  2:471335/1‑185 (MQ=37)
                                                                                                                                                                                    gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGc                                                                                    >  1:298787/1‑196 (MQ=37)
                                                                                                                                                                                     aCGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGc                                              <  2:1294415/233‑1 (MQ=21)
                                                                                                                                                                                               cGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAg                                  <  1:323396/235‑1 (MQ=37)
                                                                                                                                                                                               cGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAg                                  <  1:738020/235‑1 (MQ=37)
                                                                                                                                                                                                gTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAg                                  <  2:1080657/234‑1 (MQ=37)
                                                                                                                                                                                                    tcaCCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGACGGCGACCTGGTCAAGGt                                >  1:223512/1‑232 (MQ=37)
                                                                                                                                                                                                      aCCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCAt                          >  1:671334/1‑236 (MQ=39)
                                                                                                                                                                                                           gCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCAcc                                                                                           <  1:329109/166‑1 (MQ=38)
                                                                                                                                                                                                            ccGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACcgg                                                                                         >  2:1386744/1‑166 (MQ=39)
                                                                                                                                                                                                              ggCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCg                  <  1:161883/236‑1 (MQ=37)
                                                                                                                                                                                                               gCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCAcc                                                                            >  2:452610/1‑177 (MQ=38)
                                                                                                                                                                                                               gCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAg                                  <  1:615748/219‑1 (MQ=38)
                                                                                                                                                                                                                  aaCGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAAcc                                                                  >  2:1315219/1‑184 (MQ=39)
                                                                                                                                                                                                                       cACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCTGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCa              >  2:874701/1‑231 (MQ=37)
                                                                                                                                                                                                                             gggCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGt                                >  2:1204840/1‑207 (MQ=38)
                                                                                                                                                                                                                              ggCACTGCCACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTgg                                       <  1:683876/199‑1 (MQ=16)
                                                                                                                                                                                                                               gCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGc  >  1:719993/1‑235 (MQ=14)
                                                                                                                                                                                                                                      tACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCCAGGTCACCATCACGGCCGACCAGACTTCGtt     >  1:235089/1‑225 (MQ=25)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                   
AACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGTTGGC  >  AE015451/210860‑211317

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