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 F98 I0 R1 799 294.1 9752082 92.2% 8991419 150.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,830,311 G→T 69.5% intergenic (+395/‑2) yicT → / → setC pseudogene/putative arabinose efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,830,3110GT69.5% 201.6 / 145.3 150intergenic (+395/‑2)yicT/setCpseudogene/putative arabinose efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (27/19);  new base T (50/53);  total (78/72)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

TAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTC  >  CP009273/3830165‑3830457
                                                                                                                                                  |                                                                                                                                                  
taataaTCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAAt                                                                                                                                                 >  1:2911976/1‑150 (MQ=255)
taataaTCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAAt                                                                                                                                                 >  2:571931/1‑150 (MQ=255)
taataaTCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGaataat                                                                                                                                                 >  2:38746/1‑150 (MQ=255)
taataaTCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGaataat                                                                                                                                                 >  1:2547119/1‑150 (MQ=255)
taataaTCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGaataat                                                                                                                                                 >  1:2394509/1‑150 (MQ=255)
      tctcAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCaaaa                                                                                                                                           >  1:4647017/1‑150 (MQ=255)
      tctcAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCaaaa                                                                                                                                           >  1:159618/1‑150 (MQ=255)
             tGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCt                                                                                                                                    >  2:3601108/1‑150 (MQ=255)
             tGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCt                                                                                                                                    >  1:680320/1‑150 (MQ=255)
             tGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCt                                                                                                                                    >  1:3161890/1‑150 (MQ=255)
                tGGCCCTATAATATTGAGCTTGAATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTAcc                                                                                                                                 <  1:3600303/150‑1 (MQ=255)
                tGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTAcc                                                                                                                                 >  2:1492495/1‑150 (MQ=255)
                tGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTAcc                                                                                                                                 >  2:1492353/1‑150 (MQ=255)
                 ggCCCTATAATATTGAGCTTGCTTTGCTTTCTGGTTGTATCAGCTATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCa                                                                                                                                <  2:823376/150‑1 (MQ=255)
                      tataATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCa                                                                                                                           >  1:281389/1‑150 (MQ=255)
                      tataATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCa                                                                                                                           >  1:2712333/1‑150 (MQ=255)
                         aataTTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCa                                                                                                                        >  2:1681699/1‑150 (MQ=255)
                         aataTTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCa                                                                                                                        <  2:1831930/150‑1 (MQ=255)
                         aataTTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCa                                                                                                                        <  2:2624489/150‑1 (MQ=255)
                         aataTTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCa                                                                                                                        <  2:195297/150‑1 (MQ=255)
                          atatTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCaa                                                                                                                       >  1:1183567/1‑150 (MQ=255)
                           tatTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCaaa                                                                                                                      >  1:2644621/1‑150 (MQ=255)
                            atTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCaaaa                                                                                                                     >  2:4491409/1‑150 (MQ=255)
                                      cATTGCTTTATGGTTTTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTTTAACTCATTTATATTTTTTCGCCCGGTGTTTCAGAATAATGAAAAAAACGGCTACCACTCCATCAAAAATACTTGAtc                                                                                                           <  1:4621067/150‑1 (MQ=255)
                                       aTTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGAtct                                                                                                          >  2:777121/1‑150 (MQ=255)
                                        ttGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGAtctc                                                                                                         <  2:3194093/150‑1 (MQ=255)
                                         tggtttCTGTTTTTATCAGCTATACAAAAAAATTCTTAATACAATACTCATTCGACTTTTATTTTTTTGTAAATCATTTTTTTTTTTTCGCCCGGTTATTCAGATTAATGAAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCa                                                                                                        <  1:4086800/147‑1 (MQ=255)
                                         tGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCa                                                                                                        >  2:319945/1‑150 (MQ=255)
                                           cTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACt                                                                                                      >  1:1298474/1‑150 (MQ=255)
                                                    tGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCa                                                                                             >  1:3041200/1‑150 (MQ=255)
                                                      tATCAGCGATACTCAAAAATTCTTAATACAATACCCATTCGACTGGTACTTATTTGTAACTCCGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACTGCTAAAACTCCATCAAAAATACTTGATCTCACTGCCGCGGCAtt                                                                                           <  2:3032894/150‑1 (MQ=255)
                                                        tCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCAtttt                                                                                         >  2:4855414/1‑150 (MQ=255)
                                                        tCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCAtttt                                                                                         >  2:471708/1‑150 (MQ=255)
                                                          aGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTCATTTGTAACTCAGTTCTATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCCTCAAAAATACTTGATCTCACTGCCGCGGCATTTTTa                                                                                       >  2:1717413/1‑150 (MQ=255)
                                                            cGATACTCAAAAATTCTTAATACAATACTCATTCTAATGGTACTTATTTGTAACTCAGTAATATTTTTTCGCCCGGTGATACAGAAAAATGCAAAAAACGGCTACCACTCCATCAAAAAAACGGGAAGTCAAAGCCGCGGCAAATTTaat                                                                                     >  1:1045070/1‑148 (MQ=255)
                                                               taaaaaaaaaTTCTAAATACAATACTAATTCGACTGGTAATTATTTGTAACTCATTTATA‑TTTTTCGCCGGGTTATTAAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTc                                                                                 <  1:2452210/145‑1 (MQ=255)
                                                               tACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGt                                                                                  >  1:3240498/1‑150 (MQ=255)
                                                                aCTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTc                                                                                 <  1:4278118/150‑1 (MQ=255)
                                                                aCTAAAAAATTCTTAATACAATACTCATTCGACTGGTACTTTTTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTc                                                                                 <  1:1825846/150‑1 (MQ=255)
                                                                  taaaaaaTTCTTAATACAATACTCTTTCGACTGTTACTTTTTTGTAACTCTGTTTTATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGc                                                                               <  1:957169/148‑1 (MQ=255)
                                                                    aaaaaTTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCt                                                                             >  1:3292721/1‑150 (MQ=255)
                                                                     aaaaTTCTTAATACAATACTCATTCGACTTGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCtt                                                                            <  2:3078946/150‑1 (MQ=255)
                                                                      aaaTTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCttt                                                                           <  2:3109039/150‑1 (MQ=255)
                                                                      aaaTTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCttt                                                                           <  1:208696/150‑1 (MQ=255)
                                                                       aaTTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTc                                                                          >  2:2894302/1‑150 (MQ=255)
                                                                         ttCTTAATAAAATAATCATTCGACTGGTACTTATTTTTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTg                                                                        <  1:3725617/150‑1 (MQ=255)
                                                                          tCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGa                                                                       <  2:3370487/150‑1 (MQ=255)
                                                                              aatacaatacTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACggg                                                                   >  1:223139/1‑150 (MQ=255)
                                                                                 acaatacTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTAt                                                                <  2:2421946/150‑1 (MQ=255)
                                                                                 acaatacTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTAt                                                                <  2:2421192/150‑1 (MQ=255)
                                                                                  caatacTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTAtt                                                               >  1:4471303/1‑150 (MQ=255)
                                                                                     tacTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCg                                                            <  1:743848/150‑1 (MQ=255)
                                                                                     tacTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCg                                                            <  2:3169491/150‑1 (MQ=255)
                                                                                      acTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCgg                                                           >  2:1421546/1‑150 (MQ=255)
                                                                                        tCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGgc                                                         <  1:366340/150‑1 (MQ=255)
                                                                                          aTTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGgcgc                                                       >  1:1710802/1‑150 (MQ=255)
                                                                                          aTTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGgcgc                                                       >  1:1710418/1‑150 (MQ=255)
                                                                                          aTTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGgcgc                                                       >  2:1972593/1‑150 (MQ=255)
                                                                                           ttCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAACAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGct                                                      <  2:2737794/150‑1 (MQ=255)
                                                                                            tCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGctc                                                     <  2:584864/150‑1 (MQ=255)
                                                                                             cGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGctct                                                    <  2:3761155/150‑1 (MQ=255)
                                                                                               aCTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTc                                                  <  1:1107799/150‑1 (MQ=255)
                                                                                                  ggTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGa                                               <  2:3895078/150‑1 (MQ=255)
                                                                                                    tACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACt                                             <  2:4103681/150‑1 (MQ=255)
                                                                                                    tACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACt                                             <  2:4101691/150‑1 (MQ=255)
                                                                                                    tACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACt                                             <  2:4205366/150‑1 (MQ=255)
                                                                                                     aCTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTc                                            >  2:2647625/1‑150 (MQ=255)
                                                                                                       ttATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCt                                          >  1:1539683/1‑150 (MQ=255)
                                                                                                        tATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTa                                         <  2:159618/150‑1 (MQ=255)
                                                                                                        tATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTa                                         <  2:2201035/150‑1 (MQ=255)
                                                                                                          ttttTAACTCATTTATTTTTTTTCGCCCGGTTATTCAGAAGAATGAAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTAcc                                       <  1:4074094/150‑1 (MQ=255)
                                                                                                          tttGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTAcc                                       <  1:1085280/150‑1 (MQ=255)
                                                                                                           ttGTAACTCAGTTATATTTTTTCGCCCGGTTATTCATAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTAccc                                      <  1:1437598/150‑1 (MQ=255)
                                                                                                              tAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTaa                                   >  2:2094196/1‑150 (MQ=255)
                                                                                                               aaCTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAg                                  <  1:3566576/150‑1 (MQ=255)
                                                                                                               aaCTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAg                                  <  1:3566894/150‑1 (MQ=255)
                                                                                                                aCTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGt                                 <  2:2959947/150‑1 (MQ=255)
                                                                                                                aCTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGt                                 <  1:2036652/150‑1 (MQ=255)
                                                                                                                  tCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGtat                               >  2:106169/1‑150 (MQ=255)
                                                                                                                   cAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGtata                              >  2:871381/1‑150 (MQ=255)
                                                                                                                     gTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATAtt                            <  1:4414023/150‑1 (MQ=255)
                                                                                                                      ttATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTc                           >  1:59746/1‑150 (MQ=255)
                                                                                                                       tataTTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTcc                          >  1:2479882/1‑150 (MQ=255)
                                                                                                                       tataTTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTcc                          >  2:4577845/1‑150 (MQ=255)
                                                                                                                         tatTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTc                        <  1:4159535/150‑1 (MQ=255)
                                                                                                                            tttttCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCa                     <  2:2804672/150‑1 (MQ=255)
                                                                                                                            tttttCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCa                     >  1:4047181/1‑150 (MQ=255)
                                                                                                                             ttttCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAg                    >  2:4286463/1‑150 (MQ=255)
                                                                                                                                 cGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGa                <  2:739722/150‑1 (MQ=255)
                                                                                                                                 cGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGa                >  2:3226051/1‑150 (MQ=255)
                                                                                                                                 cGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGa                >  1:2235178/1‑150 (MQ=255)
                                                                                                                                     cGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTg            <  2:4582292/150‑1 (MQ=255)
                                                                                                                                      tgTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGa           >  1:2161748/2‑150 (MQ=255)
                                                                                                                                      tgTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGa           >  1:2161740/2‑150 (MQ=255)
                                                                                                                                      ggTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGa           <  2:680320/150‑1 (MQ=255)
                                                                                                                                       gTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGaa          >  1:1554935/1‑150 (MQ=255)
                                                                                                                                             cAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCg    >  2:863192/1‑150 (MQ=255)
                                                                                                                                             cAGAAGAATGAAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCg    <  1:571931/150‑1 (MQ=255)
                                                                                                                                              agaataATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCCGATGAACTGAAAGCCCGt   >  1:59894/1‑150 (MQ=255)
                                                                                                                                               gaataaTGCAAAAAACGGCAACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTc  >  2:1592597/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTC  >  CP009273/3830165‑3830457

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTCCTA  >  CP009273/3830165‑3830460
                                                                                                                                                  |                                                                                                                                                     
TAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAAT                                                                                                                                                    >  SRR19973081.2929923/1‑150 (MQ=60)
TAATAATCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAAT                                                                                                                                                    >  SRR19973081.2562967/1‑150 (MQ=60)
TAATAATCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAAT                                                                                                                                                    >  SRR19973081.2562992/1‑150 (MQ=60)
TAATAATCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAAT                                                                                                                                                    >  SRR19973081.2409272/1‑150 (MQ=60)
    AATCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAA                                                                                                                                                >  SRR19973081.2212905/1‑150 (MQ=60)
      TCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAA                                                                                                                                              >  SRR19973081.4675229/1‑150 (MQ=60)
      TCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAA                                                                                                                                              >  SRR19973081.2093253/1‑150 (MQ=60)
      TCTCAATTGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAA                                                                                                                                              >  SRR19973081.161093/1‑150 (MQ=60)
             TGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCT                                                                                                                                       >  SRR19973081.684910/1‑150 (MQ=60)
             TGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCT                                                                                                                                       >  SRR19973081.3181157/1‑150 (MQ=60)
             TGCTGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCT                                                                                                                                       >  SRR19973081.3181672/1‑150 (MQ=60)
                TGGCCCTATAATTTTGAGCTAGCTTTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACC                                                                                                                                    <  SRR19973081.3622113/150‑1 (MQ=60)
                TGGCCCTATAATATTGAGCTTGAATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACC                                                                                                                                    <  SRR19973081.3622132/150‑1 (MQ=60)
                TGGCCCTATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACC                                                                                                                                    >  SRR19973081.4610095/1‑150 (MQ=60)
                    CCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTC                                                                                                                                >  SRR19973081.3455607/1‑150 (MQ=60)
                      TATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCA                                                                                                                              >  SRR19973081.283563/1‑150 (MQ=60)
                      TATAATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCA                                                                                                                              >  SRR19973081.2729096/1‑150 (MQ=60)
                          ATATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAA                                                                                                                          >  SRR19973081.1191111/1‑150 (MQ=60)
                           TATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAA                                                                                                                         >  SRR19973081.4162882/1‑150 (MQ=60)
                           TATTGAGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAA                                                                                                                         >  SRR19973081.2661086/1‑150 (MQ=60)
                            ATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAA                                                                                                                        >  SRR19973081.2420635/1‑150 (MQ=60)
                                      CATTGCTTTATGGTTTTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTTTAACTCATTTATATTTTTTCGCCCGGTGTTTCAGAATAATGAAAAAAACGGCTACCACTCCATCAAAAATACTTGATC                                                                                                              <  SRR19973081.4649122/150‑1 (MQ=60)
                                         ttgtttgtgtttttatctgctataataaaaatttataAATACAATACTCATTCTACTGTTATTTATTTGTAACTATTTTATTTTTTTTCGCCGGGTGATTAATAATATTTAAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCA                                                                                                           <  SRR19973081.4111557/113‑1 (MQ=34)
                                         TGGTTTCTGTTTTTATCAGCTATACAAAAAAATTCTTAATACAATACTCATTCGACTTTTATTTTTTTGTAAATCATTTTTTTTTTTTCGCCCGGTTATTCAGATTAATGAAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCA                                                                                                           <  SRR19973081.4111598/150‑1 (MQ=58)
                                           CTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACT                                                                                                         >  SRR19973081.1307139/1‑150 (MQ=60)
                                                    TGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCA                                                                                                >  SRR19973081.3059936/1‑150 (MQ=60)
                                                            CGATACTCAAAAATTCTTAATACAATACTCATTCTAATGGTACTTATTTGTAACTCAGTAATATTTTTTCGCCCGGTGATACAGAAAAATGCAAAAAACGGCTACCACTCCATCAAAAAaacgggaagtcaaagccgcggcaaatttaat                                                                                        >  SRR19973081.1052084/1‑119 (MQ=60)
                                                               taaaaAAAAATTCTAAATACAATACTAATTCGACTGGTAATTATTTGTAACTCATTTATA‑TTTTTCGCCGGGTTATTAAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTC                                                                                    <  SRR19973081.2467356/145‑1 (MQ=60)
                                                               TACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGT                                                                                     >  SRR19973081.3260007/1‑150 (MQ=60)
                                                                ACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTC                                                                                    <  SRR19973081.4304207/150‑1 (MQ=60)
                                                                ACTAAAAAATTCTTAATACAATACTCATTCGACTGGTACTTTTTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTC                                                                                    <  SRR19973081.1837737/150‑1 (MQ=60)
                                                                  TAAAAAATTCTTAATACAATACTCTTTCGACTGTTACTTTTTTGTAACTCTGTTTTATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGC                                                                                  <  SRR19973081.963661/150‑1 (MQ=60)
                                                                   CAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCC                                                                                 >  SRR19973081.667747/1‑150 (MQ=60)
                                                                    AAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCT                                                                                >  SRR19973081.3312787/1‑150 (MQ=60)
                                                                      AAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTT                                                                              <  SRR19973081.210488/150‑1 (MQ=60)
                                                                         TTCTTAATACAATACTCCTTCGACTGGTACTTTTTTGTAACTCAGTTATATTTTTTCGCCCGGTTATTCATAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTG                                                                           <  SRR19973081.3748652/150‑1 (MQ=60)
                                                                         TTCTTAATAAAATAATCATTCGACTGGTACTTATTTTTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTG                                                                           <  SRR19973081.3748384/150‑1 (MQ=60)
                                                                              AATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGG                                                                      >  SRR19973081.225017/1‑150 (MQ=60)
                                                                               ATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGT                                                                     <  SRR19973081.1193898/150‑1 (MQ=60)
                                                                                  CAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATT                                                                  >  SRR19973081.4498157/1‑150 (MQ=60)
                                                                                  CAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATT                                                                  >  SRR19973081.1726228/1‑150 (MQ=60)
                                                                                     TACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCG                                                               <  SRR19973081.748703/150‑1 (MQ=60)
                                                                                     TACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCG                                                               <  SRR19973081.748776/150‑1 (MQ=60)
                                                                                       CTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGG                                                             >  SRR19973081.2900775/1‑150 (MQ=60)
                                                                                        TCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGC                                                            <  SRR19973081.368809/150‑1 (MQ=60)
                                                                                          ATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGC                                                          >  SRR19973081.1722026/1‑150 (MQ=60)
                                                                                          ATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGC                                                          >  SRR19973081.1721640/1‑150 (MQ=60)
                                                                                          ATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGC                                                          >  SRR19973081.1721623/1‑150 (MQ=60)
                                                                                               ACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTC                                                     <  SRR19973081.1115092/150‑1 (MQ=60)
                                                                                               ACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTC                                                     <  SRR19973081.1115020/150‑1 (MQ=60)
                                                                                               ACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTC                                                     >  SRR19973081.2121327/1‑150 (MQ=60)
                                                                                                       TTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCT                                             >  SRR19973081.1479319/1‑150 (MQ=60)
                                                                                                       TTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCT                                             >  SRR19973081.1549726/1‑150 (MQ=60)
                                                                                                         ATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTAC                                           >  SRR19973081.1554539/1‑150 (MQ=60)
                                                                                                          TTTTTAACTCATTTATTTTTTTTCGCCCGGTTATTCAGAAGAATGAAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACC                                          <  SRR19973081.4098764/150‑1 (MQ=60)
                                                                                                          TTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACC                                          <  SRR19973081.1092472/150‑1 (MQ=60)
                                                                                                           TTGTAACTCAGTTATATTTTTTCGCCCGGTTATTCATAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCC                                         <  SRR19973081.1447191/150‑1 (MQ=60)
                                                                                                               AACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAG                                     <  SRR19973081.2392775/150‑1 (MQ=60)
                                                                                                               AACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAG                                     <  SRR19973081.3588635/150‑1 (MQ=60)
                                                                                                               AACTCAGTTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAG                                     <  SRR19973081.3588315/150‑1 (MQ=60)
                                                                                                                ACTCAGTTATATTTTTTCGCCCGGTGATTAAGAATAATGAAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGT                                    <  SRR19973081.2049429/150‑1 (MQ=60)
                                                                                                                   CAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATA                                 <  SRR19973081.2671995/150‑1 (MQ=60)
                                                                                                                     GTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATT                               <  SRR19973081.4439629/150‑1 (MQ=60)
                                                                                                                     GTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATT                               <  SRR19973081.4440600/150‑1 (MQ=60)
                                                                                                                      TTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTC                              >  SRR19973081.60447/1‑150 (MQ=60)
                                                                                                                      TTATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTC                              >  SRR19973081.59829/1‑150 (MQ=60)
                                                                                                                       TATATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCC                             >  SRR19973081.2495300/1‑150 (MQ=60)
                                                                                                                         TATTTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTC                           <  SRR19973081.4184950/150‑1 (MQ=60)
                                                                                                                            TTTTTCGCCCGGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCA                        >  SRR19973081.4071629/1‑150 (MQ=60)
                                                                                                                            TTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCA                        <  SRR19973081.3622941/150‑1 (MQ=60)
                                                                                                                                 CGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGA                   >  SRR19973081.2249340/1‑150 (MQ=60)
                                                                                                                                      TGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGA              >  SRR19973081.2175508/1‑150 (MQ=60)
                                                                                                                                      TGTGATTCAGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGA              >  SRR19973081.2175516/1‑150 (MQ=60)
                                                                                                                                       GTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAA             >  SRR19973081.1565026/1‑150 (MQ=60)
                                                                                                                                             CAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCG       >  SRR19973081.557698/1‑150 (MQ=60)
                                                                                                                                             CAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCG       >  SRR19973081.3003707/1‑150 (MQ=60)
                                                                                                                                             CAGAAGAATGAAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCG       <  SRR19973081.575997/150‑1 (MQ=60)
                                                                                                                                              AGAATAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCCGATGAACTGAAAGCCCGT      >  SRR19973081.60598/1‑150 (MQ=60)
                                                                                                                                                  TAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTCCTA  >  SRR19973081.2627468/1‑150 (MQ=60)
                                                                                                                                                  TAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCATCGCAGATGAACTGAAAGCCCGTCCTA  >  SRR19973081.4337747/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                                     
TAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTCCTA  >  CP009273/3830165‑3830460

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