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 F434 I0 R1 1130 215.9 7031662 97.1% 6827743 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,984,987 C→A 100% T38T (ACG→ACT)  yecA ← UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,984,9870CA99.0% 654.8 / ‑1.8 203T38T (ACG→ACT) yecAUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (106/95);  total (107/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 205 total aligned reads displayed.

TATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCA  >  CP009273/1984845‑1985130
                                                                                                                                              |                                                                                                                                               
tATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCaa                                                                                                                                            <  1:1945182/148‑1 (MQ=255)
tATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCaa                                                                                                                                            >  2:2050232/1‑148 (MQ=255)
 aTCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAAt                                                                                                                                           >  2:829046/1‑148 (MQ=255)
    ggCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATaaa                                                                                                                                        <  2:3502642/148‑1 (MQ=255)
     gCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAAc                                                                                                                                       <  1:354507/148‑1 (MQ=255)
      ccATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAAcc                                                                                                                                      >  1:1858409/1‑148 (MQ=255)
       cATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCg                                                                                                                                     <  1:3060317/148‑1 (MQ=255)
       cATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCg                                                                                                                                     <  1:1345008/148‑1 (MQ=255)
        atatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGt                                                                                                                                    <  1:609743/148‑1 (MQ=255)
         tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTc                                                                                                                                   <  2:1799933/148‑1 (MQ=255)
          atGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTcc                                                                                                                                  >  2:1876868/1‑148 (MQ=255)
           tGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCa                                                                                                                                 >  1:3144066/1‑148 (MQ=255)
            gTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAg                                                                                                                                <  1:2590219/148‑1 (MQ=255)
             ttGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGc                                                                                                                               >  2:603034/1‑148 (MQ=255)
             ttGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGc                                                                                                                               <  1:1941532/148‑1 (MQ=255)
               gAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTc                                                                                                                             <  2:439547/148‑1 (MQ=255)
                aaaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTcc                                                                                                                            >  2:552030/1‑148 (MQ=255)
                 aaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCg                                                                                                                           >  2:930034/1‑148 (MQ=255)
                    gCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCca                                                                                                                        >  2:3238688/1‑148 (MQ=255)
                     ccAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCcac                                                                                                                       >  2:2738838/1‑148 (MQ=255)
                         aTTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCa                                                                                                                   <  2:875206/148‑1 (MQ=255)
                         aTTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCa                                                                                                                   >  2:3329239/1‑148 (MQ=255)
                          ttCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCaa                                                                                                                  >  1:721868/1‑148 (MQ=255)
                           tCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAg                                                                                                                 <  1:55530/148‑1 (MQ=255)
                             aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTAAATCTCTGGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGa                                                                                                               <  1:1219015/148‑1 (MQ=255)
                               aaaGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATg                                                                                                             >  1:177198/1‑148 (MQ=255)
                                aaGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATgg                                                                                                            >  1:1846930/1‑148 (MQ=255)
                                 aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGAGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                           <  2:2141244/148‑1 (MQ=255)
                                  gcgcGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCg                                                                                                          <  1:1925879/148‑1 (MQ=255)
                                    gcgTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTg                                                                                                        <  1:449487/148‑1 (MQ=255)
                                       tCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGGGGTc                                                                                                     >  2:540678/1‑148 (MQ=255)
                                        cATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTAAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                    >  2:302796/1‑148 (MQ=255)
                                         aTCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAg                                                                                                   >  2:2079716/1‑148 (MQ=255)
                                          tctcTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgt                                                                                                  >  1:2479650/1‑148 (MQ=255)
                                           ctctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgtg                                                                                                 >  1:1405001/1‑148 (MQ=255)
                                             ctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgtt                                                                                               >  1:2130816/1‑148 (MQ=255)
                                               ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttgt                                                                                             >  1:3236876/1‑148 (MQ=255)
                                               ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttgt                                                                                             >  1:1333068/1‑148 (MQ=255)
                                                 ctctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAt                                                                                           <  2:1892148/148‑1 (MQ=255)
                                                   ctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAttt                                                                                         >  2:933822/1‑148 (MQ=255)
                                                   ctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAttt                                                                                         <  1:2780705/148‑1 (MQ=255)
                                                      aCGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGt                                                                                      <  1:3006856/148‑1 (MQ=255)
                                                      aCGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGt                                                                                      >  2:2012289/1‑148 (MQ=255)
                                                         cccAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAg                                                                                   >  2:913751/1‑148 (MQ=255)
                                                             gcgcGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAata                                                                               <  2:1627565/148‑1 (MQ=255)
                                                             gcgcGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAata                                                                               <  1:1803628/148‑1 (MQ=255)
                                                             gcgcGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAata                                                                               <  1:1433313/148‑1 (MQ=255)
                                                                  gCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTc                                                                          <  1:1095087/148‑1 (MQ=255)
                                                                  gCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTc                                                                          <  1:3112442/148‑1 (MQ=255)
                                                                     cagaGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAAGACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAg                                                                       <  2:2245092/145‑1 (MQ=255)
                                                                       tAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGcc                                                                     >  2:502523/1‑148 (MQ=255)
                                                                       tAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGcc                                                                     >  2:2826639/1‑148 (MQ=255)
                                                                        aGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCa                                                                    <  1:747363/148‑1 (MQ=255)
                                                                         gTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCAt                                                                   >  1:2505062/1‑148 (MQ=255)
                                                                           cAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTc                                                                 <  2:1947112/148‑1 (MQ=255)
                                                                              ccccaccccaCACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTCTTTGGTCAGAATATCGTCCAGCCATTCCaa                                                              >  1:1387379/1‑148 (MQ=255)
                                                                              ccccaccccaCACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCaa                                                              <  2:359051/148‑1 (MQ=255)
                                                                              ccccaccccaCACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATACCaa                                                              >  1:1824896/1‑148 (MQ=255)
                                                                                     ccacACGGCAACCAGCCACTGTTCCGGTTCAATCTTTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACtt                                                       >  2:1543429/1‑148 (MQ=255)
                                                                                     ccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACtt                                                       >  2:1785492/1‑148 (MQ=255)
                                                                                     ccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACtt                                                       <  1:1284852/148‑1 (MQ=255)
                                                                                        cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCg                                                    >  2:126330/1‑148 (MQ=255)
                                                                                        cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCg                                                    >  2:1939989/1‑148 (MQ=255)
                                                                                         acGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGt                                                   >  1:140217/1‑148 (MQ=255)
                                                                                         acGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGt                                                   <  2:957677/148‑1 (MQ=255)
                                                                                                ccAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACgg                                            <  1:2220031/148‑1 (MQ=255)
                                                                                                ccAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACgg                                            >  1:2998510/1‑148 (MQ=255)
                                                                                                ccAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATGTGGTCAGAATATCGTCCAGCCATTCCAACTCACTGTCGGTTAACgg                                            >  1:2871676/1‑148 (MQ=255)
                                                                                                  aGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTc                                          <  1:930034/148‑1 (MQ=255)
                                                                                                  aGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTc                                          >  2:1443062/1‑148 (MQ=255)
                                                                                                  aGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTAACTTTCGTTTAACGGTc                                          <  1:1065809/148‑1 (MQ=255)
                                                                                                    ccACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTccc                                        >  2:1142491/1‑148 (MQ=255)
                                                                                                     cACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCg                                       <  2:3236876/148‑1 (MQ=255)
                                                                                                     cACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCg                                       <  2:255018/148‑1 (MQ=255)
                                                                                                      aCTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGt                                      >  1:1724461/1‑148 (MQ=255)
                                                                                                          ttCCGTTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTc                                  <  1:3238688/148‑1 (MQ=255)
                                                                                                           tgcGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCa                                 >  2:1019542/3‑148 (MQ=255)
                                                                                                            ccGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCAt                                >  1:2364415/1‑148 (MQ=255)
                                                                                                             cGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTGCGTTTAACGGTCCCGTTTTCATa                               >  2:1411316/1‑148 (MQ=255)
                                                                                                             cGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATAGCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATa                               >  2:3161160/1‑148 (MQ=255)
                                                                                                              ggTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATAc                              >  1:2799443/1‑148 (MQ=255)
                                                                                                                ttCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGc                            <  2:1857278/148‑1 (MQ=255)
                                                                                                                ttCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGc                            >  1:316315/1‑148 (MQ=255)
                                                                                                                       tctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCt                     >  2:2994669/1‑148 (MQ=255)
                                                                                                                       tctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCt                     <  1:406277/148‑1 (MQ=255)
                                                                                                                          tGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTtgt                  <  1:80951/148‑1 (MQ=255)
                                                                                                                            cGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTgg                <  1:1056362/148‑1 (MQ=255)
                                                                                                                              gagaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAt              <  1:2692618/148‑1 (MQ=255)
                                                                                                                              gagaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAt              >  2:2390284/1‑148 (MQ=255)
                                                                                                                                gaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtct            >  2:1438016/1‑148 (MQ=255)
                                                                                                                                gaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtct            >  1:3127332/1‑148 (MQ=255)
                                                                                                                                  aCTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATCTCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTGTTCATACGCTTTTCCTTGTGGATCTCa          >  2:762404/1‑148 (MQ=255)
                                                                                                                                  aCTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCa          <  2:2347712/148‑1 (MQ=255)
                                                                                                                                    tCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAc        <  1:680171/148‑1 (MQ=255)
                                                                                                                                    tCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAc        >  1:521819/1‑148 (MQ=255)
                                                                                                                                    tCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAc        <  2:2407903/148‑1 (MQ=255)
                                                                                                                                     cAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACt       >  1:2394063/1‑148 (MQ=255)
                                                                                                                                      aaCACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTc      >  2:1006365/1‑148 (MQ=255)
                                                                                                                                        cacCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGc    >  2:217950/1‑148 (MQ=255)
                                                                                                                                          ccgcagTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCa  >  2:1229796/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCA  >  CP009273/1984845‑1985130

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
TCTGCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCT  >  CP009273/1984838‑1985136
                                                                                                                                                     |                                                                                                                                                     
TCTGCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCA                                                                                                                                                       >  SRR19973065.187863/1‑150 (MQ=60)
  TGCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGT                                                                                                                                                     <  SRR19973065.2827046/150‑1 (MQ=60)
       TATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATA                                                                                                                                                <  SRR19973065.1957899/150‑1 (MQ=60)
       TATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATA                                                                                                                                                <  SRR19973065.1274806/150‑1 (MQ=60)
           GGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACC                                                                                                                                            >  SRR19973065.1870637/1‑150 (MQ=60)
            GCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCG                                                                                                                                           <  SRR19973065.357004/150‑1 (MQ=60)
              CATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTC                                                                                                                                         <  SRR19973065.3079028/150‑1 (MQ=60)
              CATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTC                                                                                                                                         <  SRR19973065.3078972/150‑1 (MQ=60)
              CATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTC                                                                                                                                         <  SRR19973065.1353903/150‑1 (MQ=60)
               ATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTTC                                                                                                                                        <  SRR19973065.613623/150‑1 (MQ=60)
               ATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCC                                                                                                                                        <  SRR19973065.244932/150‑1 (MQ=60)
               ATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCC                                                                                                                                        <  SRR19973065.2310821/150‑1 (MQ=60)
                TATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCA                                                                                                                                       >  SRR19973065.3163166/1‑150 (MQ=60)
                TATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCA                                                                                                                                       >  SRR19973065.3163140/1‑150 (MQ=60)
                   GTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCT                                                                                                                                    <  SRR19973065.2606337/150‑1 (MQ=60)
                    TTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTC                                                                                                                                   <  SRR19973065.1954230/150‑1 (MQ=60)
                    TTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTC                                                                                                                                   <  SRR19973065.1954227/150‑1 (MQ=60)
                      GAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCG                                                                                                                                 >  SRR19973065.2192892/1‑150 (MQ=60)
                               GATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAA                                                                                                                        >  SRR19973065.726698/1‑150 (MQ=60)
                                  TCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGA                                                                                                                     <  SRR19973065.54922/150‑1 (MQ=60)
                                  TCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGA                                                                                                                     <  SRR19973065.56128/150‑1 (MQ=60)
                                   CATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGAT                                                                                                                    >  SRR19973065.1956701/1‑150 (MQ=60)
                                    ATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATG                                                                                                                   >  SRR19973065.178600/1‑150 (MQ=60)
                                    ATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTAAATCTCTGGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATG                                                                                                                   <  SRR19973065.1226932/150‑1 (MQ=60)
                                     TAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGG                                                                                                                  >  SRR19973065.1859519/1‑150 (MQ=60)
                                     TAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGG                                                                                                                  >  SRR19973065.1859099/1‑150 (MQ=60)
                                       AAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCG                                                                                                                <  SRR19973065.2603575/150‑1 (MQ=60)
                                         GCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTG                                                                                                              <  SRR19973065.1938501/150‑1 (MQ=60)
                                           GCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGT                                                                                                            <  SRR19973065.2713065/150‑1 (MQ=60)
                                           GCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTGGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGT                                                                                                            <  SRR19973065.452569/150‑1 (MQ=60)
                                           GCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGT                                                                                                            <  SRR19973065.452605/150‑1 (MQ=60)
                                               CATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGT                                                                                                        >  SRR19973065.2495262/1‑150 (MQ=60)
                                                ATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTG                                                                                                       >  SRR19973065.1414216/1‑150 (MQ=60)
                                                ATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTG                                                                                                       >  SRR19973065.1414208/1‑150 (MQ=60)
                                                ATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTG                                                                                                       >  SRR19973065.91735/1‑150 (MQ=60)
                                                  CTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTT                                                                                                     >  SRR19973065.2144528/1‑150 (MQ=60)
                                                   TCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTG                                                                                                    <  SRR19973065.3182156/150‑1 (MQ=60)
                                                   TCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTG                                                                                                    <  SRR19973065.753181/150‑1 (MQ=60)
                                                    CTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGT                                                                                                   >  SRR19973065.3256363/1‑150 (MQ=60)
                                                    CTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGT                                                                                                   >  SRR19973065.3256451/1‑150 (MQ=60)
                                                    CTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGT                                                                                                   >  SRR19973065.1341892/1‑150 (MQ=60)
                                                          CTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACCCGGCAACCAGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGG                                                                                             <  SRR19973065.2797958/150‑1 (MQ=60)
                                                          CTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGG                                                                                             <  SRR19973065.2798341/150‑1 (MQ=60)
                                                           TGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGT                                                                                            >  SRR19973065.178906/1‑150 (MQ=60)
                                                             ACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCA                                                                                          <  SRR19973065.3025274/150‑1 (MQ=60)
                                                                   AGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATAT                                                                                    <  SRR19973065.1990811/150‑1 (MQ=60)
                                                                    GCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATC                                                                                   <  SRR19973065.1815728/150‑1 (MQ=60)
                                                                    GCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATC                                                                                   <  SRR19973065.1815482/150‑1 (MQ=60)
                                                                    GCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATC                                                                                   <  SRR19973065.1442624/150‑1 (MQ=60)
                                                                         GCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCA                                                                              <  SRR19973065.1102295/150‑1 (MQ=60)
                                                                         GCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCA                                                                              <  SRR19973065.3131411/150‑1 (MQ=60)
                                                                         GCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTACGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCA                                                                              <  SRR19973065.2073165/150‑1 (MQ=60)
                                                                          CACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAG                                                                             >  SRR19973065.361575/1‑150 (MQ=60)
                                                                             ATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCA                                                                          >  SRR19973065.749951/1‑150 (MQ=60)
                                                                              TAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCAT                                                                         >  SRR19973065.2520807/1‑150 (MQ=60)
                                                                              TAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCAT                                                                         >  SRR19973065.2520801/1‑150 (MQ=60)
                                                                               AGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATT                                                                        >  SRR19973065.3270027/1‑150 (MQ=60)
                                                                               AGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATT                                                                        <  SRR19973065.752361/150‑1 (MQ=60)
                                                                                   AGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTCTTTGGTCAGAATATCGTCCAGCCATTCCAA                                                                    >  SRR19973065.1396494/1‑150 (MQ=60)
                                                                                   AGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAA                                                                    >  SRR19973065.1757594/1‑150 (MQ=60)
                                                                                   AGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAA                                                                    >  SRR19973065.1396325/1‑150 (MQ=60)
                                                                                   AGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATACCAA                                                                    >  SRR19973065.1836912/1‑150 (MQ=60)
                                                                                    GCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAAC                                                                   <  SRR19973065.495050/150‑1 (MQ=60)
                                                                                     CCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACT                                                                  <  SRR19973065.3136495/150‑1 (MQ=60)
                                                                                         ACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACT                                                              <  SRR19973065.604291/150‑1 (MQ=60)
                                                                                           CCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTT                                                            <  SRR19973065.1397525/150‑1 (MQ=60)
                                                                                           CCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTT                                                            <  SRR19973065.1326148/150‑1 (MQ=60)
                                                                                           CCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTT                                                            >  SRR19973065.295569/1‑150 (MQ=60)
                                                                                            CCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTC                                                           <  SRR19973065.1293327/150‑1 (MQ=60)
                                                                                             CACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCG                                                          >  SRR19973065.1085198/1‑150 (MQ=60)
                                                                                              ACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGT                                                         <  SRR19973065.3424268/150‑1 (MQ=60)
                                                                                              ACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGT                                                         >  SRR19973065.141419/1‑150 (MQ=60)
                                                                                              ACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGT                                                         >  SRR19973065.141542/1‑150 (MQ=60)
                                                                                                  GGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAA                                                     >  SRR19973065.725572/1‑150 (MQ=60)
                                                                                                     AACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGG                                                  >  SRR19973065.3016877/1‑150 (MQ=60)
                                                                                                     AACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGG                                                  <  SRR19973065.3099626/150‑1 (MQ=60)
                                                                                                     AACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGG                                                  >  SRR19973065.3016885/1‑150 (MQ=60)
                                                                                                     AACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATGTGGTCAGAATATCGTCCAGCCATTCCAACTCACTGTCGGTTAACGG                                                  >  SRR19973065.2889276/1‑150 (MQ=60)
                                                                                                       CCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTC                                                <  SRR19973065.2234200/150‑1 (MQ=60)
                                                                                                         AGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCC                                              <  SRR19973065.936161/150‑1 (MQ=60)
                                                                                                         AGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTAACTTTCGTTTAACGGTCCC                                              <  SRR19973065.1072869/150‑1 (MQ=60)
                                                                                                         AGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCC                                              <  SRR19973065.936154/150‑1 (MQ=60)
                                                                                                           CCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGT                                            >  SRR19973065.1735903/1‑150 (MQ=60)
                                                                                                                GTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTGTCA                                       >  SRR19973065.1353145/1‑150 (MQ=60)
                                                                                                                 TTCCGTTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCAT                                      <  SRR19973065.3258191/150‑1 (MQ=60)
                                                                                                                 TTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCAT                                      >  SRR19973065.2379135/1‑150 (MQ=60)
                                                                                                                   CCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATAC                                    >  SRR19973065.2816787/1‑150 (MQ=60)
                                                                                                                     GGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGC                                  >  SRR19973065.318472/1‑150 (MQ=60)
                                                                                                                     GGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGC                                  >  SRR19973065.1842305/1‑150 (MQ=60)
                                                                                                                              TCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTG                         <  SRR19973065.409122/150‑1 (MQ=60)
                                                                                                                               CTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGT                        >  SRR19973065.3253498/1‑150 (MQ=60)
                                                                                                                                 TGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGG                      <  SRR19973065.81758/150‑1 (MQ=60)
                                                                                                                                 TGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGG                      <  SRR19973065.1870679/150‑1 (MQ=60)
                                                                                                                                   CGGAGAACTCAACACCGCAGTCAATAACCCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCATGTGGAT                    <  SRR19973065.1063430/150‑1 (MQ=60)
                                                                                                                                   CGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAT                    <  SRR19973065.2863197/150‑1 (MQ=60)
                                                                                                                                   CGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAT                    <  SRR19973065.1063367/150‑1 (MQ=60)
                                                                                                                                     GAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCT                  >  SRR19973065.3146371/1‑150 (MQ=60)
                                                                                                                                     GAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCT                  <  SRR19973065.2709336/150‑1 (MQ=60)
                                                                                                                                      AGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCCAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTC                 <  SRR19973065.1521711/150‑1 (MQ=60)
                                                                                                                                      AGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATAGCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATATC                 >  SRR19973065.2643627/1‑150 (MQ=60)
                                                                                                                                         ACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAC              >  SRR19973065.1156665/1‑150 (MQ=60)
                                                                                                                                         ACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAC              >  SRR19973065.525254/1‑150 (MQ=60)
                                                                                                                                          CTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACT             >  SRR19973065.2408993/1‑150 (MQ=60)
                                                                                                                                          CTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACT             >  SRR19973065.2409111/1‑150 (MQ=60)
                                                                                                                                           TCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTC            <  SRR19973065.684477/150‑1 (MQ=60)
                                                                                                                                           TCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTC            <  SRR19973065.684700/150‑1 (MQ=60)
                                                                                                                                                     AGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCT  <  SRR19973065.2885841/150‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
TCTGCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCT  >  CP009273/1984838‑1985136

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