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 F78 I0 R1 626 241.3 10389380 77.7% 8072548 149.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,386,008 G→T 79.9% Q35K (CAG→AAG)  ycjY ← S9 homolog non‑peptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,386,0080GT79.9% 524.5 / 140.0 255Q35K (CAG→AAG) ycjYS9 homolog non‑peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base T (82/120);  total (111/144)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

CTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGAT  >  CP009273/1385861‑1386154
                                                                                                                                                    |                                                                                                                                                  
cTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGAt                                                                                                                                                  <  2:4797061/150‑1 (MQ=255)
 tAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATa                                                                                                                                                 <  1:358877/150‑1 (MQ=255)
 tAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAAATGCCTGATa                                                                                                                                                 >  2:3796626/1‑150 (MQ=255)
   ccTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACt                                                                                                                                               <  2:3974122/149‑1 (MQ=255)
     tGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTTTTCTTTATCCCCCCCACCAGGGTGAGATAAAA‑CAAAAGCCTGATaaaga                                                                                                                                             >  1:3917034/1‑146 (MQ=255)
      gACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCGTTAACGCCCCCACCGGGGTGATATAACA‑CAAATGCCTGATACTGAc                                                                                                                                            >  2:2282229/1‑150 (MQ=255)
         gcgGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACggg                                                                                                                                         <  1:4529698/150‑1 (MQ=255)
          cgGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGt                                                                                                                                        <  2:1655780/150‑1 (MQ=255)
          cgGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGt                                                                                                                                        <  1:3096614/150‑1 (MQ=255)
          cgGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGACTGATACTGACGGGt                                                                                                                                        >  2:4110378/1‑150 (MQ=255)
             cTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCAGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTtc                                                                                                                                     >  2:4188448/1‑150 (MQ=255)
              tCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTCTTTAACCCCCCCCCCGGGGTGGGATAACC‑CAATTGCCTGAAACAGACGGGtttaa                                                                                                                                    >  2:1422827/1‑148 (MQ=255)
                gccgccACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTTGTGTTTAAAGCCCCCACCAGGGTGTGATAACA‑CAAATGCCAGAAACAGACGGGTTtcatc                                                                                                                                  >  1:1424994/1‑150 (MQ=255)
                       cTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAAtttt                                                                                                                           >  2:3247971/1‑150 (MQ=255)
                       cTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAAtttt                                                                                                                           <  1:3166334/150‑1 (MQ=255)
                        tttCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTg                                                                                                                          <  1:66701/150‑1 (MQ=255)
                           cACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAAAGCCCCCCCCAGGGTGGGATAACA‑CAAATGCCATATACAGACGGGGTTCATCAAATTTTggtg                                                                                                                       >  1:3079400/1‑148 (MQ=255)
                           cACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTATTTAAAGCCCCCACCAGGGTGGGATAACACCAATTG‑CTTATAAAGAAGGGTTTAATCAAAattttgtg                                                                                                                       >  1:183112/1‑146 (MQ=255)
                            aCCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAgg                                                                                                                      <  2:2091408/150‑1 (MQ=255)
                              cTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCCCCAGGGTGAGCTAACA‑CAATTGCCTTATAATGCCGGGTGTCATAAAATTTTGGAGGGa                                                                                                                    >  2:4042668/1‑150 (MQ=255)
                              cTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACCGGTTTCATCAAATTTTGGAGGGa                                                                                                                    >  2:4043193/1‑150 (MQ=255)
                              cTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCATTTACGCCCCCCCCGGGGAGCGATAAAA‑ACATAACCAGATACTGAACGGTTTAATAAAAATTTGGAcggg                                                                                                                    >  1:3952015/1‑149 (MQ=255)
                               ttGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCCCCAGGGTGAGATAACA‑CAAATGGCGTATAAAGACGGGTTTCATCAAATTTTGGAGGGaa                                                                                                                   >  1:2768558/1‑150 (MQ=255)
                               ttGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGaa                                                                                                                   <  2:4637879/150‑1 (MQ=255)
                                  aTAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCGGGGTGAGATAACA‑CAAATTCCTTATAATGACGGGTTTGATCAAATTTTGGTGGGAAATa                                                                                                                >  1:4068345/1‑150 (MQ=255)
                                     aGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTTTTCTGTAAAGCCCCCAACAGGGTGATGTAAAA‑CAAATGACATAATCTGTCGGGTTTAAGCAAATTTTGGAGGGAAATAGAt                                                                                                             >  1:3455203/1‑150 (MQ=255)
                                      gATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATc                                                                                                            <  1:4421560/150‑1 (MQ=255)
                                        tGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCAc                                                                                                          <  2:438111/150‑1 (MQ=255)
                                         gCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACt                                                                                                         <  2:2370225/150‑1 (MQ=255)
                                           gTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGGTCTTTAACGCCCCCCCCGGGGTGAGATAACA‑CCAGTGCCTTTAACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGc                                                                                                       >  2:2678309/1‑150 (MQ=255)
                                            tCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCa                                                                                                      <  1:1347061/150‑1 (MQ=255)
                                             cATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCGTATACTGACGGGTTTCATCAAATTTTTGAGGGCAATAGATCAcagcag                                                                                                     >  1:3791358/1‑150 (MQ=255)
                                             cATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCGGGGTGAGATACCA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAg                                                                                                     >  2:4955669/1‑150 (MQ=255)
                                              aTAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGa                                                                                                    <  1:2904509/150‑1 (MQ=255)
                                              aTAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGa                                                                                                    <  1:290784/150‑1 (MQ=255)
                                                 aGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACaa                                                                                                 <  1:4147649/150‑1 (MQ=255)
                                                   cAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCGTTAACGCCCCCACCAGGGTGAGATAACA‑CCATTGCCGTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCCGACAAag                                                                                               >  2:3561611/1‑150 (MQ=255)
                                                    aaTAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAaga                                                                                              >  2:2756814/1‑150 (MQ=255)
                                                     aTAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAagag                                                                                             <  2:4183313/150‑1 (MQ=255)
                                                     aTAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAagag                                                                                             >  2:4426203/1‑150 (MQ=255)
                                                       aGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGAt                                                                                           <  2:3096063/150‑1 (MQ=255)
                                                       aGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTGTAACGCCCCCCACAGGGTGTGATAAAA‑AAAATGACTTATACTGACGGGGTTAATCAAAGTTTGGATGGAAAAAGATCAATGCAGACAAAGAGAt                                                                                           >  2:969455/1‑150 (MQ=255)
                                                        gTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATAAAATTTTGGAGGGAAATAGATCACTGCAGACAAAgagaag                                                                                          >  2:2545677/1‑148 (MQ=255)
                                                            cAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTg                                                                                      <  2:2862765/150‑1 (MQ=255)
                                                               aTCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCAGGGTGAGATAAAA‑CAATTGCCTTATACTGACGGGTTTGATCCAATTTTGGTGGGGAATAGATAACGGCAGACAAAGAGGTGGTTGGat                                                                                   >  2:1485031/1‑150 (MQ=255)
                                                                 cccTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAAAGCCCCCACCGGGGTGAGATAACA‑CCATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATAACTGCAGACAAAGAGGGGGTTGGgtta                                                                                 >  2:2512203/1‑150 (MQ=255)
                                                                  ccTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGattat                                                                                <  2:389327/150‑1 (MQ=255)
                                                                  ccTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGattat                                                                                <  2:2836441/150‑1 (MQ=255)
                                                                   cTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACCCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACAGACGGGTTTCAACAAATTTTGGAGGGAAATATATCACTGCAGAAAAAGAGATGGTTGGattatt                                                                               >  2:1202088/1‑150 (MQ=255)
                                                                    ttttCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTg                                                                              <  1:2705246/150‑1 (MQ=255)
                                                                        ccGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGt                                                                          <  1:798153/150‑1 (MQ=255)
                                                                            cAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGt                                                                      <  2:445326/150‑1 (MQ=255)
                                                                             aGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACCCCCCCACCGGGGTGAGATAAAA‑CAATAGCCGTATAATGTCGGGTTTAAGAAAAATTTGGAGGGAAATATATAACTGCAGACACAGAGATTGTTGGATTATTGCTGTTAGTg                                                                     >  1:2526694/1‑150 (MQ=255)
                                                                                tttttGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAg                                                                  <  2:4413529/150‑1 (MQ=255)
                                                                                  tttGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCt                                                                <  1:5186910/150‑1 (MQ=255)
                                                                                   ttGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTg                                                               <  1:1970497/150‑1 (MQ=255)
                                                                                        aTAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATAACTGCAGACAAAGAGATTGTTGGATTATTGCTGTTAGTTAAGCTGAAttt                                                          >  2:2949591/1‑150 (MQ=255)
                                                                                               ccGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtc                                                   <  1:2075953/150‑1 (MQ=255)
                                                                                                  gCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatc                                                <  2:2216776/150‑1 (MQ=255)
                                                                                                   ccGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatca                                               >  1:574918/1‑150 (MQ=255)
                                                                                                     gTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCAtt                                             >  2:1809437/1‑150 (MQ=255)
                                                                                                     gTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCAtt                                             <  2:1586418/150‑1 (MQ=255)
                                                                                                     gTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCAtt                                             >  1:3785476/1‑150 (MQ=255)
                                                                                                      tttGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTc                                            <  1:4575198/150‑1 (MQ=255)
                                                                                                       ttGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTcc                                           >  1:5000622/1‑150 (MQ=255)
                                                                                                       ttGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTcc                                           >  2:3698406/1‑150 (MQ=255)
                                                                                                        tGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCg                                          >  1:1662685/1‑150 (MQ=255)
                                                                                                         gTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGt                                         <  2:1868852/150‑1 (MQ=255)
                                                                                                         gTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGt                                         <  2:1821898/150‑1 (MQ=255)
                                                                                                            cTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTcc                                      <  1:4426203/150‑1 (MQ=255)
                                                                                                            cTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTcc                                      <  2:849937/150‑1 (MQ=255)
                                                                                                             tttAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCt                                     >  1:3691917/1‑150 (MQ=255)
                                                                                                              ttAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTc                                    <  2:3213927/150‑1 (MQ=255)
                                                                                                              ttAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTc                                    <  2:555065/150‑1 (MQ=255)
                                                                                                                  cGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATgg                                >  1:5087948/1‑150 (MQ=255)
                                                                                                                      cccACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGtt                            >  1:975621/1‑150 (MQ=255)
                                                                                                                      cccACCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGtt                            >  1:2227186/1‑150 (MQ=255)
                                                                                                                         aCCTGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTc                         <  1:2987671/150‑1 (MQ=255)
                                                                                                                          ccTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCg                        <  2:3501643/150‑1 (MQ=255)
                                                                                                                            tGGGTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACCGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGt                      <  2:2147081/150‑1 (MQ=255)
                                                                                                                             gggTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGtt                     <  2:676638/150‑1 (MQ=255)
                                                                                                                             gggTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGtt                     <  1:10179/150‑1 (MQ=255)
                                                                                                                             gggTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGtt                     <  2:131633/150‑1 (MQ=255)
                                                                                                                             gggTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGtt                     <  1:3828818/150‑1 (MQ=255)
                                                                                                                              ggTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTc                    <  1:4337416/150‑1 (MQ=255)
                                                                                                                               gTGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCg                   <  1:827465/150‑1 (MQ=255)
                                                                                                                                tGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGt                  <  2:3750623/150‑1 (MQ=255)
                                                                                                                                tGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGt                  >  1:95702/1‑150 (MQ=255)
                                                                                                                                tGAGATAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGt                  <  2:4287397/150‑1 (MQ=255)
                                                                                                                                   gaTAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGtttt               >  1:4874423/1‑150 (MQ=255)
                                                                                                                                   gaTAACA‑CAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGAGGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGtttt               <  2:4221755/150‑1 (MQ=255)
                                                                                                                                             aTTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTTTTAGAGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGGTGTAGGTTCGTTTTAACGGTTgg      >  2:2166513/1‑150 (MQ=255)
                                                                                                                                             aTTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTgg      <  1:2414945/150‑1 (MQ=255)
                                                                                                                                              ttGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGt     >  1:4766645/1‑150 (MQ=255)
                                                                                                                                               tGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTg    >  2:5142294/1‑150 (MQ=255)
                                                                                                                                                gCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGa   >  1:2608458/1‑150 (MQ=255)
                                                                                                                                                 ccTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGat  >  2:1859325/1‑150 (MQ=255)
                                                                                                                                                 ccTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGat  >  1:1556535/1‑150 (MQ=255)
                                                                                                                                                 ccTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGat  <  2:1452622/150‑1 (MQ=255)
                                                                                                                                                 ccTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCAGTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGat  <  2:4686421/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACA‑CAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGAT  >  CP009273/1385861‑1386154

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ 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
TCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGATATC  >  CP009273/1385860‑1386157
                                                                                                                                                    |                                                                                                                                                       
TCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTA                                                                                                                                                        <  SRR19973083.4987118/150‑1 (MQ=60)
TCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTA                                                                                                                                                        <  SRR19973083.4988191/150‑1 (MQ=60)
  TAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA                                                                                                                                                      <  SRR19973083.360575/150‑1 (MQ=60)
       GACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGAC                                                                                                                                               <  SRR19973083.4049751/150‑1 (MQ=60)
          GCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGG                                                                                                                                            <  SRR19973083.2599515/150‑1 (MQ=60)
          GCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGG                                                                                                                                            <  SRR19973083.5185648/150‑1 (MQ=60)
          GCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGG                                                                                                                                            <  SRR19973083.4549346/150‑1 (MQ=60)
           CGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGT                                                                                                                                           <  SRR19973083.3111278/150‑1 (MQ=60)
           CGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGT                                                                                                                                           <  SRR19973083.3111453/150‑1 (MQ=60)
            GGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTT                                                                                                                                          <  SRR19973083.941486/150‑1 (MQ=60)
             GCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTT                                                                                                                                         <  SRR19973083.978913/150‑1 (MQ=60)
               TCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCA                                                                                                                                       <  SRR19973083.267578/150‑1 (MQ=60)
                CGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCAT                                                                                                                                      <  SRR19973083.3528119/150‑1 (MQ=60)
                 GCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTTGTGTTTAAAGCCCCCACCAGGGTGTGATAACACAAATGCCAGAAA‑‑CAGACGGGTTTCATC                                                                                                                                     >  SRR19973083.1432137/1‑150 (MQ=60)
                        CTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTT                                                                                                                              <  SRR19973083.3181655/150‑1 (MQ=60)
                        CTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTT                                                                                                                              <  SRR19973083.3181428/150‑1 (MQ=60)
                        CTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCGGGGTGAGATAAAACAAATGCCATATA‑‑ATGACGGGTTTCATCAAAATTT                                                                                                                              >  SRR19973083.1665852/1‑150 (MQ=60)
                        CTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTCTTTAACGCCCCCACCAGGGTGAGATAACACAAAAGCCTTATA‑‑CTGACGGGTTTAATCAAATTTT                                                                                                                              >  SRR19973083.3646373/1‑150 (MQ=60)
                         TTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTG                                                                                                                             <  SRR19973083.67192/150‑1 (MQ=60)
                            CACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAG                                                                                                                          <  SRR19973083.649971/150‑1 (MQ=60)
                            CACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAAAGCCCCCCCCAGGGTGGGATAACACAAATGCCATATA‑‑CAGACGGGGTTCATCAAATTTTGGTG                                                                                                                          >  SRR19973083.3094159/1‑150 (MQ=60)
                             AACTTGATAAGATGCGTCATAAGCAATATTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGTGG                                                                                                                         <  SRR19973083.4519466/150‑1 (MQ=60)
                                TTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCCCCAGGGTGAGATAACACAAATGGCGTATA‑‑AAGACGGGTTTCATCAAATTTTGGAGGGAA                                                                                                                      >  SRR19973083.2782324/1‑150 (MQ=60)
                                   ATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCGGGGTGAGATAACACAAATTCCTTATA‑‑ATGACGGGTTTGATCAAATTTTGGTGGGAAATA                                                                                                                   >  SRR19973083.4086401/1‑150 (MQ=60)
                                     AAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTGTTTAAAGCCCCCCCCGGGGGGAGATAACACAAATGCCATATA‑‑ATTACGGGTTTCATAAAATTTTTGAGGGAAAAAGA                                                                                                                 >  SRR19973083.132464/1‑150 (MQ=60)
                                      AGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTTTTCTGTAAAGCCCCCAACAGGGTGATGTAAAACAAATGAC‑‑ATAATCTGTCGGGTTTAAGCAAATTTTGGAGGGAAATAGAT                                                                                                                >  SRR19973083.3471046/1‑150 (MQ=60)
                                       GATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATC                                                                                                               <  SRR19973083.4440806/150‑1 (MQ=60)
                                         TGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCGGGGTGAGATAACACAATTGCCGTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGTAATAGATCAC                                                                                                             >  SRR19973083.2974108/1‑150 (MQ=60)
                                            GTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGC                                                                                                          <  SRR19973083.2205513/150‑1 (MQ=60)
                                             TCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCA                                                                                                         <  SRR19973083.1353676/150‑1 (MQ=60)
                                             TCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCA                                                                                                         <  SRR19973083.1354111/150‑1 (MQ=60)
                                              CATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCGTATA‑‑CTGACGGGTTTCATCAAATTTTTGAGGGCAATAGATCACAGCAG                                                                                                        >  SRR19973083.3808448/1‑150 (MQ=60)
                                               ATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGA                                                                                                       <  SRR19973083.2918736/150‑1 (MQ=60)
                                               ATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGA                                                                                                       <  SRR19973083.2918743/150‑1 (MQ=60)
                                               ATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGA                                                                                                       <  SRR19973083.292210/150‑1 (MQ=60)
                                               ATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGA                                                                                                       <  SRR19973083.292289/150‑1 (MQ=60)
                                                TAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGAC                                                                                                      <  SRR19973083.19074/150‑1 (MQ=60)
                                                  AGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAA                                                                                                    <  SRR19973083.4165968/150‑1 (MQ=60)
                                                                    CTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGTTCACTGCAGACAAAGAGATGGTTGGATTATT                                                                                  <  SRR19973083.3304473/150‑1 (MQ=60)
                                                                    CTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATT                                                                                  <  SRR19973083.4113392/150‑1 (MQ=60)
                                                                     TTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTG                                                                                 <  SRR19973083.2718785/150‑1 (MQ=60)
                                                                     TTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTG                                                                                 <  SRR19973083.2311706/150‑1 (MQ=60)
                                                                     TTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTG                                                                                 <  SRR19973083.2719413/150‑1 (MQ=60)
                                                                       TTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCT                                                                               <  SRR19973083.1408647/150‑1 (MQ=60)
                                                                       TTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCT                                                                               <  SRR19973083.3823168/150‑1 (MQ=60)
                                                                       TTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CAGAAGGGTTTAATCAAATTTTTGAGGGAAATATATCACAGCAGACAAAGAGATGGgtgggttatttgg                                                                               >  SRR19973083.4708969/1‑137 (MQ=60)
                                                                         CCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGT                                                                             <  SRR19973083.801171/150‑1 (MQ=60)
                                                                         CCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGT                                                                             <  SRR19973083.801994/150‑1 (MQ=60)
                                                                         CCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGT                                                                             <  SRR19973083.801685/150‑1 (MQ=60)
                                                                         CCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTATTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGT                                                                             <  SRR19973083.3864519/150‑1 (MQ=60)
                                                                         CCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTTTTATTTAACGCCCCCACCGGGGTGTGATAACACAAATGCCATATA‑‑AAGACGGGTTTAATCAAAATTTTGAGGGAAATAGAAAACAGCAGAAAAAGATATGGTTGGATTATTTCTGT                                                                             >  SRR19973083.3854646/1‑150 (MQ=45)
                                                                           GCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTA                                                                           <  SRR19973083.3088229/150‑1 (MQ=60)
                                                                              AGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACCCCCCCACCGGGGTGAGATAAAACAATAGCCGTATA‑‑ATGTCGGGTTTAAGAAAAATTTGGAGGGAAATATATAACTGCAGACACAGAGATTGTTGGATTATTGCTGTTAGTG                                                                        >  SRR19973083.2539520/1‑150 (MQ=60)
                                                                                   TTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCT                                                                   <  SRR19973083.5208875/150‑1 (MQ=60)
                                                                                    TTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTG                                                                  <  SRR19973083.1980570/150‑1 (MQ=60)
                                                                                    TTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTG                                                                  <  SRR19973083.1982865/150‑1 (MQ=60)
                                                                                     TGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGA                                                                 <  SRR19973083.4066642/150‑1 (MQ=60)
                                                                                                CCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTC                                                      <  SRR19973083.2086508/150‑1 (MQ=60)
                                                                                                  GGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCAT                                                    <  SRR19973083.2102590/150‑1 (MQ=60)
                                                                                                  GGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCAT                                                    <  SRR19973083.5313/150‑1 (MQ=60)
                                                                                                  GGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑ATGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTTCTGTTATTGAAGCTGACTTTATTATTCAT                                                    >  SRR19973083.805552/1‑150 (MQ=60)
                                                                                                    CCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCA                                                  >  SRR19973083.577798/1‑150 (MQ=60)
                                                                                                      GTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATT                                                >  SRR19973083.4828547/1‑150 (MQ=60)
                                                                                                      GTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATT                                                >  SRR19973083.3802553/1‑150 (MQ=60)
                                                                                                       TTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTC                                               <  SRR19973083.4594989/150‑1 (MQ=60)
                                                                                                        TTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCC                                              >  SRR19973083.2432364/1‑150 (MQ=60)
                                                                                                        TTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCC                                              >  SRR19973083.5022088/1‑150 (MQ=60)
                                                                                                        TTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGTTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCC                                              <  SRR19973083.3652285/150‑1 (MQ=60)
                                                                                                        TTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCC                                              <  SRR19973083.3133222/150‑1 (MQ=60)
                                                                                                        TTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCATGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCC                                              <  SRR19973083.3133555/150‑1 (MQ=60)
                                                                                                         TGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCG                                             >  SRR19973083.1670930/1‑150 (MQ=60)
                                                                                                             CTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCC                                         <  SRR19973083.41921/150‑1 (MQ=60)
                                                                                                             CTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCC                                         <  SRR19973083.4445462/150‑1 (MQ=60)
                                                                                                             CTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCC                                         <  SRR19973083.4445168/150‑1 (MQ=60)
                                                                                                              TTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCT                                        >  SRR19973083.3708742/1‑150 (MQ=60)
                                                                                                              TTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCT                                        >  SRR19973083.3089836/1‑150 (MQ=60)
                                                                                                                TAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCA                                      >  SRR19973083.286333/1‑150 (MQ=60)
                                                                                                                   CGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGG                                   >  SRR19973083.5109682/1‑150 (MQ=60)
                                                                                                                   CGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGG                                   >  SRR19973083.5109413/1‑150 (MQ=60)
                                                                                                                      CCCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGT                                >  SRR19973083.4555641/1‑150 (MQ=60)
                                                                                                                       CCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTT                               >  SRR19973083.980455/1‑150 (MQ=60)
                                                                                                                       CCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTT                               <  SRR19973083.439207/150‑1 (MQ=60)
                                                                                                                       CCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTT                               <  SRR19973083.1505235/150‑1 (MQ=60)
                                                                                                                       CCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTT                               >  SRR19973083.2238465/1‑150 (MQ=60)
                                                                                                                       CCCACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTT                               <  SRR19973083.439172/150‑1 (MQ=60)
                                                                                                                       CCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTT                               >  SRR19973083.469150/1‑150 (MQ=60)
                                                                                                                          ACCTGGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTC                            <  SRR19973083.3002144/150‑1 (MQ=60)
                                                                                                                              GGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTT                        <  SRR19973083.10260/150‑1 (MQ=60)
                                                                                                                              GGGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTT                        <  SRR19973083.3846010/150‑1 (MQ=60)
                                                                                                                               GGTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTC                       <  SRR19973083.4356222/150‑1 (MQ=60)
                                                                                                                                GTGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCG                      <  SRR19973083.831391/150‑1 (MQ=60)
                                                                                                                                 TGAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGT                     >  SRR19973083.96344/1‑150 (MQ=60)
                                                                                                                                  GAGATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTT                    >  SRR19973083.3762438/1‑150 (MQ=60)
                                                                                                                                   AGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTT                   >  SRR19973083.4973268/1‑150 (MQ=60)
                                                                                                                                    GATAACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTT                  >  SRR19973083.4895355/1‑150 (MQ=60)
                                                                                                                                       AACACAATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAAC               >  SRR19973083.3803310/1‑150 (MQ=60)
                                                                                                                                            AATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTAGGTTCGTTTTAACGGTTG          >  SRR19973083.1755451/1‑150 (MQ=60)
                                                                                                                                             ATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGG         <  SRR19973083.2427346/150‑1 (MQ=60)
                                                                                                                                             ATTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGG         >  SRR19973083.4912271/1‑150 (MQ=60)
                                                                                                                                              TTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGT        <  SRR19973083.4206883/150‑1 (MQ=60)
                                                                                                                                              TTGCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGT        >  SRR19973083.4786953/1‑150 (MQ=60)
                                                                                                                                                GCCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGA      >  SRR19973083.2558933/1‑150 (MQ=60)
                                                                                                                                                GCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGA      >  SRR19973083.2621582/1‑150 (MQ=60)
                                                                                                                                                GCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGA      <  SRR19973083.1871928/150‑1 (MQ=60)
                                                                                                                                                 CCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGAT     <  SRR19973083.2757201/150‑1 (MQ=60)
                                                                                                                                                 CCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGAT     >  SRR19973083.1565072/1‑150 (MQ=60)
                                                                                                                                                 CCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGAT     >  SRR19973083.1564356/1‑150 (MQ=60)
                                                                                                                                                 CCTTATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGAT     <  SRR19973083.3804936/150‑1 (MQ=60)
                                                                                                                                                    TATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGATATC  <  SRR19973083.3889029/150‑1 (MQ=60)
                                                                                                                                                    |                                                                                                                                                       
TCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATA‑‑CTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGAGTTGTCGGTTCGTTTTAACGGTTGGTGATATC  >  CP009273/1385860‑1386157

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