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,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,8880CA98.9% 602.0 / ‑1.9 185C98F (TGC→TTC) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (96/87);  total (97/88)
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.96e-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 185 total aligned reads displayed.

AGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATT  >  CP009273/1820744‑1821028
                                                                                                                                                |                                                                                                                                            
aGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAAt                                                                                                                                           <  1:2084289/148‑1 (MQ=255)
 gATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATa                                                                                                                                          >  2:903251/1‑148 (MQ=255)
 gATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATa                                                                                                                                          <  1:2835605/148‑1 (MQ=255)
  aTCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTAATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATaa                                                                                                                                         <  2:2704662/148‑1 (MQ=255)
   tCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAc                                                                                                                                        >  2:1573938/1‑148 (MQ=255)
   tCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAc                                                                                                                                        >  2:1166176/1‑148 (MQ=255)
    cAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGTTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACg                                                                                                                                       >  1:1883753/1‑148 (MQ=255)
     aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                      <  2:104465/148‑1 (MQ=255)
     aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                      <  1:1099922/148‑1 (MQ=255)
       ggTGCCAGCGCACATATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGttt                                                                                                                                    <  1:1415769/148‑1 (MQ=255)
        gTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTc                                                                                                                                   >  1:3017471/1‑148 (MQ=255)
        gTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTc                                                                                                                                   >  1:2022574/1‑148 (MQ=255)
         tGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCCGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTcc                                                                                                                                  <  1:154180/148‑1 (MQ=255)
         tGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACCGGTTTCTCGGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTcc                                                                                                                                  <  1:645614/148‑1 (MQ=255)
           ccAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCt                                                                                                                                <  2:2983063/148‑1 (MQ=255)
           ccAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCt                                                                                                                                <  2:2982965/148‑1 (MQ=255)
              gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGc                                                                                                                             <  2:3205011/148‑1 (MQ=255)
              gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGc                                                                                                                             <  2:1461313/148‑1 (MQ=255)
              gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGc                                                                                                                             <  1:1096954/148‑1 (MQ=255)
              gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTCGc                                                                                                                             >  2:1232064/1‑148 (MQ=255)
              gcgcACAGATTCTGTGCCCTGGTCATAAAAATCTTCCAGGCACTGGTCCAGTTCGCGCGCCCGCCAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAACCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGc                                                                                                                             <  2:434623/148‑1 (MQ=255)
                gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                           <  1:2252483/148‑1 (MQ=255)
                gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                           <  1:3062068/148‑1 (MQ=255)
                gcAAAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                           <  1:2252397/148‑1 (MQ=255)
                 cacaGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTc                                                                                                                          <  2:760403/148‑1 (MQ=255)
                 cacaGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTc                                                                                                                          >  1:2134018/1‑148 (MQ=255)
                   caGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAg                                                                                                                        >  1:670071/1‑148 (MQ=255)
                   caGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAg                                                                                                                        >  1:547025/1‑148 (MQ=255)
                       ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACCGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                    <  1:1239352/148‑1 (MQ=255)
                       ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTTCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                    >  1:3249654/1‑148 (MQ=255)
                         cTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                  <  2:1983417/148‑1 (MQ=255)
                         cTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                  <  2:262074/148‑1 (MQ=255)
                         cTTTGCCCTGGTCATAAAAATCGTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                  <  2:2066888/148‑1 (MQ=255)
                          tttGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagga                                                                                                                 <  1:527633/148‑1 (MQ=255)
                           tggCCGTGTTCATAAAAATATTCCTGGCACTGTTCCAGTTCGCGCGCCCGACATTTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCaggag                                                                                                                <  2:621560/146‑1 (MQ=255)
                               ccTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGAtt                                                                                                            <  2:3407915/148‑1 (MQ=255)
                                     cATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAg                                                                                                      >  1:2419967/1‑148 (MQ=255)
                                     cATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAg                                                                                                      >  2:1876454/1‑148 (MQ=255)
                                     cATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAg                                                                                                      >  1:658142/1‑148 (MQ=255)
                                       tAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGa                                                                                                    <  2:1346449/148‑1 (MQ=255)
                                       tAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGa                                                                                                    >  2:2803205/1‑148 (MQ=255)
                                        aaaaaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGAt                                                                                                   <  1:889094/148‑1 (MQ=255)
                                        aaaaaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGAt                                                                                                   >  2:1840726/1‑148 (MQ=255)
                                         aaaaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATc                                                                                                  >  2:2412401/1‑148 (MQ=255)
                                           aaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAc                                                                                                <  2:1724074/148‑1 (MQ=255)
                                                tCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAAAATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCa                                                                                           <  1:936607/148‑1 (MQ=255)
                                                     gCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGc                                                                                      >  1:3397899/1‑148 (MQ=255)
                                                       aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCa                                                                                    >  1:88018/1‑148 (MQ=255)
                                                        cTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCaa                                                                                   >  1:13777/1‑148 (MQ=255)
                                                        cTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCaa                                                                                   >  1:13878/1‑148 (MQ=255)
                                                         tGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAAc                                                                                  >  1:2596240/1‑148 (MQ=255)
                                                         tGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAAc                                                                                  <  1:2874343/148‑1 (MQ=255)
                                                         tGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAAc                                                                                  <  1:1123678/148‑1 (MQ=255)
                                                            tCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTa                                                                               <  1:763250/148‑1 (MQ=255)
                                                             ccAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTaa                                                                              <  2:2036461/148‑1 (MQ=255)
                                                               aGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACg                                                                            >  2:2526437/1‑148 (MQ=255)
                                                                 ttCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                                          <  2:3492336/148‑1 (MQ=255)
                                                                  tCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGa                                                                         <  2:1650033/148‑1 (MQ=255)
                                                                  tCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGa                                                                         >  2:2946218/1‑148 (MQ=255)
                                                                              cAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTgag                                                             <  1:1831891/148‑1 (MQ=255)
                                                                               aGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTgaga                                                            <  2:504935/148‑1 (MQ=255)
                                                                                tgTTTCTCCGCTTTCAGAAATTAGCTGCCAACGACCGCCAAACATTCGATTCTTATCGCTTTGGAAATACCGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGat                                                           <  2:539008/147‑1 (MQ=255)
                                                                                ggTTTCTCCGCTTTCAGCAAAGAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGat                                                           <  2:1082810/148‑1 (MQ=255)
                                                                                  tttCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGatat                                                         <  1:2729876/148‑1 (MQ=255)
                                                                                       ccGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCgcgc                                                    >  2:2981627/1‑148 (MQ=255)
                                                                                         gCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCa                                                  <  2:1920956/148‑1 (MQ=255)
                                                                                         gCTTTCAGCAAAGAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCa                                                  <  2:1920959/148‑1 (MQ=255)
                                                                                            ttCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGc                                               >  2:1959688/1‑148 (MQ=255)
                                                                                            ttCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGc                                               >  1:1546407/1‑148 (MQ=255)
                                                                                                    aTAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTc                                       >  1:1573250/1‑148 (MQ=255)
                                                                                                    aTAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTc                                       >  1:429682/1‑148 (MQ=255)
                                                                                                     tAGCAGCCAACGACCGCCAAACAGTCGATTCATATAGCAATGGAAATAACGTTTCCCTAGCTTCAGCTCAGGAGGATTCACGAGGATCACCAGCAACCGCAAAAGTAACGGGCTCTCGCCGTGAGATCTCGCGCCAAGCAAAGCGTcc                                      >  2:1603532/1‑148 (MQ=255)
                                                                                                       gCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAg                                    <  2:2531659/148‑1 (MQ=255)
                                                                                                           ccAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAtc                                >  1:1718904/1‑148 (MQ=255)
                                                                                                           ccAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAtc                                >  2:1689936/1‑148 (MQ=255)
                                                                                                           ccAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAtc                                <  2:3283117/148‑1 (MQ=255)
                                                                                                             aaCGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAtctc                              >  1:3493174/1‑148 (MQ=255)
                                                                                                              aCGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTcc                             <  2:2255009/148‑1 (MQ=255)
                                                                                                               cGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCca                            >  2:2141505/1‑148 (MQ=255)
                                                                                                                     ccAAACATTCGATTCATATCGCTATTGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCACGAGGATCACCAGCAAACGCCAACGTAACGGGATCTCGCCGTGAGATATCGAGCCAAGCAACGCGTCCAGCAACTCCACAGGTg                      >  2:1665882/1‑148 (MQ=255)
                                                                                                                        aaCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCg                   >  1:2533114/1‑148 (MQ=255)
                                                                                                                         aCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGt                  >  2:825924/1‑148 (MQ=255)
                                                                                                                         aCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGt                  >  1:820817/1‑148 (MQ=255)
                                                                                                                         aCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTACTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGt                  >  2:2666941/1‑148 (MQ=255)
                                                                                                                          cATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGtc                 >  2:2744520/1‑148 (MQ=255)
                                                                                                                            ttCGATTCATATCGCTATGGAAATAACGTTTCCCTCGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGtctc               <  1:1232064/148‑1 (MQ=255)
                                                                                                                             tCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGAATCTCCACAGGTGCCGTCTCa              >  2:1559018/1‑148 (MQ=255)
                                                                                                                              cGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCAt             <  2:2872714/148‑1 (MQ=255)
                                                                                                                                aTTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTa           >  2:127422/1‑148 (MQ=255)
                                                                                                                                 ttCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc          >  2:1885409/1‑148 (MQ=255)
                                                                                                                                 ttCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc          >  2:108616/1‑148 (MQ=255)
                                                                                                                                  tCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCATCATCTCCACAGGTGCCGTCTCATTAcc         >  1:1033442/1‑148 (MQ=255)
                                                                                                                                  tCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAcc         >  2:1633475/1‑148 (MQ=255)
                                                                                                                                   cATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCg        >  1:1177668/1‑148 (MQ=255)
                                                                                                                                   cATATCGCTATGGAAATAACGTGTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCg        <  2:1883356/148‑1 (MQ=255)
                                                                                                                                    atatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgt       >  2:2213673/1‑148 (MQ=255)
                                                                                                                                    atatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgt       >  2:2454682/1‑148 (MQ=255)
                                                                                                                                     tatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtg      <  2:123568/148‑1 (MQ=255)
                                                                                                                                     tatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtg      >  1:2159591/1‑148 (MQ=255)
                                                                                                                                      atCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtgt     <  1:2265683/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
AGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATT  >  CP009273/1820744‑1821028

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAT  >  CP009273/1820742‑1821027
                                                                                                                                                  |                                                                                                                                           
GTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAAT                                                                                                                                          <  SRR19973065.2134833/150‑1 (MQ=60)
GTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAAT                                                                                                                                          <  SRR19973065.269588/150‑1 (MQ=60)
  AGATCACGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAA                                                                                                                                        <  SRR19973065.2097651/150‑1 (MQ=60)
   GATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAC                                                                                                                                       <  SRR19973065.2853067/150‑1 (MQ=60)
    ATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGTTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACG                                                                                                                                      >  SRR19973065.1896153/1‑150 (MQ=60)
       AAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTT                                                                                                                                   <  SRR19973065.1107152/150‑1 (MQ=60)
       AAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTT                                                                                                                                   <  SRR19973065.2119836/150‑1 (MQ=60)
        AGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTC                                                                                                                                  >  SRR19973065.3035953/1‑150 (MQ=60)
        AGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTC                                                                                                                                  >  SRR19973065.2035583/1‑150 (MQ=60)
         GGTGCCAGCGCACATATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCC                                                                                                                                 <  SRR19973065.1425018/150‑1 (MQ=60)
         GGTGCCAGCGCAATGATTTTTTGCCCTGTTCATAAAAATCTTCATGGCACTGTTCCATTTGGCGGGCCGTACAGGTTTCTCCGCTTTCAGCAAATTGCTGCCAACGACCGCCAAACATCCGATTCATATCGCTATGGAAATAACGTTTCC                                                                                                                                 <  SRR19973065.1424929/150‑1 (MQ=60)
           TGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCCGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCT                                                                                                                               <  SRR19973065.155466/150‑1 (MQ=60)
           TGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCT                                                                                                                               <  SRR19973065.3281154/150‑1 (MQ=60)
           TGCCAGCGCACAGATTCTTTGCCCTGGTAATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACCGGTTTCTCGGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCT                                                                                                                               <  SRR19973065.649859/150‑1 (MQ=60)
               AGCGCACAGTTTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCT                                                                                                                           <  SRR19973065.664393/150‑1 (MQ=60)
                GCGCACTGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTT                                                                                                                          <  SRR19973065.1104250/150‑1 (MQ=60)
                GCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTT                                                                                                                          <  SRR19973065.1104174/150‑1 (MQ=60)
                 CGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTC                                                                                                                         >  SRR19973065.2147755/1‑150 (MQ=60)
                  GCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCA                                                                                                                        <  SRR19973065.3080788/150‑1 (MQ=60)
                  GCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCA                                                                                                                        <  SRR19973065.2266756/150‑1 (MQ=60)
                  GCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCA                                                                                                                        <  SRR19973065.3079762/150‑1 (MQ=60)
                  GCAAAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCA                                                                                                                        <  SRR19973065.2266670/150‑1 (MQ=60)
                   CACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAG                                                                                                                       >  SRR19973065.2449896/1‑150 (MQ=60)
                   CACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAG                                                                                                                       <  SRR19973065.3149568/150‑1 (MQ=60)
                   CACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAG                                                                                                                       >  SRR19973065.550554/1‑150 (MQ=60)
                   CACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAG                                                                                                                       >  SRR19973065.674536/1‑150 (MQ=60)
                       GATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTTCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCA                                                                                                                   >  SRR19973065.3269292/1‑150 (MQ=60)
                         TTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACCGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGG                                                                                                                 <  SRR19973065.1247491/150‑1 (MQ=60)
                         TTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGG                                                                                                                 <  SRR19973065.1247472/150‑1 (MQ=60)
                            TTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGG                                                                                                              <  SRR19973065.531093/150‑1 (MQ=60)
                                     GTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAG                                                                                                     >  SRR19973065.2435226/1‑150 (MQ=60)
                                     GTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAG                                                                                                     >  SRR19973065.662487/1‑150 (MQ=60)
                                        ATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGAT                                                                                                  >  SRR19973065.195348/1‑150 (MQ=60)
                                          AAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCA                                                                                                <  SRR19973065.894739/150‑1 (MQ=60)
                                                  TCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAAAATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAC                                                                                        <  SRR19973065.942795/150‑1 (MQ=60)
                                                     AGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGC                                                                                     >  SRR19973065.3418673/1‑150 (MQ=60)
                                                       GCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCA                                                                                   >  SRR19973065.88867/1‑150 (MQ=60)
                                                       GCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCA                                                                                   >  SRR19973065.89488/1‑150 (MQ=60)
                                                       GCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCA                                                                                   >  SRR19973065.88862/1‑150 (MQ=60)
                                                        CACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAA                                                                                  >  SRR19973065.1195801/1‑150 (MQ=60)
                                                        CACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAA                                                                                  >  SRR19973065.13984/1‑150 (MQ=60)
                                                        CACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAA                                                                                  >  SRR19973065.14086/1‑150 (MQ=60)
                                                         ACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAAC                                                                                 >  SRR19973065.2612295/1‑150 (MQ=60)
                                                         ACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAAC                                                                                 >  SRR19973065.2612377/1‑150 (MQ=60)
                                                           TGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGT                                                                               <  SRR19973065.1130997/150‑1 (MQ=60)
                                                           TGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGT                                                                               <  SRR19973065.2891953/150‑1 (MQ=60)
                                                              TCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAAC                                                                            <  SRR19973065.768237/150‑1 (MQ=60)
                                                              TCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAAC                                                                            <  SRR19973065.768329/150‑1 (MQ=60)
                                                                 AGTTCGCGCGCCCGACAGTTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGG                                                                         >  SRR19973065.646641/1‑150 (MQ=60)
                                                                 AGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGG                                                                         >  SRR19973065.2269291/1‑150 (MQ=60)
                                                                         CGCCCGACAGGTTTCTCCGATTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCC                                                                 <  SRR19973065.2862579/150‑1 (MQ=60)
                                                                               ACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATAAGATTCATATCGCTATGGAAAAAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGACGTGAGA                                                           >  SRR19973065.3463604/1‑150 (MQ=60)
                                                                                CAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGAT                                                          <  SRR19973065.1843958/150‑1 (MQ=60)
                                                                                    TTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCG                                                      <  SRR19973065.2746860/150‑1 (MQ=60)
                                                                                         CCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCA                                                 >  SRR19973065.139277/1‑150 (MQ=60)
                                                                                            CTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGC                                              >  SRR19973065.1556640/1‑150 (MQ=60)
                                                                                                CAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACG                                          >  SRR19973065.491434/1‑150 (MQ=60)
                                                                                                   CAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGT                                       >  SRR19973065.2624504/1‑150 (MQ=60)
                                                                                                    AAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTC                                      >  SRR19973065.1563004/1‑150 (MQ=60)
                                                                                                    AAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTC                                      >  SRR19973065.432663/1‑150 (MQ=60)
                                                                                                    AAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTC                                      >  SRR19973065.1583721/1‑150 (MQ=60)
                                                                                                      ATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCA                                    <  SRR19973065.1556508/150‑1 (MQ=60)
                                                                                                      ATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCA                                    <  SRR19973065.2941278/150‑1 (MQ=60)
                                                                                                        AGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGC                                  >  SRR19973065.1645292/1‑150 (MQ=60)
                                                                                                        AGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGC                                  >  SRR19973065.1110837/1‑150 (MQ=60)
                                                                                                           TGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATC                               >  SRR19973065.930349/1‑150 (MQ=60)
                                                                                                           TGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATC                               >  SRR19973065.1730327/1‑150 (MQ=60)
                                                                                                            GCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCT                              >  SRR19973065.206879/1‑150 (MQ=60)
                                                                                                             CCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTC                             >  SRR19973065.3514411/1‑150 (MQ=60)
                                                                                                               AACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCA                           <  SRR19973065.1852867/150‑1 (MQ=60)
                                                                                                                ACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCAC                          <  SRR19973065.1436344/150‑1 (MQ=60)
                                                                                                                  GACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAG                        <  SRR19973065.550418/150‑1 (MQ=60)
                                                                                                                  GACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAG                        <  SRR19973065.2224463/150‑1 (MQ=60)
                                                                                                                        CAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCG                  >  SRR19973065.2549014/1‑150 (MQ=60)
                                                                                                                         AAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGT                 >  SRR19973065.826209/1‑150 (MQ=60)
                                                                                                                            CATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTC              <  SRR19973065.652433/150‑1 (MQ=60)
                                                                                                                            CATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTC              <  SRR19973065.2427590/150‑1 (MQ=60)
                                                                                                                              TTCGATTCATATCGCTATGGAAATAACGTTTCCCTCGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCAT            <  SRR19973065.1240129/150‑1 (MQ=60)
                                                                                                                                CGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTA          <  SRR19973065.1173621/150‑1 (MQ=60)
                                                                                                                                  ATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCATCATCTCCACAGGTGCCGTCTCATTACC        >  SRR19973065.1040284/1‑150 (MQ=60)
                                                                                                                                   TTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCG       >  SRR19973065.1185145/1‑150 (MQ=60)
                                                                                                                                     CATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTG     >  SRR19973065.2173473/1‑150 (MQ=60)
                                                                                                                                     CATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTG     >  SRR19973065.2173443/1‑150 (MQ=60)
                                                                                                                                        ATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAT  <  SRR19973065.2280002/150‑1 (MQ=60)
                                                                                                                                                  |                                                                                                                                           
GTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAT  >  CP009273/1820742‑1821027

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