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 F361 I0 R1 1353 250.7 8230166 96.0% 7900959 148.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,820,888 C→A 100% C98F (TGC→TTC)  astE ← succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,820,8880CA99.0% 643.1 / ‑2.0 195C98F (TGC→TTC) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (97/96);  total (99/96)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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 198 total aligned reads displayed.

GTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGC  >  CP009273/1820742‑1821033
                                                                                                                                                  |                                                                                                                                                 
gTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGaaa                                                                                                                                                 <  2:1481872/149‑1 (MQ=255)
gTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGaaa                                                                                                                                                 >  1:1817670/1‑149 (MQ=255)
   gATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGAAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATaa                                                                                                                                              <  1:2887786/149‑1 (MQ=255)
    aTCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAc                                                                                                                                             <  1:3017556/149‑1 (MQ=255)
    aTCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAc                                                                                                                                             <  2:159343/149‑1 (MQ=255)
      cAAGGTGCCAGCGAACAGATTCTATGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGAAAGGTTTCTCCGCTTACAGCAAAAAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                           >  2:3369275/1‑149 (MQ=255)
       aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTCCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGtt                                                                                                                                          <  2:2627381/149‑1 (MQ=255)
         ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTc                                                                                                                                        <  2:562134/149‑1 (MQ=255)
         ggTGCCAGCGCACAGATTCTTTGCCCTGGTAATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTc                                                                                                                                        <  1:2545433/149‑1 (MQ=255)
             ccAGCGCACATATCCTTTGCCCTGGTCATAAAAATCATCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCtt                                                                                                                                    <  2:2848714/149‑1 (MQ=255)
             ccAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCtt                                                                                                                                    <  1:1730880/149‑1 (MQ=255)
                gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCt                                                                                                                                 <  1:1018295/149‑1 (MQ=255)
                  gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTc                                                                                                                               <  2:1821215/149‑1 (MQ=255)
                       gATTCTTTGCCCTGGTCATAAAAATCTTCCCGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAgcgc                                                                                                                          <  1:1206111/149‑1 (MQ=255)
                       gATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAgcgc                                                                                                                          >  2:1948295/1‑149 (MQ=255)
                       gATTCTTTGCCCTGGTCAAAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAgcgc                                                                                                                          <  1:1777992/149‑1 (MQ=255)
                         ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCag                                                                                                                        <  2:1934276/149‑1 (MQ=255)
                         ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCACCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCag                                                                                                                        >  1:1673268/1‑149 (MQ=255)
                            tttGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCaggag                                                                                                                     >  1:1624997/1‑149 (MQ=255)
                             ttGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCaggagg                                                                                                                    <  1:558293/149‑1 (MQ=255)
                                cccTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGAtt                                                                                                                 <  2:3996885/149‑1 (MQ=255)
                                cccTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGAtt                                                                                                                 <  2:669830/149‑1 (MQ=255)
                                    ggTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCg                                                                                                             <  1:1558893/149‑1 (MQ=255)
                                     gTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCACCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGa                                                                                                            <  2:1673268/149‑1 (MQ=255)
                                        aTAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGa                                                                                                         >  1:2616991/1‑149 (MQ=255)
                                         tAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGAt                                                                                                        >  2:4107137/1‑149 (MQ=255)
                                         tAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGAt                                                                                                        <  2:3711622/149‑1 (MQ=255)
                                         tAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGAt                                                                                                        >  2:219867/1‑149 (MQ=255)
                                          aaaaaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCGGCGCAGGAGGATTCCCGAGGATc                                                                                                       <  1:1163047/149‑1 (MQ=255)
                                          aaaaaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATc                                                                                                       <  1:3854136/149‑1 (MQ=255)
                                           aaaaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCa                                                                                                      >  2:1269556/1‑149 (MQ=255)
                                           aaaaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCa                                                                                                      >  2:2740659/1‑149 (MQ=255)
                                            aaaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAc                                                                                                     <  1:2215693/149‑1 (MQ=255)
                                             aaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAcc                                                                                                    <  1:979465/149‑1 (MQ=255)
                                              aTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTATCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCa                                                                                                   <  2:3797816/149‑1 (MQ=255)
                                                 ttCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCa                                                                                                >  2:815697/1‑149 (MQ=255)
                                                   ccAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAc                                                                                              <  1:3308913/149‑1 (MQ=255)
                                                     aGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCg                                                                                            >  2:2285230/1‑149 (MQ=255)
                                                         aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCCGCAACCGCCaa                                                                                        >  2:1050236/1‑149 (MQ=255)
                                                         aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCaa                                                                                        >  2:1052249/1‑149 (MQ=255)
                                                          cTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAAc                                                                                       >  2:273806/1‑149 (MQ=255)
                                                          cTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAAc                                                                                       >  2:1002292/1‑149 (MQ=255)
                                                              tCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTaa                                                                                   <  1:274257/149‑1 (MQ=255)
                                                              tCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTaa                                                                                   <  1:2273521/149‑1 (MQ=255)
                                                                cAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACg                                                                                 >  1:848396/1‑149 (MQ=255)
                                                                 aGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACgg                                                                                >  2:1434504/1‑149 (MQ=255)
                                                                  gTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCACCGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                                               <  2:3234886/149‑1 (MQ=255)
                                                                  gTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                                               <  2:3234768/149‑1 (MQ=255)
                                                                    tCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAt                                                                             >  2:4065653/1‑149 (MQ=255)
                                                                      gcgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtct                                                                           >  1:4113485/1‑149 (MQ=255)
                                                                       cgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtctc                                                                          <  2:814311/149‑1 (MQ=255)
                                                                               aCAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTgag                                                                  <  1:691406/149‑1 (MQ=255)
                                                                                cAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTgaga                                                                 <  2:1653127/149‑1 (MQ=255)
                                                                                 aGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGat                                                                >  1:3643935/1‑149 (MQ=255)
                                                                                  ggTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGata                                                               <  1:814284/149‑1 (MQ=255)
                                                                                  ggTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGAAACGGGATCTCGCCGTGAGata                                                               >  2:1992497/1‑149 (MQ=255)
                                                                                    tttCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATc                                                             >  2:2526727/1‑149 (MQ=255)
                                                                                     ttCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATct                                                            <  2:3905994/149‑2 (MQ=255)
                                                                                      tctcCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATcgc                                                           <  2:53306/149‑1 (MQ=255)
                                                                                      tctcCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATcgc                                                           <  2:2317966/149‑1 (MQ=255)
                                                                                       ctcCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATcgcg                                                          >  2:3524269/1‑149 (MQ=255)
                                                                                             tttCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGc                                                    <  1:3196201/149‑1 (MQ=255)
                                                                                                    aaaTAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGt                                             >  2:1745389/1‑149 (MQ=255)
                                                                                                      aTAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTcc                                           <  2:3952851/149‑1 (MQ=255)
                                                                                                       tAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCa                                          <  2:2812793/149‑1 (MQ=255)
                                                                                                         gCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCGACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGc                                        <  2:888874/149‑1 (MQ=255)
                                                                                                              cAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAtctc                                   >  1:3333111/1‑149 (MQ=255)
                                                                                                                aCGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCca                                 <  2:699463/149‑1 (MQ=255)
                                                                                                                 cGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCcac                                >  2:2392253/1‑149 (MQ=255)
                                                                                                                    ccGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAgg                             >  1:3763563/1‑149 (MQ=255)
                                                                                                                     cGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGt                            <  1:3685478/149‑1 (MQ=255)
                                                                                                                      gCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTg                           >  2:2124143/1‑149 (MQ=255)
                                                                                                                      gCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTg                           >  1:4087344/1‑149 (MQ=255)
                                                                                                                       ccAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGc                          <  1:3183280/149‑1 (MQ=255)
                                                                                                                        cAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGcc                         >  2:1399420/1‑149 (MQ=255)
                                                                                                                        cAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGcc                         >  1:1631202/1‑149 (MQ=255)
                                                                                                                        cAAACATTCGAATCATATCGCTATGGAAATAACGTATCCCATGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTCAGAGATCCCGCCAAGCAACGCCTCCAGCATCTACACAGGTGcc                         >  2:3447495/1‑149 (MQ=255)
                                                                                                                           aCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGtc                      <  1:936491/149‑1 (MQ=255)
                                                                                                                            cATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGtct                     >  1:1609429/1‑149 (MQ=255)
                                                                                                                            cATTCGATTCAGATCGCTATGGAAATAACGTGTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCCTCTCCACAGGTGACGtct                     <  2:2142192/149‑1 (MQ=255)
                                                                                                                             aTTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGtctc                    >  2:3092973/1‑149 (MQ=255)
                                                                                                                             aTTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGtctc                    >  2:3164464/1‑149 (MQ=255)
                                                                                                                              ttCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCa                   >  2:3009057/1‑149 (MQ=255)
                                                                                                                                cGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCAtt                 >  1:4048174/1‑149 (MQ=255)
                                                                                                                                 gATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTa                >  1:2618818/1‑149 (MQ=255)
                                                                                                                                  aTTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc               >  1:1067984/1‑149 (MQ=255)
                                                                                                                                  aTTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc               >  1:1355228/1‑149 (MQ=255)
                                                                                                                                  aTTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc               >  1:3120084/1‑149 (MQ=255)
                                                                                                                                  aTTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc               <  2:1609264/149‑1 (MQ=255)
                                                                                                                                   ttCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAcc              >  2:1142539/1‑149 (MQ=255)
                                                                                                                                     cATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgt            >  1:2207542/1‑149 (MQ=255)
                                                                                                                                      atatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtg           >  2:3987498/1‑149 (MQ=255)
                                                                                                                                       tatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCTTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtgt          >  1:2988825/1‑149 (MQ=255)
                                                                                                                                       tatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtgt          >  2:632861/1‑149 (MQ=255)
                                                                                                                                        atCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTa         >  1:630075/1‑149 (MQ=255)
                                                                                                                                         tCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAt        >  1:1591791/1‑149 (MQ=255)
                                                                                                                                          cGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAtt       <  1:1948295/149‑1 (MQ=255)
                                                                                                                                           gCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTc      <  2:770697/149‑1 (MQ=255)
                                                                                                                                            cTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTcc     <  1:3973653/149‑1 (MQ=255)
                                                                                                                                              aTGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCg   <  2:1695960/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGC  >  CP009273/1820742‑1821033

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
TATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGA  >  CP009273/1820739‑1821036
                                                                                                                                                     |                                                                                                                                                    
TATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGA                                                                                                                                                      >  SRR19973068.775913/1‑150 (MQ=60)
 ATGTAGATCACGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAA                                                                                                                                                     <  SRR19973068.2831805/150‑1 (MQ=60)
  TGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAA                                                                                                                                                    >  SRR19973068.1829614/1‑150 (MQ=60)
      GATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAC                                                                                                                                                <  SRR19973068.2905799/150‑1 (MQ=60)
      GATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAC                                                                                                                                                <  SRR19973068.994198/150‑1 (MQ=60)
      GATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGAAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAC                                                                                                                                                <  SRR19973068.2905761/150‑1 (MQ=60)
       ATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACG                                                                                                                                               <  SRR19973068.3036141/150‑1 (MQ=60)
       ATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACG                                                                                                                                               <  SRR19973068.1889966/150‑1 (MQ=60)
            GGTGCCAGCGCACAGATTCTTTGCCCTGGTAATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCC                                                                                                                                          <  SRR19973068.2561591/150‑1 (MQ=60)
                CCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTG                                                                                                                                      <  SRR19973068.1741854/150‑1 (MQ=60)
                 CAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGC                                                                                                                                     >  SRR19973068.3210186/1‑150 (MQ=60)
                   GCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTT                                                                                                                                   <  SRR19973068.1024906/150‑1 (MQ=60)
                        CAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCG                                                                                                                              <  SRR19973068.2292864/150‑1 (MQ=60)
                          GATTCTTTGCCCTGGTCATAAAAATCTTCCCGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCA                                                                                                                            <  SRR19973068.1214158/150‑1 (MQ=60)
                          GATTCTTTGCCCTGGTCAAAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCA                                                                                                                            <  SRR19973068.1789547/150‑1 (MQ=60)
                           ATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAG                                                                                                                           >  SRR19973068.3387498/1‑150 (MQ=60)
                           ATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAG                                                                                                                           >  SRR19973068.1619805/1‑150 (MQ=60)
                           ATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCACCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAG                                                                                                                           >  SRR19973068.1684060/1‑150 (MQ=60)
                           ATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCACCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAG                                                                                                                           >  SRR19973068.1684070/1‑150 (MQ=60)
                              CTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAG                                                                                                                        >  SRR19973068.1635602/1‑150 (MQ=60)
                              CTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAG                                                                                                                        >  SRR19973068.1635600/1‑150 (MQ=60)
                                TTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGA                                                                                                                      <  SRR19973068.562138/150‑1 (MQ=60)
                                   CCCTGGTCATAAAAATCTTCCAGGCACTGTTCATGTTCGCGCGCCCGACAGGTTTCTCCGCTTACAGCAAATAGCTGCCCACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTC                                                                                                                   <  SRR19973068.3153989/150‑1 (MQ=60)
                                       GGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGA                                                                                                               <  SRR19973068.1569154/150‑1 (MQ=60)
                                        GTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAG                                                                                                              >  SRR19973068.2830288/1‑150 (MQ=60)
                                        GTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAG                                                                                                              >  SRR19973068.1763635/1‑150 (MQ=60)
                                          CATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGA                                                                                                            >  SRR19973068.2633472/1‑150 (MQ=60)
                                          CATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGA                                                                                                            >  SRR19973068.2632729/1‑150 (MQ=60)
                                          CATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGA                                                                                                            >  SRR19973068.325969/1‑150 (MQ=60)
                                          CATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGA                                                                                                            >  SRR19973068.1423673/1‑150 (MQ=60)
                                             AAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCGGCGCAGGAGGATTCCCGAGGATCA                                                                                                         <  SRR19973068.1170821/150‑1 (MQ=60)
                                             AAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCA                                                                                                         <  SRR19973068.1170846/150‑1 (MQ=60)
                                             AAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCA                                                                                                         <  SRR19973068.3877407/150‑1 (MQ=60)
                                               AAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACC                                                                                                       <  SRR19973068.2230179/150‑1 (MQ=60)
                                                AATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCA                                                                                                      <  SRR19973068.985934/150‑1 (MQ=60)
                                                   CTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCA                                                                                                   <  SRR19973068.1829418/150‑1 (MQ=60)
                                                      CCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACC                                                                                                <  SRR19973068.3329147/150‑1 (MQ=60)
                                                                 TCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAAC                                                                                     <  SRR19973068.2288263/150‑1 (MQ=60)
                                                                 TCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAAC                                                                                     <  SRR19973068.2288014/150‑1 (MQ=60)
                                                                 TCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAAC                                                                                     <  SRR19973068.3713685/150‑1 (MQ=60)
                                                                 TCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAAC                                                                                     <  SRR19973068.276162/150‑1 (MQ=60)
                                                                  CCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACG                                                                                    >  SRR19973068.854172/1‑150 (MQ=60)
                                                                        CGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCT                                                                              >  SRR19973068.4138284/1‑150 (MQ=60)
                                                                        CGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCT                                                                              >  SRR19973068.4138088/1‑150 (MQ=60)
                                                                                  ACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGA                                                                    <  SRR19973068.695757/150‑1 (MQ=60)
                                                                                  ACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGA                                                                    >  SRR19973068.1638074/1‑150 (MQ=60)
                                                                                   CAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGAT                                                                   <  SRR19973068.785328/150‑1 (MQ=60)
                                                                                   CAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGAT                                                                   >  SRR19973068.3666008/1‑150 (MQ=60)
                                                                                     GGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATAT                                                                 <  SRR19973068.819829/150‑1 (MQ=60)
                                                                                     GGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATAT                                                                 <  SRR19973068.829215/150‑1 (MQ=60)
                                                                                            CCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCA                                                          >  SRR19973068.611788/1‑150 (MQ=60)
                                                                                                TTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCA                                                      <  SRR19973068.3215867/150‑1 (MQ=60)
                                                                                                     GCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCG                                                 <  SRR19973068.4131713/150‑1 (MQ=60)
                                                                                                     GCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCG                                                 <  SRR19973068.3018719/150‑1 (MQ=60)
                                                                                                         ATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCA                                             <  SRR19973068.972722/150‑1 (MQ=60)
                                                                                                               GCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCT                                       <  SRR19973068.1120984/150‑1 (MQ=60)
                                                                                                               GCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCT                                       <  SRR19973068.3910725/150‑1 (MQ=60)
                                                                                                                CCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTC                                      >  SRR19973068.3353449/1‑150 (MQ=60)
                                                                                                                CCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTC                                      >  SRR19973068.3353197/1‑150 (MQ=60)
                                                                                                                  ATCGACCGCCAAACATTCGATTCATATCTCTATGGAAATAACGTTTACCTTGCTTCAGCGCAGGAGGATTCCCGTGGATCACCTGCTACCGCCATCGTAACGGGATCTCGCCGTGTGATATCGCGCCAAGGAACGCGTCCAGCATCTCCA                                    <  SRR19973068.1356905/150‑1 (MQ=60)
                                                                                                                      ACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGG                                >  SRR19973068.3786042/1‑150 (MQ=60)
                                                                                                                        CGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTG                              <  SRR19973068.3707707/150‑1 (MQ=60)
                                                                                                                        CGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTG                              >  SRR19973068.4110931/1‑150 (MQ=60)
                                                                                                                        CGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTG                              >  SRR19973068.4111853/1‑150 (MQ=60)
                                                                                                                          CCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCC                            >  SRR19973068.1641837/1‑150 (MQ=60)
                                                                                                                          CCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCC                            >  SRR19973068.1641914/1‑150 (MQ=60)
                                                                                                                          CCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCC                            <  SRR19973068.3202876/150‑1 (MQ=60)
                                                                                                                              ACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCT                        >  SRR19973068.3676112/1‑150 (MQ=60)
                                                                                                                              ACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCT                        >  SRR19973068.1619970/1‑150 (MQ=60)
                                                                                                                              ACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCT                        <  SRR19973068.942779/150‑1 (MQ=60)
                                                                                                                              ACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCT                        <  SRR19973068.942570/150‑1 (MQ=60)
                                                                                                                              ACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGAGTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCT                        <  SRR19973068.3816862/150‑1 (MQ=60)
                                                                                                                                  TCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATT                    >  SRR19973068.4072544/1‑150 (MQ=60)
                                                                                                                                   CGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTA                   >  SRR19973068.2635305/1‑150 (MQ=60)
                                                                                                                                   CGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTA                   >  SRR19973068.2635248/1‑150 (MQ=60)
                                                                                                                                    GATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAC                  >  SRR19973068.1074996/1‑150 (MQ=60)
                                                                                                                                    GATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAC                  >  SRR19973068.3139174/1‑150 (MQ=60)
                                                                                                                                    GATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAC                  >  SRR19973068.1363932/1‑150 (MQ=60)
                                                                                                                                     ATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACC                 <  SRR19973068.3108261/150‑1 (MQ=60)
                                                                                                                                       TCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGT               >  SRR19973068.2221969/1‑150 (MQ=60)
                                                                                                                                       TCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGT               >  SRR19973068.2221994/1‑150 (MQ=60)
                                                                                                                                        CATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTG              >  SRR19973068.1257904/1‑150 (MQ=60)
                                                                                                                                         ATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCTTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGT             >  SRR19973068.3007309/1‑150 (MQ=60)
                                                                                                                                         ATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGT             <  SRR19973068.2698925/150‑1 (MQ=60)
                                                                                                                                          TATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTA            >  SRR19973068.634364/1‑150 (MQ=60)
                                                                                                                                          TATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTA            >  SRR19973068.634233/1‑150 (MQ=60)
                                                                                                                                           ATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAT           >  SRR19973068.1602234/1‑150 (MQ=60)
                                                                                                                                           ATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAT           >  SRR19973068.1602200/1‑150 (MQ=60)
                                                                                                                                             CGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTC         <  SRR19973068.1961234/150‑1 (MQ=60)
                                                                                                                                               CTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCC       <  SRR19973068.3997685/150‑1 (MQ=60)
                                                                                                                                                  TGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCT    <  SRR19973068.4054628/150‑1 (MQ=60)
                                                                                                                                                  TGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCT    <  SRR19973068.2005641/150‑1 (MQ=60)
                                                                                                                                                    GAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGA  >  SRR19973068.637197/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                    
TATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGA  >  CP009273/1820739‑1821036

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: