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
A4 F1 I1 R1 256 520.0 22514234 99.3% 22356634 151.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,947,645 G→C Q153H (CAG→CAC)  PA14_RS22670 → DNA polymerase Y family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,947,6450GC99.7% 1387.9 / NA 356Q153H (CAG→CAC) PA14_RS22670DNA polymerase Y family protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (207/148);  total (208/148)
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 356 total aligned reads displayed.

CCCTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAAGGT‑CAT‑GCGCACTACACAGGACAACATTTCCATCGAGGACGAGACAGCGCTGCAACACCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTGCG  >  NC_008463/4947521‑4947761
                                                                                                                              |                                                                                                                    
cccTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGAc                                                                                            >  2:2999879/1‑151 (MQ=255)
 ccTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACa                                                                                           >  2:3721859/1‑151 (MQ=255)
 ccTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACa                                                                                           <  2:10308446/151‑1 (MQ=255)
 ccTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCGGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACa                                                                                           <  2:8254640/151‑1 (MQ=255)
  cTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAg                                                                                          >  1:49202/1‑151 (MQ=255)
    gTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAgcg                                                                                        >  1:2902996/1‑151 (MQ=255)
     ttCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAgcgc                                                                                       >  2:10412285/1‑151 (MQ=255)
        ggCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCGG‑GCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGc                                                                                    <  2:4977329/151‑1 (MQ=255)
         gCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCa                                                                                   >  2:9135472/1‑151 (MQ=255)
          ccccTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCaa                                                                                  >  2:8581071/1‑151 (MQ=255)
          ccccTGGCCACGCTTCGAGGCCCGCGTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCac                                                                                  <  1:1154569/151‑2 (MQ=255)
           cccTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAac                                                                                 <  2:7403807/151‑1 (MQ=255)
            ccTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGACCTCACGGCGCTGGGCTTCCGCCATCGCCGGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGCCGGCCTGGCCATCGAGGACGAGACAGCGCTGCaccg                                                                                >  2:11074807/1‑148 (MQ=255)
            ccTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAacg                                                                                <  1:138788/151‑2 (MQ=255)
              tGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCAG‑GCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGcc                                                                              >  1:7697867/1‑151 (MQ=255)
              tGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGcc                                                                              >  1:3318310/1‑151 (MQ=255)
               ggCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACCGCGCTGCAACGCCa                                                                             >  2:3915123/1‑151 (MQ=255)
               ggCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCa                                                                             >  1:2117446/1‑151 (MQ=255)
                 ccACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGt                                                                           >  1:10762652/1‑151 (MQ=255)
                 ccACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGt                                                                           <  1:6791528/151‑1 (MQ=255)
                   aCGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGACGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTg                                                                         >  2:396826/1‑151 (MQ=255)
                     gCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTc                                                                       <  1:8128916/151‑1 (MQ=255)
                     gCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTc                                                                       <  1:6866154/151‑1 (MQ=255)
                     gCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTc                                                                       >  2:6313630/1‑151 (MQ=255)
                       ttCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCg                                                                     >  1:4454482/1‑151 (MQ=255)
                        tCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGc                                                                    >  2:2806648/1‑151 (MQ=255)
                         cGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGcc                                                                   >  2:1313219/1‑151 (MQ=255)
                         cGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGcc                                                                   <  2:4728907/151‑1 (MQ=255)
                          gAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGccc                                                                  >  2:2082823/1‑151 (MQ=255)
                          gAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGccc                                                                  <  1:7257796/151‑1 (MQ=255)
                           aGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCt                                                                 >  2:4399862/1‑151 (MQ=255)
                            ggCCCGCTTGCGCAAGGAACTCATGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCtt                                                                <  2:7941872/151‑1 (MQ=255)
                            ggCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCtt                                                                >  2:7101950/1‑151 (MQ=255)
                             gCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTc                                                               <  1:2606839/151‑1 (MQ=255)
                             gCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTc                                                               >  2:5006448/1‑151 (MQ=255)
                              cccGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTcc                                                              >  2:10868192/1‑151 (MQ=255)
                               ccGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCg                                                             >  1:9587581/1‑151 (MQ=255)
                                cGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGNCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGc                                                            >  1:5880742/1‑151 (MQ=255)
                                 gCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCt                                                           <  2:4094736/151‑1 (MQ=255)
                                 gCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCt                                                           >  1:2768751/1‑151 (MQ=255)
                                  cTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCCTCGCCTGGGCCTGCCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGCCGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTg                                                          <  1:10539793/151‑1 (MQ=255)
                                    tGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGa                                                        >  1:1602945/1‑151 (MQ=255)
                                      cgcAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGAcc                                                      <  2:2912591/151‑1 (MQ=255)
                                       gcAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCg                                                     >  2:6055707/1‑151 (MQ=255)
                                       gcAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCg                                                     >  2:7000748/1‑151 (MQ=255)
                                        cAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGc                                                    >  2:8640448/1‑151 (MQ=255)
                                         aaGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGcc                                                   >  2:5520300/1‑151 (MQ=255)
                                         aaGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGcc                                                   >  1:7107583/1‑151 (MQ=255)
                                          aGGCACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCt                                                  <  1:7936716/151‑1 (MQ=255)
                                          aGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCt                                                  <  1:6645761/151‑1 (MQ=255)
                                            gAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTgg                                                <  1:8372826/151‑1 (MQ=255)
                                            gAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTgg                                                >  2:6710795/1‑151 (MQ=255)
                                                tCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGAct                                            <  2:7820619/151‑1 (MQ=255)
                                                 cACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGActc                                           >  1:6228788/1‑151 (MQ=255)
                                                 cACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑ACCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGActc                                           <  1:439023/151‑1 (MQ=255)
                                                  aCGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGActct                                          >  2:1148264/1‑151 (MQ=255)
                                                   cGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACtctc                                         >  1:4798616/1‑151 (MQ=255)
                                                   cGGCGCTGGGCGTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACtctc                                         <  1:6388033/151‑1 (MQ=255)
                                                    ggCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTcg                                        >  2:8085339/1‑151 (MQ=255)
                                                    ggCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTcg                                        <  1:5650036/151‑1 (MQ=255)
                                                     gcgcTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTcgc                                       >  2:2180857/1‑151 (MQ=255)
                                                      cgcTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTcgcg                                      >  1:2383796/1‑151 (MQ=255)
                                                       gcTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCgg                                     >  1:10127516/1‑151 (MQ=255)
                                                         tGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGa                                   <  2:8319326/151‑1 (MQ=255)
                                                         tGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGa                                   >  2:10605287/1‑151 (MQ=255)
                                                          gggCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAg                                  <  1:4593905/151‑1 (MQ=255)
                                                           ggCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGc                                 >  1:3198434/1‑151 (MQ=255)
                                                           ggCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGc                                 >  2:6161224/1‑151 (MQ=255)
                                                            gCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCt                                >  2:8902683/1‑151 (MQ=255)
                                                               tCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGc                             <  1:10035531/151‑1 (MQ=255)
                                                               tCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGc                             <  1:6170396/151‑1 (MQ=255)
                                                                 cGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCg                           >  1:2491125/1‑151 (MQ=255)
                                                                 cGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCg                           >  1:9793947/1‑151 (MQ=255)
                                                                 cGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCg                           >  2:4584737/1‑151 (MQ=255)
                                                                   ccATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTc                         >  2:1730263/1‑151 (MQ=255)
                                                                   ccATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTc                         <  2:7164536/151‑1 (MQ=255)
                                                                         ccTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTc                   <  2:1773975/151‑1 (MQ=255)
                                                                         ccTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTc                   <  2:9519938/151‑1 (MQ=255)
                                                                         ccTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTc                   <  1:7350156/151‑1 (MQ=255)
                                                                          cTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCg                  <  2:5717384/151‑1 (MQ=255)
                                                                          cTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCg                  <  1:5503814/151‑1 (MQ=255)
                                                                           tgtgCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCgg                 >  1:2755840/1‑151 (MQ=255)
                                                                            gtgCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGa                >  1:8540520/1‑151 (MQ=255)
                                                                            gtgCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGa                >  1:411602/1‑151 (MQ=255)
                                                                            gtgCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGa                >  2:269179/1‑151 (MQ=255)
                                                                            gtgCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGa                <  1:1031844/151‑1 (MQ=255)
                                                                             tgCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGAc               <  1:9126833/151‑1 (MQ=255)
                                                                              gCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACg              <  2:3248573/151‑1 (MQ=255)
                                                                              gCCTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCAGCAACGCCAGTTGTCCCCCCATCCGCTGGACCGCCTGGGACACTCGCGGGAGCTGGCCGTCGCCCTCGGACg              >  2:10761405/1‑151 (MQ=255)
                                                                                cTGGCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCa            >  1:8186924/1‑151 (MQ=255)
                                                                                  ggCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATg          <  1:8908407/151‑1 (MQ=255)
                                                                                  ggCGCCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATg          >  2:3632044/1‑151 (MQ=255)
                                                                                   gcgcCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATgg         >  2:10204707/1‑151 (MQ=255)
                                                                                     gcCGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGa       >  2:9831198/1‑151 (MQ=255)
                                                                                      ccGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGAc      >  1:8710655/1‑151 (MQ=255)
                                                                                      ccGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGAc      >  1:6982153/1‑151 (MQ=255)
                                                                                      ccGAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGAc      >  1:11037019/1‑151 (MQ=255)
                                                                                        gAACCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTg    >  1:3528662/1‑151 (MQ=255)
                                                                                          aCCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTgcg  >  1:787081/1‑151 (MQ=255)
                                                                                          aCCCGGTGGCGGCGAGGGT‑CCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTgcg  >  2:8869001/1‑151 (MQ=255)
                                                                                                                              |                                                                                                                    
CCCTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAAGGT‑CAT‑GCGCACTACACAGGACAACATTTCCATCGAGGACGAGACAGCGCTGCAACACCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTGCG  >  NC_008463/4947521‑4947761

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAAGGTCAT‑GCGCACTACACAGGACAACATTTCCATCGAGGACGAGACAGCGCTGCAACACCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCG  >  NC_008463/4947547‑4947739
                                                                                                   |                                                                                              
GAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCC                                            <  SRR18135074.7285933/151‑1 (MQ=60)
  GGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTT                                          <  SRR18135074.10871728/151‑1 (MQ=60)
  GGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTT                                          <  SRR18135074.465145/151‑1 (MQ=60)
   GCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTC                                         <  SRR18135074.2617428/151‑1 (MQ=60)
   GCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTC                                         <  SRR18135074.6742761/151‑1 (MQ=60)
   GCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTC                                         >  SRR18135074.1601304/1‑151 (MQ=60)
    CCCGCTTGCGCAAGGAACTCCCGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCC                                        <  SRR18135074.2448943/151‑1 (MQ=60)
     CCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCG                                       >  SRR18135074.9623542/1‑151 (MQ=60)
     CCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCG                                       >  SRR18135074.3957615/1‑151 (MQ=60)
      CGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGNCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGC                                      >  SRR18135074.5904426/1‑151 (MQ=60)
      CGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGC                                      <  SRR18135074.5868867/151‑1 (MQ=60)
      CGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGC                                      >  SRR18135074.779268/1‑151 (MQ=60)
       GCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCT                                     <  SRR18135074.10599287/151‑1 (MQ=60)
       GCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCT                                     >  SRR18135074.2779956/1‑151 (MQ=60)
       GCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCT                                     >  SRR18135074.957820/1‑151 (MQ=60)
        CTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCCTCGCCTGGGCCTGCCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGCCGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTG                                    <  SRR18135074.10578988/151‑1 (MQ=47)
        CTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTG                                    <  SRR18135074.7520580/151‑1 (MQ=60)
        CTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTG                                    >  SRR18135074.8802511/1‑151 (MQ=60)
         TTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGG                                   >  SRR18135074.8409815/1‑151 (MQ=60)
         TTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGG                                   >  SRR18135074.6869442/1‑151 (MQ=60)
          TGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGA                                  >  SRR18135074.1609590/1‑151 (MQ=60)
            CGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACC                                >  SRR18135074.3239596/1‑151 (MQ=60)
               AAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCC                             >  SRR18135074.7135274/1‑151 (MQ=60)
               AAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCC                             >  SRR18135074.7135262/1‑151 (MQ=60)
                AGGCACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCT                            <  SRR18135074.7967067/151‑1 (MQ=60)
                AGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCT                            <  SRR18135074.6671893/151‑1 (MQ=60)
                  GAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGG                          <  SRR18135074.8404598/151‑1 (MQ=60)
                   AACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGG                         >  SRR18135074.3416693/1‑151 (MQ=60)
                     CTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGAC                       >  SRR18135074.2095962/1‑151 (MQ=60)
                       CACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTC                     >  SRR18135074.6253642/1‑151 (MQ=60)
                       CACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTC                     >  SRR18135074.1628490/1‑151 (MQ=60)
                       CACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑ACCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTC                     <  SRR18135074.440963/151‑1 (MQ=60)
                        ACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCT                    >  SRR18135074.4434326/1‑151 (MQ=60)
                         CGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTC                   >  SRR18135074.10125760/1‑151 (MQ=60)
                         CGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTC                   <  SRR18135074.819617/151‑1 (MQ=60)
                         CGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTC                   >  SRR18135074.4818013/1‑151 (MQ=60)
                         CGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTC                   >  SRR18135074.4817961/1‑151 (MQ=60)
                         CGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTC                   <  SRR18135074.1425044/151‑1 (MQ=60)
                         CGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTC                   >  SRR18135074.4817962/1‑151 (MQ=60)
                         CGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTC                   >  SRR18135074.264563/1‑151 (MQ=60)
                         CGGCGCTGGGCGTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTC                   <  SRR18135074.6413348/151‑1 (MQ=60)
                          GGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCG                  >  SRR18135074.10434617/1‑151 (MQ=60)
                          GGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCG                  <  SRR18135074.5672849/151‑1 (MQ=60)
                          GGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCG                  <  SRR18135074.5672806/151‑1 (MQ=60)
                            CGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCG                >  SRR18135074.2393523/1‑151 (MQ=60)
                             GCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGG               >  SRR18135074.10165252/1‑151 (MQ=60)
                             GCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGG               <  SRR18135074.11015974/151‑1 (MQ=60)
                             GCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGG               >  SRR18135074.7186626/1‑151 (MQ=60)
                             GCTGGGCTTCCGCCAGCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGG               <  SRR18135074.10618060/151‑1 (MQ=60)
                              CTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGG              >  SRR18135074.4424794/1‑151 (MQ=60)
                                GGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAG            <  SRR18135074.4612373/151‑1 (MQ=60)
                                 GGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGC           >  SRR18135074.3211407/1‑151 (MQ=60)
                                  GCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCT          <  SRR18135074.8164068/151‑1 (MQ=60)
                                  GCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCT          >  SRR18135074.2311988/1‑151 (MQ=60)
                                   CTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTG         <  SRR18135074.1866904/151‑1 (MQ=60)
                                    TTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGG        >  SRR18135074.7250755/1‑151 (MQ=60)
                                     TCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGC       >  SRR18135074.8360156/1‑151 (MQ=60)
                                     TCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGC       <  SRR18135074.6195090/151‑1 (MQ=60)
                                     TCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGC       <  SRR18135074.10072990/151‑1 (MQ=60)
                                     TCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTAGC       >  SRR18135074.3301176/1‑151 (MQ=60)
                                       CGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCG     >  SRR18135074.2497311/1‑151 (MQ=60)
                                       CGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCG     >  SRR18135074.2501274/1‑151 (MQ=60)
                                       CGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCG     >  SRR18135074.9830543/1‑151 (MQ=60)
                                       CGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCG     >  SRR18135074.9830589/1‑151 (MQ=60)
                                       CGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCG     >  SRR18135074.9830633/1‑151 (MQ=60)
                                        GCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGT    >  SRR18135074.130219/1‑151 (MQ=60)
                                        GCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGT    >  SRR18135074.8521538/1‑151 (MQ=60)
                                        GCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACGCTCGCGGGAGCTGGCCGT    >  SRR18135074.130224/1‑151 (MQ=60)
                                          CATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAAC‑AGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCG  >  SRR18135074.9331010/1‑151 (MQ=60)
                                                                                                   |                                                                                              
GAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAAGGTCAT‑GCGCACTACACAGGACAACATTTCCATCGAGGACGAGACAGCGCTGCAACACCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCG  >  NC_008463/4947547‑4947739

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