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 F1 I1 R1 343 268.7 10873516 97.0% 10547310 150.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,957,233 2 bp→GT coding (666‑667/756 nt) RPPX_21945 → ATPase AAA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,957,2330TG99.5% 495.0 / NA 205R222R (CGT→CGGRPPX_21945ATPase AAA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (101/103);  total (102/103)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*CP0099744,957,2340CT99.5% 503.8 / NA 205L223L (CTG→TTG) RPPX_21945ATPase AAA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (102/102);  minor base G (0/1);  total (102/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

TACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTCACCAGCCAGATGCCGGTGGACAAATGGCATGCGCTGATCGGCGATCCGACCTTGGGCGACGCGATCCTCGACCGTCTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGA  >  CP009974/4957089‑4957375
                                                                                                                                                ||                                                                                                                                             
tACCGCAGCGCAACGCCGCGAGATGCTCGAACTGCTGGACGACCGCTACGGCAGCCGCTCGACACGGGTGACCAGCTAGATGCCGGTGGATAAATGGCATGCGCTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGGCCGGTTGGTg                                                                                                                                          <  2:4392034/151‑1 (MQ=11)
 aCCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGTTTGGTGc                                                                                                                                         >  1:246630/1‑151 (MQ=11)
 aCCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGc                                                                                                                                         <  1:1522477/151‑1 (MQ=11)
  ccGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGca                                                                                                                                        >  1:1748512/1‑151 (MQ=11)
   cGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGcac                                                                                                                                       >  2:5376444/1‑151 (MQ=11)
   cGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGcac                                                                                                                                       <  1:2594287/151‑1 (MQ=11)
   cGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGcac                                                                                                                                       >  1:969431/1‑151 (MQ=11)
    gCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGcaca                                                                                                                                      >  1:3171524/1‑151 (MQ=11)
    gCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGcaca                                                                                                                                      >  1:3005618/1‑151 (MQ=11)
    gCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGcaca                                                                                                                                      >  1:1781870/1‑151 (MQ=11)
     cggcggAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACaa                                                                                                                                     <  2:1925814/149‑1 (MQ=11)
      aGCGCGCCGCCGGGAGATGCGCGAACTGCTGGACGCCCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAAc                                                                                                                                    <  2:2723143/151‑1 (MQ=11)
       gcgcAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACg                                                                                                                                   <  1:2945499/151‑1 (MQ=11)
            aCGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTac                                                                                                                              >  2:679263/1‑150 (MQ=11)
            aCGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTAt                                                                                                                              <  2:3562897/151‑1 (MQ=11)
            aCGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTAt                                                                                                                              <  1:1575432/151‑1 (MQ=11)
             cgccgcGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATc                                                                                                                             <  1:1778899/151‑1 (MQ=11)
              gccgcGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCg                                                                                                                            <  1:2050262/151‑1 (MQ=11)
                                            gCTACGGGAACCGCGGGACAGTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCaa                                                                                              <  2:3618643/151‑1 (MQ=11)
                                            gCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCaa                                                                                              <  1:5078944/151‑1 (MQ=11)
                                              tACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGc                                                                                            >  2:4637018/1‑151 (MQ=11)
                                               aCGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCg                                                                                           >  1:1223379/1‑151 (MQ=11)
                                                cGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGt                                                                                          >  1:3375107/1‑151 (MQ=11)
                                                cGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGt                                                                                          <  1:2832001/151‑1 (MQ=11)
                                                 ggCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTg                                                                                         <  1:1348072/151‑1 (MQ=11)
                                                  gCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGc                                                                                        >  2:4549912/1‑151 (MQ=11)
                                                  gCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGc                                                                                        >  2:4799528/1‑151 (MQ=11)
                                                  gCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGc                                                                                        <  1:4717109/151‑1 (MQ=11)
                                                  gCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGc                                                                                        >  2:3486355/1‑151 (MQ=11)
                                                    caCCGCTGGACACTGGTGACCAGCGAGATGCCGGTGGATAAATGGAATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATACTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcga                                                                                      <  2:2309267/150‑1 (MQ=11)
                                                    aaCCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcga                                                                                      >  1:4556800/1‑151 (MQ=11)
                                                    aaCCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcga                                                                                      >  1:2440923/1‑151 (MQ=11)
                                                     aCCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcgac                                                                                     <  2:1757485/151‑1 (MQ=11)
                                                     aCCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcgac                                                                                     <  1:2915411/151‑1 (MQ=11)
                                                     aCCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcgac                                                                                     >  1:3567928/1‑151 (MQ=11)
                                                        gCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGaa                                                                                  <  2:1200151/151‑1 (MQ=11)
                                                        gCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGaa                                                                                  <  2:4966805/151‑1 (MQ=14)
                                                         cTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGaaa                                                                                 >  2:2957711/1‑151 (MQ=11)
                                                         cTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGaaa                                                                                 >  1:3755854/1‑151 (MQ=14)
                                                         cTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGaaa                                                                                 >  1:4108282/1‑151 (MQ=14)
                                                            gACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTg                                                                              >  1:1672887/1‑151 (MQ=11)
                                                             acacTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGa                                                                             >  1:1174772/1‑151 (MQ=11)
                                                               acTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGAcg                                                                           >  1:1011018/1‑151 (MQ=11)
                                                                  ggtgaCCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGAcggcg                                                                        <  2:789084/147‑1 (MQ=14)
                                                                    tgaCCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACggcggc                                                                      >  1:1764369/3‑151 (MQ=14)
                                                                       ccagccagATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAgg                                                                   <  1:3437614/151‑1 (MQ=14)
                                                                       ccagccagATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAgg                                                                   <  1:5259944/151‑1 (MQ=14)
                                                                        cagccagATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAggg                                                                  >  2:516987/1‑151 (MQ=14)
                                                                          gccagATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGAc                                                                <  1:1669684/151‑1 (MQ=14)
                                                                           ccagATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGAcg                                                               >  2:2097484/1‑150 (MQ=14)
                                                                             agATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTc                                                             <  1:1912427/151‑1 (MQ=14)
                                                                              gATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCa                                                            >  1:2647206/1‑151 (MQ=14)
                                                                                    ggTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTaa                                                      <  2:168090/151‑1 (MQ=14)
                                                                                     gTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAAc                                                     <  1:4643308/151‑1 (MQ=14)
                                                                                     gTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAAc                                                     >  1:3539128/1‑151 (MQ=14)
                                                                                      tGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACa                                                    <  1:2402768/151‑1 (MQ=14)
                                                                                      tGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACa                                                    <  1:2187479/151‑1 (MQ=14)
                                                                                      tGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACa                                                    <  2:1119484/151‑1 (MQ=14)
                                                                                      tGGATAAATGGCATGCACTGACCGGCGATCCGACCTTGTGCGATGTGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCCAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTacaa                                                    >  2:5294264/1‑148 (MQ=14)
                                                                                       ggATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACaa                                                   >  1:4939937/1‑151 (MQ=14)
                                                                                         ataaaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATg                                                 >  1:3033622/3‑151 (MQ=14)
                                                                                          taaaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGc                                                <  2:3758411/150‑1 (MQ=255)
                                                                                           aagtGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCa                                               <  2:4205772/148‑1 (MQ=14)
                                                                                           aaaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCa                                               <  2:1535545/151‑1 (MQ=14)
                                                                                           aaaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCa                                               >  1:4568834/1‑151 (MQ=14)
                                                                                            aaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCaa                                              <  1:5118845/151‑1 (MQ=14)
                                                                                            aaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCaa                                              <  2:1748512/151‑1 (MQ=14)
                                                                                              tGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAAcc                                            <  2:676415/151‑1 (MQ=255)
                                                                                               ggCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAAccc                                           >  1:1571131/1‑151 (MQ=14)
                                                                                                     cacTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTc                                     <  1:3260497/149‑1 (MQ=255)
                                                                                                      acTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCg                                    <  1:127683/150‑1 (MQ=255)
                                                                                                      acTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCg                                    >  1:3954123/2‑151 (MQ=255)
                                                                                                      acTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCg                                    <  1:1309492/150‑1 (MQ=255)
                                                                                                        tGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCt                                  <  2:2469334/151‑1 (MQ=255)
                                                                                                        tGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCt                                  <  2:1599354/151‑1 (MQ=255)
                                                                                                         gATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCTCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTg                                 <  1:252702/151‑1 (MQ=255)
                                                                                                         gATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTg                                 >  1:3583038/1‑151 (MQ=255)
                                                                                                          aTCGGCGATCCGACCGTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTgc                                <  1:5011631/151‑1 (MQ=255)
                                                                                                           tCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTgcg                               >  1:3276447/1‑151 (MQ=255)
                                                                                                            cGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTgcgc                              >  1:3781537/1‑151 (MQ=255)
                                                                                                             ggCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTgcgcc                             >  2:253039/1‑150 (MQ=255)
                                                                                                             ggCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCt                             <  2:5144022/151‑1 (MQ=255)
                                                                                                               cGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTcc                           >  1:5376139/1‑151 (MQ=255)
                                                                                                                gATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCg                          <  1:4483819/151‑1 (MQ=255)
                                                                                                                   ccGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATTGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACt                       >  2:1402888/1‑151 (MQ=255)
                                                                                                                   ccGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACt                       <  2:1331585/151‑1 (MQ=255)
                                                                                                                    cGACGTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTg                      <  2:3322231/151‑1 (MQ=255)
                                                                                                                       ccTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCt                   <  1:400419/151‑1 (MQ=255)
                                                                                                                        cTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTg                  <  2:3268529/151‑1 (MQ=255)
                                                                                                                        cTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTg                  <  1:735003/151‑1 (MQ=255)
                                                                                                                        cTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTg                  >  2:1883830/1‑151 (MQ=255)
                                                                                                                            ggCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCa              >  1:3581475/1‑151 (MQ=255)
                                                                                                                              cGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCCGCGCTCCGACTCCCTGTCCaac            >  2:1661003/1‑150 (MQ=14)
                                                                                                                                atgcgATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTg          >  1:1974645/3‑151 (MQ=255)
                                                                                                                                  gcgATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGa        >  1:3774470/1‑151 (MQ=255)
                                                                                                                                    gATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACg      >  2:55867/1‑151 (MQ=255)
                                                                                                                                     aTCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGt     <  1:679263/151‑1 (MQ=255)
                                                                                                                                     aTCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGt     <  1:2040769/151‑1 (MQ=255)
                                                                                                                                        cTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGa  >  2:3868145/1‑151 (MQ=255)
                                                                                                                                        cTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGa  >  1:603180/1‑151 (MQ=255)
                                                                                                                                                ||                                                                                                                                             
TACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTCACCAGCCAGATGCCGGTGGACAAATGGCATGCGCTGATCGGCGATCCGACCTTGGGCGACGCGATCCTCGACCGTCTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGA  >  CP009974/4957089‑4957375

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTCACCAGCCAGATGCCGGTGGACAAATGGCATGCGCTGATCGGCGATCCGACCTTGGGCGACGCGATCCTCGACCGTCTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGATCAG  >  CP009974/4957083‑4957379
                                                                                                                                                      |                                                                                                                                                  
CCCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGG                                                                                                                                                    <  SRR10918230.5414450/151‑1 (MQ=19)
ACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGG                                                                                                                                                    >  SRR10918230.910448/1‑151 (MQ=19)
       ACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGTTTGGTGC                                                                                                                                             >  SRR10918230.247466/1‑151 (MQ=19)
       ACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGC                                                                                                                                             <  SRR10918230.1529204/151‑1 (MQ=19)
        CCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCA                                                                                                                                            >  SRR10918230.1756080/1‑151 (MQ=18)
        CCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCA                                                                                                                                            >  SRR10918230.5103127/1‑151 (MQ=18)
         cgtgtgctcttccgatctCATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCAC                                                                                                                                           <  SRR10918230.4989614/133‑1 (MQ=0)
         CGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCAC                                                                                                                                           >  SRR10918230.973348/1‑151 (MQ=16)
         CGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCAC                                                                                                                                           <  SRR10918230.2605571/151‑1 (MQ=16)
         CGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCAC                                                                                                                                           <  SRR10918230.2605736/151‑1 (MQ=16)
         CGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCAC                                                                                                                                           <  SRR10918230.3668044/151‑1 (MQ=16)
          GCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACA                                                                                                                                          >  SRR10918230.3185323/1‑151 (MQ=15)
          GCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACA                                                                                                                                          >  SRR10918230.3018684/1‑151 (MQ=15)
          GCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACA                                                                                                                                          >  SRR10918230.3578312/1‑151 (MQ=15)
          GCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACA                                                                                                                                          >  SRR10918230.3017743/1‑151 (MQ=15)
          GCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACA                                                                                                                                          >  SRR10918230.569152/1‑151 (MQ=15)
          GCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACA                                                                                                                                          >  SRR10918230.1789556/1‑151 (MQ=15)
            AGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAAC                                                                                                                                        >  SRR10918230.2102677/1‑151 (MQ=8)
             GCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACG                                                                                                                                       <  SRR10918230.2958305/151‑1 (MQ=6)
                 AACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTA                                                                                                                                   >  SRR10918230.4799539/1‑151 (MQ=6)
                  ACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTAT                                                                                                                                  >  SRR10918230.1140368/1‑151 (MQ=4)
                  ACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTAT                                                                                                                                  <  SRR10918230.1582314/151‑1 (MQ=4)
                  ACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTAT                                                                                                                                  <  SRR10918230.1582348/151‑1 (MQ=4)
                   CGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATC                                                                                                                                 <  SRR10918230.1786582/151‑1 (MQ=2)
                   CGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATC                                                                                                                                 >  SRR10918230.4997529/1‑151 (MQ=2)
                   CGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATC                                                                                                                                 <  SRR10918230.3772421/151‑1 (MQ=2)
                    GCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCG                                                                                                                                <  SRR10918230.2058913/151‑1 (MQ=1)
                    GCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCG                                                                                                                                >  SRR10918230.2569080/1‑151 (MQ=1)
                     CCGCGACATGGTCGATCGGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGG                                                                                                                               <  SRR10918230.2495620/151‑1 (MQ=0)
                     CCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGG                                                                                                                               <  SRR10918230.2858449/151‑1 (MQ=0)
                      CGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGA                                                                                                                              <  SRR10918230.5230322/151‑1 (MQ=0)
                       GCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGAT                                                                                                                             >  SRR10918230.3911711/1‑151 (MQ=0)
                        CGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATC                                                                                                                            <  SRR10918230.3869964/151‑1 (MQ=0)
                                 CGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAG                                                                                                                   >  SRR10918230.4427337/1‑151 (MQ=0)
                                    ACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGC                                                                                                                <  SRR10918230.4936046/151‑1 (MQ=0)
                                     CTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCG                                                                                                               >  SRR10918230.3161385/1‑151 (MQ=0)
                                       GCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAA                                                                                                             >  SRR10918230.4884466/1‑151 (MQ=0)
                                               ACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCG                                                                                                     <  SRR10918230.4504976/151‑1 (MQ=0)
                                                 CGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCA                                                                                                   <  SRR10918230.2646463/151‑1 (MQ=0)
                                                  GCTACGGCAACCGCTCGACACTGGTGATCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAA                                                                                                  <  SRR10918230.796174/151‑1 (MQ=0)
                                                  GCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAA                                                                                                  <  SRR10918230.5102232/151‑1 (MQ=0)
                                                   CTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAG                                                                                                 >  SRR10918230.3830934/1‑151 (MQ=1)
                                                   CTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAG                                                                                                 >  SRR10918230.3032753/1‑151 (MQ=1)
                                                     ACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCG                                                                                               >  SRR10918230.1228722/1‑151 (MQ=4)
                                                      CGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGT                                                                                              <  SRR10918230.2400084/151‑1 (MQ=6)
                                                      CGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGT                                                                                              >  SRR10918230.3389720/1‑151 (MQ=6)
                                                      CGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGT                                                                                              <  SRR10918230.2844371/151‑1 (MQ=6)
                                                      CGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGT                                                                                              <  SRR10918230.1613788/151‑1 (MQ=6)
                                                       GGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTG                                                                                             <  SRR10918230.1354198/151‑1 (MQ=8)
                                                       GGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTG                                                                                             <  SRR10918230.1354887/151‑1 (MQ=8)
                                                        GCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGC                                                                                            <  SRR10918230.4738518/151‑1 (MQ=10)
                                                         CAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCG                                                                                           <  SRR10918230.3589134/151‑1 (MQ=15)
                                                          AACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGA                                                                                          >  SRR10918230.2451398/1‑151 (MQ=16)
                                                          AACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGA                                                                                          >  SRR10918230.4577187/1‑151 (MQ=16)
                                                          AACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGA                                                                                          >  SRR10918230.4577399/1‑151 (MQ=16)
                                                          AACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGA                                                                                          >  SRR10918230.5167665/1‑151 (MQ=16)
                                                           cgtgtgctcttccgatctCCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGAC                                                                                         <  SRR10918230.1761517/133‑1 (MQ=19)
                                                           ACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGAC                                                                                         >  SRR10918230.3583378/1‑151 (MQ=18)
                                                           ACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGAC                                                                                         <  SRR10918230.2928092/151‑1 (MQ=18)
                                                           ACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGAC                                                                                         <  SRR10918230.3435816/151‑1 (MQ=18)
                                                           ACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGAC                                                                                         <  SRR10918230.3057026/151‑1 (MQ=18)
                                                               CTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAA                                                                                     >  SRR10918230.3772172/1‑151 (MQ=26)
                                                               CTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAA                                                                                     >  SRR10918230.4126594/1‑151 (MQ=26)
                                                                TCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAAT                                                                                    <  SRR10918230.1774019/151‑1 (MQ=27)
                                                                  GACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTG                                                                                  >  SRR10918230.1680218/1‑151 (MQ=27)
                                                                   ACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGA                                                                                 >  SRR10918230.1179827/1‑151 (MQ=27)
                                                                    CACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGAC                                                                                >  SRR10918230.1337609/1‑151 (MQ=27)
                                                                     ACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACG                                                                               >  SRR10918230.1015163/1‑151 (MQ=27)
                                                                        GGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGATCGGTTGGTGCACAACGCTTATCGGTTCGAACTGAAGGGCGAATCGATGCGCAAGCGTCCGACGAAATTGACGctg                                                                            >  SRR10918230.1195046/1‑148 (MQ=25)
                                                                        GGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCG                                                                            <  SRR10918230.242631/151‑1 (MQ=28)
                                                                          TGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGC                                                                          >  SRR10918230.1772004/1‑151 (MQ=30)
                                                                             CCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGG                                                                       <  SRR10918230.3452454/151‑1 (MQ=36)
                                                                             CCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGG                                                                       <  SRR10918230.5284335/151‑1 (MQ=36)
                                                                                GCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGAC                                                                    >  SRR10918230.2949119/1‑151 (MQ=41)
                                                                                GCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGAC                                                                    <  SRR10918230.1676996/151‑1 (MQ=41)
                                                                                GCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGAC                                                                    <  SRR10918230.2000760/151‑1 (MQ=41)
                                                                                   AGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTC                                                                 <  SRR10918230.1920540/151‑1 (MQ=47)
                                                                                    GATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCA                                                                >  SRR10918230.4585998/1‑151 (MQ=47)
                                                                                    GATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCA                                                                >  SRR10918230.2658744/1‑151 (MQ=47)
                                                                                     ATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAG                                                               >  SRR10918230.4060474/1‑151 (MQ=47)
                                                                                      TGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGA                                                              >  SRR10918230.5076422/1‑151 (MQ=47)
                                                                                       GCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTTACGGCGGCAGGGACTTCAGAC                                                             >  SRR10918230.33795/1‑151 (MQ=46)
                                                                                        CCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACT                                                            <  SRR10918230.4691154/151‑1 (MQ=47)
                                                                                         CGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTA                                                           >  SRR10918230.231409/1‑151 (MQ=47)
                                                                                          GGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAA                                                          <  SRR10918230.719542/151‑1 (MQ=47)
                                                                                           GTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAAC                                                         <  SRR10918230.4664327/151‑1 (MQ=47)
                                                                                           GTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAAC                                                         >  SRR10918230.3554406/1‑151 (MQ=47)
                                                                                            TGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACA                                                        >  SRR10918230.2460119/1‑151 (MQ=47)
                                                                                            TGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACA                                                        <  SRR10918230.2413083/151‑1 (MQ=47)
                                                                                            TGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACA                                                        <  SRR10918230.2196744/151‑1 (MQ=47)
                                                                                             GGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAA                                                       >  SRR10918230.4962394/1‑151 (MQ=47)
                                                                                             GGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAA                                                       <  SRR10918230.3481887/151‑1 (MQ=47)
                                                                                              GATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAAT                                                      >  SRR10918230.1735267/1‑151 (MQ=48)
                                                                                              GATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAAT                                                      >  SRR10918230.4478093/1‑151 (MQ=48)
                                                                                               ATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATG                                                     >  SRR10918230.3046825/1‑151 (MQ=48)
                                                                                                 AAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCA                                                   <  SRR10918230.1846928/151‑1 (MQ=48)
                                                                                                 AAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCA                                                   >  SRR10918230.4589479/1‑151 (MQ=48)
                                                                                                  AATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAA                                                  <  SRR10918230.5142349/151‑1 (MQ=48)
                                                                                                     GGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCC                                               >  SRR10918230.1578027/1‑151 (MQ=54)
                                                                                                           CACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTC                                         <  SRR10918230.3274667/151‑1 (MQ=59)
                                                                                                            ACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCG                                        <  SRR10918230.1315403/151‑1 (MQ=59)
                                                                                                            ACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCG                                        <  SRR10918230.128122/151‑1 (MQ=59)
                                                                                                            ACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCG                                        >  SRR10918230.3971471/1‑151 (MQ=59)
                                                                                                               GATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCTCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTG                                     <  SRR10918230.253556/151‑1 (MQ=58)
                                                                                                               GATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTG                                     >  SRR10918230.3598244/1‑151 (MQ=60)
                                                                                                               GATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTG                                     >  SRR10918230.3598553/1‑151 (MQ=60)
                                                                                                                ATCGGCGATCCGACCGTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGC                                    <  SRR10918230.5034566/151‑1 (MQ=58)
                                                                                                                 TCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCG                                   >  SRR10918230.3290694/1‑151 (MQ=60)
                                                                                                                  CGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGC                                  >  SRR10918230.3797983/1‑151 (MQ=60)
                                                                                                                   GGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCT                                 >  SRR10918230.1538219/1‑151 (MQ=60)
                                                                                                                   GGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCT                                 >  SRR10918230.3797965/1‑151 (MQ=60)
                                                                                                                    GCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAAgatcggaagagcacacgtctgaa                                >  SRR10918230.897466/1‑128 (MQ=52)
                                                                                                                     CGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCC                               >  SRR10918230.5401380/1‑151 (MQ=60)
                                                                                                                     CGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCC                               <  SRR10918230.4600478/151‑1 (MQ=60)
                                                                                                                      GATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCG                              <  SRR10918230.4504112/151‑1 (MQ=60)
                                                                                                                       ATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGA                             <  SRR10918230.5287284/151‑1 (MQ=60)
                                                                                                                        TCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGAC                            <  SRR10918230.608166/151‑1 (MQ=60)
                                                                                                                          CGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTG                          >  SRR10918230.1922344/1‑151 (MQ=60)
                                                                                                                            ACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCC                        <  SRR10918230.3236384/151‑1 (MQ=60)
                                                                                                                             CCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCT                       <  SRR10918230.401778/151‑1 (MQ=60)
                                                                                                                              CTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTG                      <  SRR10918230.1308570/151‑1 (MQ=60)
                                                                                                                              CTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTG                      <  SRR10918230.737733/151‑1 (MQ=60)
                                                                                                                               TTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGT                     <  SRR10918230.3276188/151‑1 (MQ=60)
                                                                                                                                  GGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCA                  >  SRR10918230.3596986/1‑151 (MQ=60)
                                                                                                                                   GCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAA                 <  SRR10918230.1789208/151‑1 (MQ=60)
                                                                                                                                    CGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAG                >  SRR10918230.2798044/1‑151 (MQ=60)
                                                                                                                                      ATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTG              >  SRR10918230.1680203/1‑151 (MQ=60)
                                                                                                                                      ATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTG              >  SRR10918230.4017212/1‑151 (MQ=60)
                                                                                                                                      ATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTG              >  SRR10918230.1982986/1‑151 (MQ=60)
                                                                                                                                        GCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGA            >  SRR10918230.3790883/1‑151 (MQ=60)
                                                                                                                                          GATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACG          >  SRR10918230.3711692/1‑151 (MQ=60)
                                                                                                                                           ATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGT         <  SRR10918230.1254439/151‑1 (MQ=60)
                                                                                                                                           ATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGT         <  SRR10918230.681726/151‑1 (MQ=60)
                                                                                                                                           ATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGT         <  SRR10918230.2049517/151‑1 (MQ=60)
                                                                                                                                           ATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGT         <  SRR10918230.2049385/151‑1 (MQ=60)
                                                                                                                                              CTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGA      >  SRR10918230.605328/1‑151 (MQ=60)
                                                                                                                                                 GACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGATCA   <  SRR10918230.3747412/151‑1 (MQ=60)
                                                                                                                                                 GACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGATCA   >  SRR10918230.3414125/1‑151 (MQ=60)
                                                                                                                                                  ACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGATCAG  <  SRR10918230.2701946/151‑1 (MQ=60)
                                                                                                                                                      |                                                                                                                                                  
ACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTCACCAGCCAGATGCCGGTGGACAAATGGCATGCGCTGATCGGCGATCCGACCTTGGGCGACGCGATCCTCGACCGTCTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAGTGGACGTGGATCAG  >  CP009974/4957083‑4957379

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