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
A2 F1 I1 R1 404 265.1 10478538 96.8% 10143224 150.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 5,105,402 G→A F150F (TTC→TTT RPPX_22650 ← benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,105,4020GA91.9% 339.8 / 19.8 173F150F (TTC→TTTRPPX_22650benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base A (76/82);  total (83/90)
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.93e-01

Only 100 of 173 total aligned reads displayed.

GGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTTG  >  CP009974/5105264‑5105549
                                                                                                                                          |                                                                                                                                                   
agctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCAcc                                                                                                                                         <  2:3887493/144‑1 (MQ=1)
 gctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGACGCCCCGGTACGACTCGAATCGTGCCACCt                                                                                                                                        <  1:1196736/145‑1 (MQ=1)
  ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCtt                                                                                                                                       <  1:3582024/146‑1 (MQ=1)
  ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCtt                                                                                                                                       >  1:3533153/6‑151 (MQ=1)
  ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCtt                                                                                                                                       <  1:792928/146‑1 (MQ=1)
                       cgtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTc                                                                                                                  >  1:4543423/6‑151 (MQ=0)
                        gtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCg                                                                                                                 >  1:484598/5‑151 (MQ=0)
                         tccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGc                                                                                                                <  1:1442260/148‑1 (MQ=1)
                           gggggACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAg                                                                                                              <  2:3613312/147‑1 (MQ=0)
                                  tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                                                                       >  2:445610/1‑151 (MQ=11)
                                   gATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCt                                                                                                      >  2:3963779/1‑151 (MQ=11)
                                    aTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTg                                                                                                     >  1:1745460/1‑151 (MQ=11)
                                      cggcgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTc                                                                                                   <  1:2794826/148‑1 (MQ=11)
                                       gacgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCg                                                                                                  >  1:2962260/1‑151 (MQ=11)
                                         cgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCt                                                                                                      <  2:610701/145‑1 (MQ=11)
                                          gaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgggg                                                                                               >  1:4409847/1‑151 (MQ=11)
                                           aTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGt                                                                                              <  2:4176349/151‑1 (MQ=11)
                                           aTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGCGt                                                                                              >  2:1855308/1‑151 (MQ=11)
                                            tCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTa                                                                                             <  1:2460073/151‑1 (MQ=11)
                                            tCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTa                                                                                             <  1:1779292/151‑1 (MQ=11)
                                            tCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTa                                                                                             <  1:1057282/151‑1 (MQ=11)
                                             cATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAg                                                                                            >  2:1802813/1‑151 (MQ=11)
                                             cATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAg                                                                                            >  1:1465485/1‑151 (MQ=11)
                                              aTGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGc                                                                                           >  1:3658969/1‑151 (MQ=11)
                                               tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                                                                          <  2:1504569/151‑1 (MQ=11)
                                               tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                                                                          >  1:3363408/1‑151 (MQ=11)
                                               tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                                                                          >  1:3363511/1‑151 (MQ=11)
                                                gTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCg                                                                                         <  1:4467167/151‑1 (MQ=11)
                                                 tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                                                        <  2:2343252/151‑1 (MQ=11)
                                                 tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                                                        <  1:5035734/151‑1 (MQ=11)
                                                 tCGATGATCCTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                                                        <  2:4231015/151‑1 (MQ=11)
                                                    atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                                                     >  2:922196/1‑151 (MQ=11)
                                                    atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                                                     >  2:3041955/1‑151 (MQ=11)
                                                    atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                                                     >  1:2468429/1‑151 (MQ=11)
                                                    atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                                                     >  1:345807/1‑151 (MQ=11)
                                                     tgatCTTCGCCGCCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                                                    <  1:4716245/151‑1 (MQ=11)
                                                     tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                                                    <  1:162535/151‑1 (MQ=11)
                                                     tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                                                    <  1:4715976/151‑1 (MQ=11)
                                                            ggCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTgggg                                                                             <  2:1005725/150‑1 (MQ=11)
                                                            cGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTgggg                                                                             >  2:1180023/1‑151 (MQ=11)
                                                             gCCGGCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                                                            <  2:2213650/151‑1 (MQ=11)
                                                             gCCGAGTGGCGCGGGTGCGCGACCGGCGGCTTGACATCGGCGTTCAGGCGGCCGAAGAAAAAGCCCCGGTACGACTCAAAACGGGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                                                            <  2:4385034/151‑1 (MQ=2)
                                                             gCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                                                            <  2:2213272/151‑1 (MQ=11)
                                                             gCCGACTCGCCCAGGTGCGCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                                                            <  2:1596594/151‑1 (MQ=11)
                                                               cGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGCCCCGTCGCCGTTGAACCTGTCCGCGTCGCCGGCGTTGCTGGGGGCt                                                                          >  2:2743638/1‑151 (MQ=11)
                                                                gACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCtt                                                                         <  2:1280656/151‑1 (MQ=11)
                                                                 ccTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCttt                                                                        <  2:5070714/150‑1 (MQ=11)
                                                                  cTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTc                                                                       >  1:4064556/1‑151 (MQ=11)
                                                                      cccAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCt                                                                   <  1:3104213/151‑1 (MQ=11)
                                                                         aGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTca                                                                <  1:1275797/151‑1 (MQ=14)
                                                                          ggTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcag                                                               <  2:1578873/151‑1 (MQ=14)
                                                                            tGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcagca                                                             >  2:3927121/1‑151 (MQ=14)
                                                                             gCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcagcag                                                            >  1:4013846/1‑151 (MQ=14)
                                                                              cTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGt                                                           >  1:2374371/1‑151 (MQ=14)
                                                                               tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt                                                          >  2:1380660/1‑151 (MQ=14)
                                                                                  aCCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGc                                                       <  2:1922833/151‑1 (MQ=14)
                                                                                   ccAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGcc                                                      <  2:5075057/151‑1 (MQ=14)
                                                                                   ccAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGcc                                                      <  2:4544586/151‑1 (MQ=14)
                                                                                   ccAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGcc                                                      <  2:1729317/151‑1 (MQ=14)
                                                                                    cAGCGGCTTCGCATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCg                                                     <  2:3483636/151‑1 (MQ=14)
                                                                                       cggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtg                                                  <  1:3024940/147‑1 (MQ=14)
                                                                                        ggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgt                                                 >  1:4850394/4‑151 (MQ=14)
                                                                                          cttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgttg                                               <  2:2367358/150‑1 (MQ=14)
                                                                                           ttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgttgt                                              <  2:1234176/151‑1 (MQ=14)
                                                                                             cacaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGttgttg                                            >  2:3643176/1‑151 (MQ=14)
                                                                                             cacaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGttgttg                                            >  1:5202230/1‑151 (MQ=14)
                                                                                              acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa                                           <  1:2686384/151‑1 (MQ=14)
                                                                                              acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa                                           >  1:4830295/1‑151 (MQ=14)
                                                                                              acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa                                           >  2:379036/1‑151 (MQ=14)
                                                                                              acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa                                           >  1:4608444/1‑151 (MQ=14)
                                                                                              acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTGGTTGa                                           >  1:4579950/1‑151 (MQ=14)
                                                                                               caTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGaa                                          <  2:4282377/151‑1 (MQ=14)
                                                                                                  gggCGTTCAGGCTGCCGAGCAGAAAGCCCCGGGACGACTCAGAACGTGCCACGTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAGGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTGCGCCTTCAGCAGTTTGCCGCTGTTGTTGAACGt                                       <  2:4163975/150‑1 (MQ=14)
                                                                                                   ggCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTc                                      <  1:2287929/151‑1 (MQ=14)
                                                                                                    gcgttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTcc                                     <  2:578502/148‑1 (MQ=14)
                                                                                                    gcgttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTcc                                     >  1:3188074/4‑151 (MQ=14)
                                                                                                    gcgttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTcc                                     >  2:626267/4‑151 (MQ=14)
                                                                                                      gttCAGGCTGCCGAACGAAACGCCCCGGTACGACTCAACCCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAg                                   <  2:3266067/150‑1 (MQ=14)
                                                                                                      gttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAg                                   >  1:5133834/2‑151 (MQ=14)
                                                                                                      gttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAg                                   >  2:4511089/2‑151 (MQ=14)
                                                                                                       ttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGc                                  >  1:4969316/1‑151 (MQ=14)
                                                                                                        tCAGGCTGCCGAACAAAAAGCCCCGGTCCGACTCAAAACGTGCCACCTTGGTCAGGCCGCGGGCGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCCGTCTGCCGCTGCTGTTGAACCTCCAGcc                                 >  2:4643434/1‑151 (MQ=14)
                                                                                                          aggcTGCCGAACAAAGAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGt                               <  2:3681940/148‑1 (MQ=14)
                                                                                                              tGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGaa                           <  1:2005523/151‑1 (MQ=255)
                                                                                                                   aacaaaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGc                      <  1:2248149/145‑1 (MQ=255)
                                                                                                                   aacaaaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGc                      <  1:2154909/145‑1 (MQ=255)
                                                                                                                    acaaaaaGCCCCGGTACGACTCAAAACGTGCAACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACATTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCa                     <  1:634750/146‑1 (MQ=14)
                                                                                                                     caaaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAg                    <  2:2003300/147‑1 (MQ=255)
                                                                                                                       aaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGt                  >  2:2467424/3‑151 (MQ=255)
                                                                                                                        aaaGCCCCGGTGCGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGta                 <  2:3413220/150‑1 (MQ=255)
                                                                                                                         aaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGtat                >  2:4640279/1‑151 (MQ=255)
                                                                                                                          aGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGCAGCGCCAGGtata               >  2:5108679/1‑151 (MQ=17)
                                                                                                                              ccGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTg           >  1:843804/1‑151 (MQ=255)
                                                                                                                               cGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCTGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTgg          >  2:3347112/1‑151 (MQ=255)
                                                                                                                                ggTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGa         >  1:3735662/1‑151 (MQ=255)
                                                                                                                                ggTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGa         <  1:4311240/151‑1 (MQ=255)
                                                                                                                                 gTACGACTCAAGACGTGCGACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGaa        <  2:4884392/151‑1 (MQ=255)
                                                                                                                                 gTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGaa        >  1:2066213/1‑151 (MQ=255)
                                                                                                                                   aCGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACg      <  2:3364968/151‑1 (MQ=255)
                                                                                                                                       cTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTTg  <  2:678597/151‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
GGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTTG  >  CP009974/5105264‑5105549

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTTGC  >  CP009974/5105252‑5105550
                                                                                                                                                      |                                                                                                                                                    
GATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCG                                                                                                                                                      <  SRR10918231.4254325/151‑1 (MQ=0)
    taggagctagctccaCGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATC                                                                                                                                                  <  SRR10918231.1668704/136‑1 (MQ=0)
     ACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACG                                                                                                                                                 <  SRR10918231.3002559/151‑1 (MQ=0)
         gctagctccaCGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCC                                                                                                                                             <  SRR10918231.4928093/141‑1 (MQ=0)
         gctagctccaCGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCC                                                                                                                                             >  SRR10918231.4929732/11‑151 (MQ=0)
          ctagctccaCGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCA                                                                                                                                            <  SRR10918231.4252598/142‑1 (MQ=0)
           tagctccaCGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCAC                                                                                                                                           >  SRR10918231.1879884/9‑151 (MQ=0)
              ctccaCGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTT                                                                                                                                        >  SRR10918231.3549268/6‑151 (MQ=0)
              CGCCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTT                                                                                                                                        <  SRR10918231.4105927/151‑1 (MQ=0)
                CCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGG                                                                                                                                      >  SRR10918231.2738863/1‑151 (MQ=0)
                 CGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGT                                                                                                                                     >  SRR10918231.1586177/1‑151 (MQ=0)
                 CGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGT                                                                                                                                     >  SRR10918231.4699795/1‑151 (MQ=0)
                 CGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGT                                                                                                                                     >  SRR10918231.1021897/1‑151 (MQ=0)
                    GCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAG                                                                                                                                  >  SRR10918231.3019929/1‑151 (MQ=0)
                     CAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCTGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGG                                                                                                                                 >  SRR10918231.4726651/1‑151 (MQ=0)
                      AGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGT                                                                                                                                >  SRR10918231.69817/1‑151 (MQ=0)
                       GCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTC                                                                                                                               >  SRR10918231.2999228/1‑151 (MQ=0)
                          CTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGACAGGTCGTG                                                                                                                            >  SRR10918231.3363849/1‑151 (MQ=0)
                           TTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGG                                                                                                                           >  SRR10918231.3457699/1‑151 (MQ=0)
                                GGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCC                                                                                                                      >  SRR10918231.540539/1‑151 (MQ=0)
                                 GCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCG                                                                                                                     <  SRR10918231.4726588/151‑1 (MQ=0)
                                  CCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGT                                                                                                                    <  SRR10918231.3520127/151‑1 (MQ=0)
                                  CCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGT                                                                                                                    <  SRR10918231.2494679/151‑1 (MQ=0)
                                  CCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGTTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGGGCCGT                                                                                                                    >  SRR10918231.2219691/1‑151 (MQ=0)
                                   CTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTC                                                                                                                   >  SRR10918231.3761790/1‑151 (MQ=0)
                                    gtccGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCG                                                                                                                  >  SRR10918231.486339/5‑151 (MQ=0)
                                     TCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGC                                                                                                                 <  SRR10918231.1449014/151‑1 (MQ=0)
                                     TCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGC                                                                                                                 <  SRR10918231.1220014/151‑1 (MQ=0)
                                     TCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGC                                                                                                                 >  SRR10918231.1456101/1‑151 (MQ=0)
                                      CCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCA                                                                                                                >  SRR10918231.3738253/1‑151 (MQ=0)
                                      CCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCA                                                                                                                <  SRR10918231.3314588/151‑1 (MQ=0)
                                          CGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTG                                                                                                            >  SRR10918231.5052514/1‑151 (MQ=0)
                                          CGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTG                                                                                                            >  SRR10918231.4796656/1‑151 (MQ=0)
                                           GACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGA                                                                                                           >  SRR10918231.5105079/1‑151 (MQ=0)
                                                ATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTG                                                                                                      >  SRR10918231.1753468/1‑151 (MQ=4)
                                                  CGGCGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTC                                                                                                    <  SRR10918231.2807641/151‑1 (MQ=5)
                                                  CGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTC                                                                                                    <  SRR10918231.253729/151‑1 (MQ=5)
                                                   GACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCG                                                                                                   >  SRR10918231.2975842/1‑151 (MQ=5)
                                                     CGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTagatcg                                                                                                 >  SRR10918231.612973/1‑145 (MQ=2)
                                                     CGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGG                                                                                                 >  SRR10918231.2125521/1‑151 (MQ=5)
                                                      GATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGG                                                                                                >  SRR10918231.4430762/1‑151 (MQ=5)
                                                        TCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTA                                                                                              <  SRR10918231.1062060/151‑1 (MQ=5)
                                                        TCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTA                                                                                              <  SRR10918231.2471327/151‑1 (MQ=5)
                                                        TCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTA                                                                                              <  SRR10918231.1787412/151‑1 (MQ=5)
                                                         CATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAG                                                                                             >  SRR10918231.1472331/1‑151 (MQ=5)
                                                          ATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGC                                                                                            >  SRR10918231.3675706/1‑151 (MQ=5)
                                                           TGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCC                                                                                           >  SRR10918231.3378756/1‑151 (MQ=5)
                                                           TGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCC                                                                                           >  SRR10918231.3378652/1‑151 (MQ=5)
                                                            GTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCG                                                                                          <  SRR10918231.4488394/151‑1 (MQ=5)
                                                             TCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGG                                                                                         <  SRR10918231.5060232/151‑1 (MQ=5)
                                                                ATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTTGCAGTTGAAGCTGTCGGGGTAGCCGGCGT                                                                                      >  SRR10918231.347096/1‑151 (MQ=5)
                                                                ATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGT                                                                                      >  SRR10918231.347061/1‑151 (MQ=5)
                                                                ATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGT                                                                                      >  SRR10918231.2479722/1‑151 (MQ=5)
                                                                 TGATCTTCGCCGCCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTT                                                                                     <  SRR10918231.4738796/151‑1 (MQ=5)
                                                                 TGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTT                                                                                     <  SRR10918231.4738526/151‑1 (MQ=5)
                                                                 TGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTT                                                                                     <  SRR10918231.4365094/151‑1 (MQ=5)
                                                                 TGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTT                                                                                     <  SRR10918231.163126/151‑1 (MQ=5)
                                                                     CTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTG                                                                                 <  SRR10918231.2954971/151‑1 (MQ=7)
                                                                           CGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCT                                                                           >  SRR10918231.201905/1‑151 (MQ=17)
                                                                            GACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTT                                                                          <  SRR10918231.314629/151‑1 (MQ=17)
                                                                             ACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTT                                                                         >  SRR10918231.1783967/1‑151 (MQ=17)
                                                                              CTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTC                                                                        >  SRR10918231.4083662/1‑151 (MQ=19)
                                                                               TCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCA                                                                       >  SRR10918231.4599006/1‑151 (MQ=20)
                                                                                  CCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCT                                                                    <  SRR10918231.3118394/151‑1 (MQ=25)
                                                                                   CCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTT                                                                   <  SRR10918231.5000203/151‑1 (MQ=25)
                                                                                     AGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCA                                                                 <  SRR10918231.1281937/151‑1 (MQ=29)
                                                                                       GTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGC                                                               >  SRR10918231.2371720/1‑151 (MQ=32)
                                                                                         GCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAG                                                             >  SRR10918231.4032685/1‑151 (MQ=34)
                                                                                          CTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGT                                                            >  SRR10918231.2385190/1‑151 (MQ=35)
                                                                                            CGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGGGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCCGTTT                                                          >  SRR10918231.4666733/1‑151 (MQ=30)
                                                                                             GACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTG                                                         >  SRR10918231.3117157/1‑151 (MQ=40)
                                                                                              ACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGC                                                        <  SRR10918231.4921066/151‑1 (MQ=42)
                                                                                                CAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCG                                                      <  SRR10918231.3417311/151‑1 (MQ=43)
                                                                                                   cggcTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTG                                                   <  SRR10918231.3038788/147‑1 (MQ=45)
                                                                                                    GGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGT                                                  >  SRR10918231.4873725/1‑151 (MQ=44)
                                                                                                     GCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTT                                                 <  SRR10918231.4865998/151‑1 (MQ=44)
                                                                                                       TTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGT                                               >  SRR10918231.3646085/1‑151 (MQ=45)
                                                                                                         CACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTG                                             >  SRR10918231.5227901/1‑151 (MQ=45)
                                                                                                          ACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGA                                            >  SRR10918231.4853525/1‑151 (MQ=45)
                                                                                                          ACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGA                                            >  SRR10918231.4630475/1‑151 (MQ=45)
                                                                                                          ACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGA                                            <  SRR10918231.2698750/151‑1 (MQ=45)
                                                                                                          ACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTGGTTGA                                            >  SRR10918231.4601801/1‑151 (MQ=43)
                                                                                                           CATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAA                                           <  SRR10918231.667582/151‑1 (MQ=45)
                                                                                                             TCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACG                                         <  SRR10918231.843978/151‑1 (MQ=45)
                                                                                                               GGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTC                                       <  SRR10918231.2298247/151‑1 (MQ=45)
                                                                                                                GCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCC                                      >  SRR10918231.3202591/1‑151 (MQ=45)
                                                                                                                  gttcaggctgccgaacaaaAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAG                                    <  SRR10918231.1341402/132‑1 (MQ=53)
                                                                                                                  gttcaggctgccgaacaaaAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAG                                    >  SRR10918231.5159007/20‑151 (MQ=53)
                                                                                                                   TTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGC                                   >  SRR10918231.4993350/1‑151 (MQ=46)
                                                                                                                      gggctgccgaacaaaagGCCCCGGTACGTCTCAAATCGTGCCACCGTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGT                                <  SRR10918231.1279171/134‑1 (MQ=56)
                                                                                                                         CTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGA                             >  SRR10918231.1877754/1‑151 (MQ=57)
                                                                                                                          TGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAA                            <  SRR10918231.2014535/151‑1 (MQ=58)
                                                                                                                            ccgaacaaaAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACG                          <  SRR10918231.2082192/142‑1 (MQ=60)
                                                                                                                             cgaacaaaAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGG                         >  SRR10918231.4476659/9‑151 (MQ=60)
                                                                                                                               AACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGC                       <  SRR10918231.2258268/151‑1 (MQ=59)
                                                                                                                               AACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGC                       <  SRR10918231.2769706/151‑1 (MQ=59)
                                                                                                                               AACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGC                       <  SRR10918231.2164563/151‑1 (MQ=59)
                                                                                                                                acaaaAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCA                      <  SRR10918231.3259442/146‑1 (MQ=60)
                                                                                                                                acaaaAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCA                      >  SRR10918231.511810/6‑151 (MQ=60)
                                                                                                                                acaaaAAGCCCCGGTACGACTCAAAACGTGCAACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACATTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCA                      <  SRR10918231.637144/146‑1 (MQ=59)
                                                                                                                                 caaaAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAG                     <  SRR10918231.3726231/147‑1 (MQ=60)
                                                                                                                                  aaaAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGG                    >  SRR10918231.3920584/4‑151 (MQ=60)
                                                                                                                                     AAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTAT                 <  SRR10918231.4276624/151‑1 (MQ=60)
                                                                                                                                      AGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATA                >  SRR10918231.4628842/1‑151 (MQ=60)
                                                                                                                                          CCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTG            >  SRR10918231.847314/1‑151 (MQ=60)
                                                                                                                                            GGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGA          >  SRR10918231.3752806/1‑151 (MQ=60)
                                                                                                                                            GGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGA          <  SRR10918231.4331696/151‑1 (MQ=60)
                                                                                                                                             GTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAA         >  SRR10918231.2075462/1‑151 (MQ=60)
                                                                                                                                              TACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAAC        >  SRR10918231.2413653/1‑151 (MQ=60)
                                                                                                                                                 GACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGT     <  SRR10918231.444355/151‑1 (MQ=60)
                                                                                                                                                 GACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGT     <  SRR10918231.4739851/151‑1 (MQ=60)
                                                                                                                                                   CTCAAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTTG   <  SRR10918231.1601707/151‑1 (MQ=60)
                                                                                                                                                   CTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTTG   >  SRR10918231.2998412/1‑151 (MQ=60)
                                                                                                                                                    TCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTTGC  >  SRR10918231.2557127/1‑151 (MQ=60)
                                                                                                                                                      |                                                                                                                                                    
GATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTTGC  >  CP009974/5105252‑5105550

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