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 F394 I0 R1 1342 215.0 7134094 96.2% 6862998 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,829,319 C→A 100% I85I (ATC→ATA ydjZ → TVP38/TMEM64 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,829,3190CA96.9% 593.2 / 7.7 196I85I (ATC→ATAydjZTVP38/TMEM64 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base A (88/102);  total (89/107)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
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 196 total aligned reads displayed.

TCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTG  >  CP009273/1829174‑1829463
                                                                                                                                                 |                                                                                                                                                
tCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAc                                                                                                                                                <  2:692348/148‑1 (MQ=255)
 cAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAcc                                                                                                                                               <  2:1834932/148‑1 (MQ=255)
 cAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAcc                                                                                                                                               <  2:1938257/148‑1 (MQ=255)
 cAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAcc                                                                                                                                               <  1:1928276/148‑1 (MQ=255)
 cAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAcc                                                                                                                                               <  2:2541982/148‑1 (MQ=255)
  aaTCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCt                                                                                                                                              <  2:2164183/148‑1 (MQ=255)
   aTCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCCACAGGGTATAGAACGCTTTATTCAGTCTTAAGGCGCACTGGCGGCAGTTGTCTCGTTCTTGGTGATGATTTTGCAGGCCAATGCTGCACCGCAACCTGCGTTTTTGATAACCtt                                                                                                                                             >  1:1099191/1‑148 (MQ=255)
   aTCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCtt                                                                                                                                             <  2:1393393/148‑1 (MQ=255)
   aTCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCtt                                                                                                                                             >  2:1272003/1‑148 (MQ=255)
     cGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTg                                                                                                                                           >  2:1799142/1‑148 (MQ=255)
     cGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTg                                                                                                                                           >  2:134007/1‑148 (MQ=255)
     cGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTg                                                                                                                                           >  2:132943/1‑148 (MQ=255)
     cGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGGTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTg                                                                                                                                           >  2:1799905/1‑148 (MQ=255)
      gCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGc                                                                                                                                          <  2:2906646/148‑1 (MQ=255)
       cAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCCGTTGTCTCGTTCTTGTTGATGTTTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTGTGATAACCTTTGcc                                                                                                                                         >  2:2727717/1‑148 (MQ=255)
       cAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAACTTTGcc                                                                                                                                         >  1:3199897/1‑148 (MQ=255)
         gCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCaa                                                                                                                                       <  1:91042/148‑1 (MQ=255)
            ttGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGc                                                                                                                                    >  2:2117996/1‑148 (MQ=255)
              gcggcgTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGt                                                                                                                                  <  2:486264/148‑1 (MQ=255)
                   gTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTg                                                                                                                             <  1:84754/148‑1 (MQ=255)
                   gTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTg                                                                                                                             <  1:603526/148‑1 (MQ=255)
                    tttGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGt                                                                                                                            >  2:1494871/1‑148 (MQ=255)
                     ttGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGtt                                                                                                                           <  1:2564284/148‑1 (MQ=255)
                           ccgTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGgcgc                                                                                                                     >  2:2306247/1‑148 (MQ=255)
                           ccgTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGGTCTTGTTGATGATTTTGCAGACCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGgggc                                                                                                                     >  1:2690003/1‑148 (MQ=255)
                             gTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGt                                                                                                                   <  2:360103/148‑1 (MQ=255)
                                 aCCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTg                                                                                                               >  2:845954/1‑148 (MQ=255)
                                 aCCAACAGGGTATAGAACGCTTTATTCAGTCATACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTg                                                                                                               <  1:249083/148‑1 (MQ=255)
                                  ccAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTgg                                                                                                              <  1:2833880/148‑1 (MQ=255)
                                   cAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTggg                                                                                                             >  1:1481973/1‑148 (MQ=255)
                                     aCAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggc                                                                                                           >  2:3524642/1‑148 (MQ=255)
                                     aCAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGAGTTCTGGggc                                                                                                           >  2:2523255/1‑148 (MQ=255)
                                        gggTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggcggc                                                                                                        >  1:2154496/1‑148 (MQ=255)
                                        gggTATAGAACGCTGTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggcggc                                                                                                        >  2:1427481/1‑148 (MQ=255)
                                         ggTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCt                                                                                                       <  2:972308/148‑1 (MQ=255)
                                           tataGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTg                                                                                                     >  1:1571307/1‑148 (MQ=255)
                                           tataGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTg                                                                                                     <  2:2042199/148‑1 (MQ=255)
                                           tataGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTg                                                                                                     <  2:540649/148‑1 (MQ=255)
                                           tataGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTg                                                                                                     >  2:125464/1‑148 (MQ=255)
                                             taGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCt                                                                                                   <  2:598454/148‑1 (MQ=255)
                                             taGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCt                                                                                                   <  2:1469391/148‑1 (MQ=255)
                                               gAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGt                                                                                                 >  2:2263599/1‑148 (MQ=255)
                                                 aCGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCg                                                                                               <  2:3258530/148‑1 (MQ=255)
                                                 aCGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCg                                                                                               >  2:1915519/1‑148 (MQ=255)
                                                 aCGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCg                                                                                               >  2:1915526/1‑148 (MQ=255)
                                                       tATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGAcc                                                                                         <  2:3185283/148‑1 (MQ=255)
                                                          tCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGt                                                                                      <  1:2235019/148‑1 (MQ=255)
                                                           cAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGtt                                                                                     >  1:3108677/1‑148 (MQ=255)
                                                           cAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGtt                                                                                     <  2:2125722/148‑1 (MQ=255)
                                                            aGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCGTGCTGTCGTGGACCAGGTc                                                                                    >  1:2988111/1‑148 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  >  2:2402793/1‑148 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  >  1:1776171/1‑148 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  >  1:3537772/1‑148 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  <  2:3551077/148‑1 (MQ=255)
                                                               cTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGGTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGAt                                                                                 >  1:1315672/1‑148 (MQ=255)
                                                                 tACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATgg                                                                               >  1:2596282/1‑148 (MQ=255)
                                                                    ggCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcg                                                                            <  1:1897564/148‑1 (MQ=255)
                                                                    ggCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcg                                                                            <  2:2424395/148‑1 (MQ=255)
                                                                    ggCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcg                                                                            <  2:2424701/148‑1 (MQ=255)
                                                                      cgcACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcggc                                                                          >  1:2220975/1‑148 (MQ=255)
                                                                      cgcACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcggc                                                                          >  1:864866/1‑148 (MQ=255)
                                                                            ggcggcAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGAGTCGCTGGTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcggcgcggcg                                                                    >  2:2501784/1‑148 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCtt                                                            >  1:3485000/1‑148 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCtt                                                            >  2:2717429/1‑148 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCtt                                                            <  1:2014938/148‑1 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCtt                                                            <  1:1088628/148‑1 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGAATGCCGGCGCGGCGCAGTGCtt                                                            >  1:1964491/1‑148 (MQ=255)
                                                                                     tGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCttt                                                           <  2:2567180/148‑1 (MQ=255)
                                                                                          cGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTTGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTAt                                                      >  1:1246614/1‑148 (MQ=255)
                                                                                          cGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTAt                                                      >  1:255435/1‑148 (MQ=255)
                                                                                          cGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTAt                                                      <  1:2994265/148‑1 (MQ=255)
                                                                                             tCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGc                                                   >  2:3511352/1‑148 (MQ=255)
                                                                                             tCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGc                                                   >  1:3169861/1‑148 (MQ=255)
                                                                                             tCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGc                                                   >  2:2440492/1‑148 (MQ=255)
                                                                                                  ttgATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCagag                                              >  1:2004034/1‑148 (MQ=255)
                                                                                                    gatgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTg                                            <  2:334540/148‑1 (MQ=255)
                                                                                                    gatgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTg                                            <  2:727333/148‑1 (MQ=255)
                                                                                                      tgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGAt                                          >  2:220396/1‑148 (MQ=255)
                                                                                                          tttGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGc                                      <  1:2144019/148‑1 (MQ=255)
                                                                                                          tttGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCGGGTTGGGGCGTGCTGGGGCGGCGTGGGGTGGTGGTCCCGGTCGATGGGCGGCGCGGCGCGGTGGTTTTTTTTAGCCCGAGGGGTGGGc                                      >  1:6235/1‑148 (MQ=255)
                                                                                                                ggCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGaa                                <  1:2402774/148‑1 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               <  1:1799905/148‑1 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               <  1:1799142/148‑1 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               >  2:1624418/1‑148 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               <  2:3503745/148‑1 (MQ=255)
                                                                                                                  ccATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCGGTGCTTTTTTATCGCCAGAGGGATGGGCCGCGCAgt                              >  1:1840177/1‑148 (MQ=255)
                                                                                                                    aTTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtgg                            >  2:2682979/1‑148 (MQ=255)
                                                                                                                     ttGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggt                           <  1:1069098/148‑1 (MQ=255)
                                                                                                                     ttGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggt                           <  1:331490/148‑1 (MQ=255)
                                                                                                                     ttGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggt                           <  1:332674/148‑1 (MQ=255)
                                                                                                                           cACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGaaaa                     <  2:1739307/148‑1 (MQ=255)
                                                                                                                             ccGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAAt                   >  2:2056314/1‑148 (MQ=255)
                                                                                                                                cTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTaa                <  2:3251787/148‑1 (MQ=255)
                                                                                                                                 tACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAAc               <  1:906044/148‑1 (MQ=255)
                                                                                                                                 tACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAAc               <  1:1157928/148‑1 (MQ=255)
                                                                                                                                    cTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCCGAGTGATGGGCCGCGCAGTGGTGGAAAAATTAAcccg            >  2:1206594/1‑146 (MQ=255)
                                                                                                                                      gCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATAGCCACAGTGATGGGCCGCGCAGGGGTGGTAAAAATAACAGGCa          >  1:2710720/1‑148 (MQ=255)
                                                                                                                                         tttttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCaaaa       <  2:2154496/148‑1 (MQ=255)
                                                                                                                                            ttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCg    <  2:39666/148‑1 (MQ=255)
                                                                                                                                            ttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCg    <  2:3169915/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTG  >  CP009273/1829174‑1829463

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
ATCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGA  >  CP009273/1829173‑1829468
                                                                                                                                                  |                                                                                                                                                     
ATCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACC                                                                                                                                                    <  SRR19973066.2290305/150‑1 (MQ=60)
  CAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCCACAGGGTATAGAACGCTTTATTCAGTCTTAAGGCGCACTGGCGGCAGTTGTCTCGTTCTTGGTGATGATTTTGCAGGCCAATGCTGCACCGCAACCTGCGTTTTTGATAACCTT                                                                                                                                                  >  SRR19973066.1107182/1‑150 (MQ=60)
  CAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTT                                                                                                                                                  <  SRR19973066.1941982/150‑1 (MQ=60)
  CAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTT                                                                                                                                                  <  SRR19973066.1941933/150‑1 (MQ=60)
     TCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGC                                                                                                                                               <  SRR19973066.2863384/150‑1 (MQ=60)
      CGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAACTTTGCC                                                                                                                                              >  SRR19973066.3221143/1‑150 (MQ=60)
         AGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAAT                                                                                                                                           <  SRR19973066.3528830/150‑1 (MQ=60)
         AGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGGTTTTGTTAACCTTTGCCAAT                                                                                                                                           >  SRR19973066.3336139/1‑150 (MQ=60)
          GCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATG                                                                                                                                          <  SRR19973066.92030/150‑1 (MQ=60)
          GCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATG                                                                                                                                          <  SRR19973066.207404/150‑1 (MQ=60)
                    GTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTT                                                                                                                                <  SRR19973066.1427236/150‑1 (MQ=60)
                    GTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTT                                                                                                                                <  SRR19973066.607686/150‑1 (MQ=60)
                    GTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTT                                                                                                                                <  SRR19973066.85698/150‑1 (MQ=60)
                      TTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTG                                                                                                                              <  SRR19973066.2581643/150‑1 (MQ=60)
                      TTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTG                                                                                                                              <  SRR19973066.2581689/150‑1 (MQ=60)
                      TTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTG                                                                                                                              <  SRR19973066.2581706/150‑1 (MQ=60)
                          CGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGGTCTTGTTGATGATTTTGCAGACCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGGGC                                                                                                                          >  SRR19973066.2707958/1‑150 (MQ=60)
                                GGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTG                                                                                                                    <  SRR19973066.1837123/150‑1 (MQ=60)
                                  ACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGG                                                                                                                  >  SRR19973066.1492385/1‑150 (MQ=60)
                                  ACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGG                                                                                                                  <  SRR19973066.251007/150‑1 (MQ=60)
                                  ACCAACAGGGTATAGAACGCTTTATTCAGTCATACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGG                                                                                                                  <  SRR19973066.251002/150‑1 (MQ=60)
                                   CCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGG                                                                                                                 <  SRR19973066.2852927/150‑1 (MQ=60)
                                       CAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGC                                                                                                             >  SRR19973066.2170290/1‑150 (MQ=60)
                                       CAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGC                                                                                                             >  SRR19973066.2169457/1‑150 (MQ=60)
                                         GGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTT                                                                                                           <  SRR19973066.3359510/150‑1 (MQ=60)
                                          GGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTG                                                                                                          >  SRR19973066.1582637/1‑150 (MQ=60)
                                             ATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTG                                                                                                       <  SRR19973066.306730/150‑1 (MQ=60)
                                                          TTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTT                                                                                          >  SRR19973066.3129460/1‑150 (MQ=60)
                                                          TTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTT                                                                                          >  SRR19973066.3129441/1‑150 (MQ=60)
                                                           TCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTC                                                                                         <  SRR19973066.1953802/150‑1 (MQ=60)
                                                           TCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTC                                                                                         <  SRR19973066.2250439/150‑1 (MQ=60)
                                                           TCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTC                                                                                         >  SRR19973066.3007893/1‑150 (MQ=60)
                                                           TCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCGTGCTGTCGTGGACCAGGTC                                                                                         >  SRR19973066.3007876/1‑150 (MQ=60)
                                                             AGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGA                                                                                       >  SRR19973066.1788875/1‑150 (MQ=60)
                                                             AGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGA                                                                                       >  SRR19973066.3561287/1‑150 (MQ=60)
                                                              GTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGGTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGAT                                                                                      >  SRR19973066.1325211/1‑150 (MQ=60)
                                                                CTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGG                                                                                    >  SRR19973066.2613862/1‑150 (MQ=60)
                                                                     GGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGC                                                                               >  SRR19973066.871079/1‑150 (MQ=60)
                                                                     GGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGC                                                                               >  SRR19973066.2236338/1‑150 (MQ=60)
                                                                     GGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGC                                                                               >  SRR19973066.2238728/1‑150 (MQ=60)
                                                                     GGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGC                                                                               <  SRR19973066.1911100/150‑1 (MQ=60)
                                                                     GGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGC                                                                               <  SRR19973066.1911095/150‑1 (MQ=60)
                                                                             GGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCT                                                                       <  SRR19973066.914914/150‑1 (MQ=60)
                                                                               CGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGT                                                                     >  SRR19973066.1283670/1‑150 (MQ=60)
                                                                                GGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTG                                                                    >  SRR19973066.1101124/1‑150 (MQ=60)
                                                                                   AGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTT                                                                 <  SRR19973066.87056/150‑1 (MQ=60)
                                                                                   AGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTT                                                                 >  SRR19973066.3508249/1‑150 (MQ=60)
                                                                                   AGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTT                                                                 <  SRR19973066.84767/150‑1 (MQ=60)
                                                                                   AGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGAATGCCGGCGCGGCGCAGTGCTT                                                                 >  SRR19973066.1978375/1‑150 (MQ=60)
                                                                                     TTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTT                                                               <  SRR19973066.1096569/150‑1 (MQ=60)
                                                                                     TTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTT                                                               <  SRR19973066.2029036/150‑1 (MQ=60)
                                                                                     TTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTT                                                               <  SRR19973066.2029049/150‑1 (MQ=60)
                                                                                       GTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTT                                                             <  SRR19973066.625744/150‑1 (MQ=60)
                                                                                        TCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTA                                                            <  SRR19973066.1246059/150‑1 (MQ=60)
                                                                                         CTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTTGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTAGATGGCCGGCGCGGCGCTGTGCTTTTTTAT                                                           >  SRR19973066.1255559/1‑150 (MQ=60)
                                                                                         CTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTAT                                                           >  SRR19973066.257371/1‑150 (MQ=60)
                                                                                           CGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCG                                                         <  SRR19973066.3014118/150‑1 (MQ=60)
                                                                                           CGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCCGCGCTGTGCTTTTTTATAG                                                         >  SRR19973066.1993861/1‑150 (MQ=60)
                                                                                            GTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGC                                                        >  SRR19973066.3191030/1‑150 (MQ=60)
                                                                                            GTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGC                                                        >  SRR19973066.3190976/1‑150 (MQ=60)
                                                                                             TTCTTGTTGATGATTTTGCAGGCCAGTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCC                                                       <  SRR19973066.533478/150‑1 (MQ=60)
                                                                                                 TGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAG                                                   >  SRR19973066.2018103/1‑150 (MQ=60)
                                                                                                  GTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGT                                                  >  SRR19973066.1160133/1‑150 (MQ=60)
                                                                                                      ATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATG                                              <  SRR19973066.1524376/150‑1 (MQ=60)
                                                                                                         ATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCGGGTTGGGGCGTGCTGGGGCGGCgtggggtggtggtcccggtcgatgggcggcgcggcgcggtggttttttttagcccgaggggtgggc                                           >  SRR19973066.6345/1‑84 (MQ=60)
                                                                                                           TTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATACCCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCG                                         <  SRR19973066.2159060/150‑1 (MQ=60)
                                                                                                           TTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCG                                         <  SRR19973066.2158902/150‑1 (MQ=60)
                                                                                                                 GGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGT                                   <  SRR19973066.2419033/150‑1 (MQ=60)
                                                                                                                 GGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCGGTGCTTTTTTATCGCCAGAGGGATGGGCCGCGCAGT                                   >  SRR19973066.1853327/1‑150 (MQ=60)
                                                                                                                 GGCCAGTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGT                                   <  SRR19973066.2419014/150‑1 (MQ=60)
                                                                                                                  GCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTG                                  <  SRR19973066.1811930/150‑1 (MQ=60)
                                                                                                                  GCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTG                                  <  SRR19973066.1812698/150‑1 (MQ=60)
                                                                                                                  GCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTG                                  <  SRR19973066.1811952/150‑1 (MQ=60)
                                                                                                                      TTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGG                              <  SRR19973066.1076933/150‑1 (MQ=60)
                                                                                                                      TTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGG                              <  SRR19973066.335216/150‑1 (MQ=60)
                                                                                                                      TTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGG                              <  SRR19973066.334133/150‑1 (MQ=60)
                                                                                                                      TTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGG                              <  SRR19973066.334020/150‑1 (MQ=60)
                                                                                                                               CGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTAGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAA                     >  SRR19973066.2305685/1‑150 (MQ=60)
                                                                                                                                 CTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCGGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACC                   >  SRR19973066.1707363/1‑150 (MQ=60)
                                                                                                                                  TACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCG                  <  SRR19973066.1166156/150‑1 (MQ=60)
                                                                                                                                  TACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCG                  <  SRR19973066.912438/150‑1 (MQ=60)
                                                                                                                                  TACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCG                  <  SRR19973066.912441/150‑1 (MQ=60)
                                                                                                                                     CTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCA               >  SRR19973066.1790306/1‑150 (MQ=60)
                                                                                                                                     CTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATAGCCACAGTGATGGGCCGCGCAGGGGTGGTAAAAATAACAGGCA               >  SRR19973066.2728916/1‑150 (MQ=60)
                                                                                                                                     CTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCGGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCA               >  SRR19973066.2728866/1‑150 (MQ=60)
                                                                                                                                               GATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCT     >  SRR19973066.2070395/1‑150 (MQ=60)
                                                                                                                                               GATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCT     <  SRR19973066.2742901/150‑1 (MQ=60)
                                                                                                                                                  AACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGGTTTTTTTTAGCCAGAGGGATGGGCCGCGAAGaggggggaaaaataaaccgcaaaaacgggcttta  >  SRR19973066.20756/1‑116 (MQ=60)
                                                                                                                                                  |                                                                                                                                                     
ATCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGA  >  CP009273/1829173‑1829468

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: