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 F368 I0 R1 1472 214.4 7116750 95.9% 6824963 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,886,725 G→T 100% I256I (ATC→ATA yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,886,7250GT96.3% 569.4 / 10.6 187I256I (ATC→ATAyoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base G (5/2);  new base T (92/88);  total (97/90)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 187 total aligned reads displayed.

GAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCT  >  CP009273/1886580‑1886869
                                                                                                                                                 |                                                                                                                                                
gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTAt                                                                                                                                                <  1:2977574/148‑1 (MQ=255)
    ttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTg                                                                                                                                            <  1:3214340/148‑1 (MQ=255)
    ttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTg                                                                                                                                            <  1:2516876/148‑1 (MQ=255)
          aTTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCt                                                                                                                                      >  2:2614552/1‑148 (MQ=255)
          aTTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCt                                                                                                                                      >  2:2614562/1‑148 (MQ=255)
            tAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCttt                                                                                                                                    >  1:2594948/1‑148 (MQ=255)
            tAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCttt                                                                                                                                    >  1:881881/1‑148 (MQ=255)
            tAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCttt                                                                                                                                    <  1:3010270/148‑1 (MQ=255)
            tAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCttt                                                                                                                                    <  1:1237894/148‑1 (MQ=255)
              gCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTg                                                                                                                                  <  1:506453/148‑1 (MQ=255)
               cTAACAGCTCACGCAGGTTACCGCGATCACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTGTGc                                                                                                                                 <  2:2143358/148‑1 (MQ=255)
                 aaCAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCa                                                                                                                               <  1:2560059/148‑1 (MQ=255)
                  aCAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAg                                                                                                                              >  1:292649/1‑148 (MQ=255)
                  aCAGCTCACGCAGGGTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAg                                                                                                                              <  2:3338644/148‑1 (MQ=255)
                   cAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAgg                                                                                                                             <  1:2866272/148‑1 (MQ=255)
                   cAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAgg                                                                                                                             <  2:1826237/148‑1 (MQ=255)
                    aGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGt                                                                                                                            >  1:906838/1‑148 (MQ=255)
                      cTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCg                                                                                                                          >  1:3369223/1‑148 (MQ=255)
                      cTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCg                                                                                                                          >  2:1247535/1‑148 (MQ=255)
                            cAGGTTACCGCGATAACCTGGCTCACCGAGTTGCATACGAACATCCTGCGCACTCAGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCTGTGCGCTAGCCGATGGTGATGTCTTTTGCCAGGTCGAGCAGt                                                                                                                    >  2:3520780/1‑148 (MQ=255)
                            cAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGt                                                                                                                    >  1:1893290/1‑148 (MQ=255)
                            cAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGt                                                                                                                    >  1:1127593/1‑148 (MQ=255)
                             aGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGtt                                                                                                                   >  1:2133169/1‑148 (MQ=255)
                              ggTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTg                                                                                                                  <  1:780301/148‑1 (MQ=255)
                              ggTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTg                                                                                                                  <  1:2842394/148‑1 (MQ=255)
                               gTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGt                                                                                                                 >  2:1542250/1‑148 (MQ=255)
                                  aCCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGa                                                                                                              <  1:3490605/148‑1 (MQ=255)
                                    cgcgATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACt                                                                                                            <  2:1876996/148‑1 (MQ=255)
                                     gcgATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTg                                                                                                           <  1:3384554/148‑1 (MQ=255)
                                        ttAACCTGGCTCAGCGAGTTGCAGACGAATATCCTGCGCACTCTGGGCAAGACTATCAGCGCACTTTTGTTACTGCTGGGTGTCTTTTAATTCGGTGCTGTATGCTATGGTGTTGTCTTTTGCAAGGTTGAGCAGTTGttgagttgag                                                                                                        <  2:2942008/147‑1 (MQ=255)
                                        aTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgag                                                                                                        <  1:3309500/148‑1 (MQ=255)
                                        aTAACCTGGCTCACCGAGGTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgag                                                                                                        <  2:2011058/148‑1 (MQ=255)
                                         tAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgaga                                                                                                       >  1:2415430/1‑148 (MQ=255)
                                           aCCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGt                                                                                                     >  2:3030520/1‑148 (MQ=255)
                                            ccTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtg                                                                                                    <  1:2546804/148‑1 (MQ=255)
                                            ccTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtg                                                                                                    >  1:3079724/1‑148 (MQ=255)
                                             cTGGCTCACCGAGTTGCAGACGAAACTCCTGCGCACTCTGGGCAAGAAGATCAGCGCACTTTTGTAACTGCTGGGTGTATTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtga                                                                                                   <  2:257290/148‑1 (MQ=255)
                                             cTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtga                                                                                                   <  2:2859980/148‑1 (MQ=255)
                                             cTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtga                                                                                                   >  2:237516/1‑148 (MQ=255)
                                              tGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgac                                                                                                  <  2:1936122/148‑1 (MQ=255)
                                                gCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgactg                                                                                                >  2:2775691/1‑148 (MQ=255)
                                                  tCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgactgac                                                                                              <  2:3238229/148‑1 (MQ=255)
                                                         gTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATa                                                                                       >  2:2080602/1‑148 (MQ=255)
                                                         gTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATa                                                                                       <  1:1486468/148‑1 (MQ=255)
                                                          ttGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATAc                                                                                      >  1:2382388/1‑148 (MQ=255)
                                                             cAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTTCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActgg                                                                                   >  1:2693967/1‑148 (MQ=255)
                                                              aGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggc                                                                                  <  1:272814/148‑1 (MQ=255)
                                                              aGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggc                                                                                  <  2:2127983/148‑1 (MQ=255)
                                                               gACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggct                                                                                 >  1:2967524/1‑148 (MQ=255)
                                                                aCGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggctg                                                                                >  2:428399/1‑148 (MQ=255)
                                                                aCGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggctg                                                                                >  1:1349229/1‑148 (MQ=255)
                                                                aCGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggctg                                                                                >  1:1452651/1‑148 (MQ=255)
                                                                 cGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggctgg                                                                               <  2:2415430/148‑1 (MQ=255)
                                                                  gAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACtggctggc                                                                              >  2:768373/1‑148 (MQ=255)
                                                                       tCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGt                                                                         <  2:697296/148‑1 (MQ=255)
                                                                       tCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTGTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGt                                                                         <  2:2413030/148‑1 (MQ=255)
                                                                        ccTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTc                                                                        >  2:626935/1‑148 (MQ=255)
                                                                          tGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCg                                                                      >  1:1670499/1‑148 (MQ=255)
                                                                             gcACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTa                                                                   <  1:2043881/148‑1 (MQ=255)
                                                                                 tctGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCtg                                                               >  2:987691/1‑148 (MQ=255)
                                                                                    gggCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCtggtg                                                            >  1:1611028/1‑148 (MQ=255)
                                                                                       cAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGc                                                         <  1:1046580/148‑1 (MQ=255)
                                                                                       cAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGc                                                         <  1:1047685/148‑1 (MQ=255)
                                                                                             gATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTtcgtc                                                   <  2:2672825/148‑1 (MQ=255)
                                                                                               tCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGa                                                 <  2:3476614/148‑1 (MQ=255)
                                                                                                   cgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGAt                                             <  1:595184/148‑1 (MQ=255)
                                                                                                   cgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGAt                                             <  2:2426360/148‑1 (MQ=255)
                                                                                                   cgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGAt                                             >  1:2454186/1‑148 (MQ=255)
                                                                                                     cACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCa                                           >  1:2600052/1‑148 (MQ=255)
                                                                                                     cACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCa                                           >  2:416343/1‑148 (MQ=255)
                                                                                                       cTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATg                                         >  2:1171103/1‑148 (MQ=255)
                                                                                                        ttttGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGa                                        <  2:1639247/148‑1 (MQ=255)
                                                                                                           tGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGt                                     <  1:636273/148‑1 (MQ=255)
                                                                                                           tGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGt                                     <  2:906838/148‑1 (MQ=255)
                                                                                                           tGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGt                                     <  2:3012654/148‑1 (MQ=255)
                                                                                                             tAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTcc                                   >  1:1214699/1‑148 (MQ=255)
                                                                                                               aCTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGc                                 >  1:223575/1‑148 (MQ=255)
                                                                                                               aCTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGc                                 <  1:2393364/148‑1 (MQ=255)
                                                                                                               aCTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGc                                 <  1:2393457/148‑1 (MQ=255)
                                                                                                                ctgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCt                                >  2:18857/1‑148 (MQ=255)
                                                                                                                  gctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTc                              >  2:104773/1‑148 (MQ=255)
                                                                                                                      ggTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCggg                          >  2:2071516/1‑148 (MQ=255)
                                                                                                                       gtgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGa                         <  2:1859296/148‑1 (MQ=255)
                                                                                                                       gtgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGa                         <  2:2577961/148‑1 (MQ=255)
                                                                                                                        tgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGAt                        >  2:2901457/1‑148 (MQ=255)
                                                                                                                         gtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATc                       >  2:157188/1‑148 (MQ=255)
                                                                                                                         gtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATc                       >  2:2288638/1‑148 (MQ=255)
                                                                                                                         gtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATc                       >  2:1377946/1‑148 (MQ=255)
                                                                                                                         gtCTTTAAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATc                       >  2:2788302/1‑148 (MQ=255)
                                                                                                                            ttttAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGt                    <  2:3416387/148‑1 (MQ=255)
                                                                                                                                aaTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGc                <  2:2515726/148‑1 (MQ=255)
                                                                                                                                 aTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGcc               >  1:2654287/1‑148 (MQ=255)
                                                                                                                                 aTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGcc               >  2:2371915/1‑148 (MQ=255)
                                                                                                                                 aTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGcc               >  2:1492818/1‑148 (MQ=255)
                                                                                                                                   tCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATGCTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCaa             >  1:2635965/1‑148 (MQ=255)
                                                                                                                                     ggTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAAt           <  1:2425598/148‑1 (MQ=255)
                                                                                                                                       tGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATcc         >  2:2734673/1‑148 (MQ=255)
                                                                                                                                           gTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCCctc     >  2:56866/1‑148 (MQ=255)
                                                                                                                                            tAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActct    >  1:375177/1‑148 (MQ=255)
                                                                                                                                              ggCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActctct  <  2:2197572/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCT  >  CP009273/1886580‑1886869

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAAC  >  CP009273/1886576‑1886874
                                                                                                                                                     |                                                                                                                                                     
CGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCT                                                                                                                                                       >  SRR19973067.2916601/1‑150 (MQ=60)
CGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAAGTCGGTGCGGTAGGCT                                                                                                                                                       >  SRR19973067.2220848/1‑150 (MQ=60)
 GCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTA                                                                                                                                                      >  SRR19973067.2594762/1‑150 (MQ=60)
 GCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTA                                                                                                                                                      >  SRR19973067.2594370/1‑150 (MQ=60)
  CTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTAT                                                                                                                                                     <  SRR19973067.1140662/150‑1 (MQ=60)
    GAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGG                                                                                                                                                   <  SRR19973067.1213749/150‑1 (MQ=60)
    GAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGG                                                                                                                                                   <  SRR19973067.257760/150‑1 (MQ=60)
    GAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGG                                                                                                                                                   <  SRR19973067.2996367/150‑1 (MQ=60)
     AATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGT                                                                                                                                                  <  SRR19973067.3090374/150‑1 (MQ=60)
        TTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGAT                                                                                                                                               <  SRR19973067.3006888/150‑1 (MQ=60)
        TTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGAT                                                                                                                                               <  SRR19973067.3234589/150‑1 (MQ=60)
        TTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGAT                                                                                                                                               <  SRR19973067.3234572/150‑1 (MQ=60)
        TTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGAT                                                                                                                                               <  SRR19973067.2533385/150‑1 (MQ=60)
        TTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGAT                                                                                                                                               <  SRR19973067.2533378/150‑1 (MQ=60)
              ATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTT                                                                                                                                         >  SRR19973067.2611833/1‑150 (MQ=60)
              ATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTT                                                                                                                                         >  SRR19973067.887755/1‑150 (MQ=60)
                TAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTG                                                                                                                                       <  SRR19973067.1246203/150‑1 (MQ=60)
                TAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTG                                                                                                                                       <  SRR19973067.3029419/150‑1 (MQ=60)
                  GCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCC                                                                                                                                     <  SRR19973067.509917/150‑1 (MQ=60)
                    TAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAG                                                                                                                                   >  SRR19973067.294530/1‑150 (MQ=60)
                     AACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGG                                                                                                                                  <  SRR19973067.2576789/150‑1 (MQ=60)
                      ACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGT                                                                                                                                 >  SRR19973067.912839/1‑150 (MQ=60)
                      ACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGT                                                                                                                                 >  SRR19973067.727486/1‑150 (MQ=60)
                       CAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTC                                                                                                                                <  SRR19973067.2884629/150‑1 (MQ=60)
                       CAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTC                                                                                                                                <  SRR19973067.2704924/150‑1 (MQ=60)
                       CAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTC                                                                                                                                <  SRR19973067.868171/150‑1 (MQ=60)
                        AGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCG                                                                                                                               >  SRR19973067.3390623/1‑150 (MQ=60)
                              CGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGT                                                                                                                         >  SRR19973067.1135206/1‑150 (MQ=60)
                              CGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGT                                                                                                                         >  SRR19973067.1906211/1‑150 (MQ=60)
                               GCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTT                                                                                                                        >  SRR19973067.2147417/1‑150 (MQ=60)
                               GCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTT                                                                                                                        <  SRR19973067.628504/150‑1 (MQ=60)
                               GCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTT                                                                                                                        <  SRR19973067.2426036/150‑1 (MQ=60)
                                  GGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTC                                                                                                                     <  SRR19973067.2575499/150‑1 (MQ=60)
                                  GGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTC                                                                                                                     <  SRR19973067.2860643/150‑1 (MQ=60)
                                  GGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTC                                                                                                                     <  SRR19973067.785701/150‑1 (MQ=60)
                                      ACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACT                                                                                                                 <  SRR19973067.3512762/150‑1 (MQ=60)
                                         GCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGA                                                                                                              <  SRR19973067.3406060/150‑1 (MQ=60)
                                           GATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGA                                                                                                            >  SRR19973067.2431249/1‑150 (MQ=60)
                                            ATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAG                                                                                                           <  SRR19973067.3330393/150‑1 (MQ=60)
                                              AACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTG                                                                                                         >  SRR19973067.3099324/1‑150 (MQ=60)
                                                CCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGAC                                                                                                       <  SRR19973067.2563466/150‑1 (MQ=60)
                                                  TGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTG                                                                                                     >  SRR19973067.2663472/1‑150 (MQ=60)
                                                            AGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATAC                                                                                           >  SRR19973067.2397857/1‑150 (MQ=60)
                                                             GTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATAAT                                                                                          <  SRR19973067.1496314/150‑1 (MQ=60)
                                                               TGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTTCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGG                                                                                        >  SRR19973067.2711293/1‑150 (MQ=60)
                                                                 CAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCT                                                                                      >  SRR19973067.2985478/1‑150 (MQ=60)
                                                                 CAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCT                                                                                      >  SRR19973067.2986258/1‑150 (MQ=60)
                                                                 CAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCT                                                                                      >  SRR19973067.2986289/1‑150 (MQ=60)
                                                                 CAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCT                                                                                      >  SRR19973067.339038/1‑150 (MQ=60)
                                                                  AGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTG                                                                                     >  SRR19973067.1461769/1‑150 (MQ=60)
                                                                  AGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTG                                                                                     <  SRR19973067.274634/150‑1 (MQ=60)
                                                                  AGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTG                                                                                     >  SRR19973067.1358453/1‑150 (MQ=60)
                                                                  AGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTG                                                                                     >  SRR19973067.1462378/1‑150 (MQ=60)
                                                                  AGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTG                                                                                     >  SRR19973067.3186705/1‑150 (MQ=60)
                                                                          ATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTC                                                                             <  SRR19973067.2870703/150‑1 (MQ=60)
                                                                            CCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCG                                                                           >  SRR19973067.1682018/1‑150 (MQ=60)
                                                                                CGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAG                                                                       <  SRR19973067.1490926/150‑1 (MQ=60)
                                                                                 GCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGC                                                                      <  SRR19973067.2057546/150‑1 (MQ=60)
                                                                                   ACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTG                                                                    >  SRR19973067.1223323/1‑150 (MQ=60)
                                                                                   ACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTG                                                                    >  SRR19973067.2125954/1‑150 (MQ=60)
                                                                                      CTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTG                                                                 >  SRR19973067.1622090/1‑150 (MQ=60)
                                                                                           CAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTT                                                            <  SRR19973067.538210/150‑1 (MQ=60)
                                                                                           CAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTT                                                            <  SRR19973067.1054878/150‑1 (MQ=60)
                                                                                           CAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTT                                                            <  SRR19973067.1053770/150‑1 (MQ=60)
                                                                                                CGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCG                                                       <  SRR19973067.166180/150‑1 (MQ=60)
                                                                                                CGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCG                                                       >  SRR19973067.2455407/1‑150 (MQ=60)
                                                                                                    CAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGA                                                   >  SRR19973067.973727/1‑150 (MQ=60)
                                                                                                     AGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGAT                                                  <  SRR19973067.2055485/150‑1 (MQ=60)
                                                                                                     AGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGAT                                                  >  SRR19973067.2470318/1‑150 (MQ=60)
                                                                                                       CGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCA                                                <  SRR19973067.598961/150‑1 (MQ=60)
                                                                                                       CGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCA                                                >  SRR19973067.2616962/1‑150 (MQ=60)
                                                                                                        GCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCAT                                               <  SRR19973067.2713382/150‑1 (MQ=60)
                                                                                                               TGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCC                                        >  SRR19973067.1222742/1‑150 (MQ=60)
                                                                                                               TGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCC                                        <  SRR19973067.640551/150‑1 (MQ=60)
                                                                                                                GTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCG                                       >  SRR19973067.250699/1‑150 (MQ=60)
                                                                                                                 TAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGC                                      >  SRR19973067.225201/1‑150 (MQ=60)
                                                                                                                 TAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGC                                      >  SRR19973067.225142/1‑150 (MQ=60)
                                                                                                                   ACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTT                                    <  SRR19973067.2409035/150‑1 (MQ=60)
                                                                                                                   ACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTT                                    <  SRR19973067.2408942/150‑1 (MQ=60)
                                                                                                                   ACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTT                                    >  SRR19973067.2276454/1‑150 (MQ=60)
                                                                                                                             GTCTTTTAAGTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAG                          <  SRR19973067.2404334/150‑1 (MQ=60)
                                                                                                                                TTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTC                       <  SRR19973067.1778131/150‑1 (MQ=60)
                                                                                                                                 TTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCG                      <  SRR19973067.1114389/150‑1 (MQ=60)
                                                                                                                                 TTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCG                      <  SRR19973067.509382/150‑1 (MQ=60)
                                                                                                                                   TAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCC                    >  SRR19973067.2671386/1‑150 (MQ=60)
                                                                                                                                     ATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATGCTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAA                  >  SRR19973067.2653018/1‑150 (MQ=60)
                                                                                                                                         GGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCC              <  SRR19973067.2441512/150‑1 (MQ=60)
                                                                                                                                         GGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCC              <  SRR19973067.2440957/150‑1 (MQ=60)
                                                                                                                                             CGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTC          >  SRR19973067.3555965/1‑150 (MQ=60)
                                                                                                                                              GGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCT         >  SRR19973067.2620098/1‑150 (MQ=60)
                                                                                                                                              GGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCT         >  SRR19973067.377787/1‑150 (MQ=60)
                                                                                                                                                TAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCT       <  SRR19973067.1255954/150‑1 (MQ=60)
                                                                                                                                                TAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCT       <  SRR19973067.3510407/150‑1 (MQ=60)
                                                                                                                                                 AGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTT      >  SRR19973067.2658348/1‑150 (MQ=60)
                                                                                                                                                 AGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTT      >  SRR19973067.866510/1‑150 (MQ=60)
                                                                                                                                                 AGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTT      <  SRR19973067.284080/150‑1 (MQ=60)
                                                                                                                                                     TATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAAC  <  SRR19973067.1552711/150‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
CGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAAC  >  CP009273/1886576‑1886874

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