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 F1 I1 R1 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,149,302 T→C Q459R (CAG→CGG)  ppc ← phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,149,3020TC93.4% 661.2 / 24.4 212Q459R (CAG→CGG) ppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base C (86/112);  total (94/118)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

AGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGG  >  NC_000913/4149212‑4149391
                                                                                          |                                                                                         
agcagGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCgg                                                                                         >  1:2711100/1‑93 (MQ=255)
agcagGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCgg                                                                                         <  2:2669034/93‑1 (MQ=255)
 gcagGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGc                                                                                        <  1:2934913/93‑1 (MQ=255)
  cagGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCa                                                                                       <  2:2308664/93‑1 (MQ=255)
   agGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAg                                                                                      >  2:4154554/1‑93 (MQ=255)
   agGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAg                                                                                      <  2:363519/93‑1 (MQ=255)
    ggTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAgg                                                                                     >  1:3248279/1‑93 (MQ=255)
    ggTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAgg                                                                                     >  2:5739559/1‑93 (MQ=255)
     gTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGt                                                                                    <  2:4162512/93‑1 (MQ=255)
     gTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGt                                                                                    >  1:1781747/1‑93 (MQ=255)
      tGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTa                                                                                   <  1:5236565/93‑1 (MQ=255)
      tGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTa                                                                                   >  1:5797542/1‑93 (MQ=255)
        gACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATc                                                                                 <  1:4964535/93‑1 (MQ=255)
        gACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATc                                                                                 <  2:246838/93‑1 (MQ=255)
           aGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAg                                                                              <  1:1037793/93‑1 (MQ=255)
           aGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAg                                                                              >  1:1496759/1‑93 (MQ=255)
           aGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGACTTGCGGTGCTTCGGCAATCACCCGGCAGGTATGGAg                                                                              >  1:2356231/1‑93 (MQ=255)
            gCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGc                                                                             >  2:2713142/1‑93 (MQ=255)
            gCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGc                                                                             <  1:3317874/93‑1 (MQ=255)
            gCCAGTACGTCGGACGGCGTTTTCGCCATCGAAATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGc                                                                             <  1:5966870/93‑1 (MQ=255)
             ccAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCa                                                                            <  1:625120/93‑1 (MQ=255)
             ccAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCa                                                                            <  1:4248974/93‑1 (MQ=255)
              cAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCAc                                                                           <  2:1145931/93‑1 (MQ=255)
                gTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACtt                                                                         <  1:1608883/93‑1 (MQ=255)
                gTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACtt                                                                         <  2:2641272/93‑1 (MQ=255)
                 tACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTc                                                                        <  2:3298094/93‑1 (MQ=255)
                 tACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTc                                                                        >  1:4829318/1‑93 (MQ=255)
                    gTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTcgcg                                                                     >  2:3917501/1‑93 (MQ=255)
                    gTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGGATCGAGCACTTcgcg                                                                     >  2:2996814/1‑93 (MQ=255)
                     tCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTcgcgc                                                                    >  1:5897780/1‑93 (MQ=255)
                         gcGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGttt                                                                <  2:5801487/92‑1 (MQ=255)
                             cGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgc                                                            >  2:31591/1‑93 (MQ=255)
                             cGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgc                                                            <  2:5299473/93‑1 (MQ=255)
                              gTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcg                                                           <  2:776548/93‑1 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGTGTTTCGgcgc                                                          >  1:868184/1‑93 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcgc                                                          <  1:5424890/93‑1 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcgc                                                          <  2:2588015/93‑1 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcgc                                                          <  1:4163364/93‑1 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcgc                                                          <  1:1449939/93‑1 (MQ=255)
                                   cGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTgg                                                      <  2:1476009/93‑1 (MQ=255)
                                   cGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTgg                                                      >  1:3265720/1‑93 (MQ=255)
                                    gCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGt                                                     >  1:4549399/1‑93 (MQ=255)
                                    gCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGt                                                     <  2:2393900/93‑1 (MQ=255)
                                     ccATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGtt                                                    <  1:4524994/93‑1 (MQ=255)
                                     ccATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGtt                                                    <  2:5448643/93‑1 (MQ=255)
                                       aTCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGc                                                  >  2:6061557/1‑93 (MQ=255)
                                        tCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGcc                                                 >  1:2468374/1‑93 (MQ=255)
                                        tCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGcc                                                 >  2:5587808/1‑93 (MQ=255)
                                          gagaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAg                                               >  2:4496738/1‑93 (MQ=255)
                                          gagaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAg                                               <  1:4356494/93‑1 (MQ=255)
                                           agaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATTACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTCGTTGCCAGt                                              >  1:3888609/1‑93 (MQ=255)
                                            gaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCCGtt                                             >  2:1464710/1‑93 (MQ=255)
                                            gaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGtt                                             <  2:3067748/93‑1 (MQ=255)
                                             aTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTg                                            >  1:2052163/1‑93 (MQ=255)
                                             aTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTAGCGCGTTTCGGCGCTTGGTTGCCAGTTg                                            >  1:21956/1‑93 (MQ=255)
                                              tCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTgc                                           >  1:3107161/1‑93 (MQ=255)
                                              tCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTgc                                           >  1:138735/1‑93 (MQ=255)
                                               cACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTgcg                                          <  2:5228130/93‑1 (MQ=255)
                                                aCGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTgcgc                                         >  1:3048036/1‑93 (MQ=255)
                                                 cGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGGATCGAGCACTTCGCGCGTTTCGGCGCGTGGTTGCCAGTTGcgcg                                        >  2:3033165/1‑93 (MQ=255)
                                                  gTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGTCAGTTGCGCgg                                       >  1:2826195/1‑93 (MQ=255)
                                                   tAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGc                                      >  2:2742889/1‑93 (MQ=255)
                                                    aGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCa                                     >  2:3651228/1‑93 (MQ=255)
                                                       gggcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGaa                                  <  2:4756640/92‑1 (MQ=255)
                                                       cggcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGaa                                  <  2:1860964/93‑1 (MQ=255)
                                                        ggcAATGGAGCCTTGCGGTGCTTCGGCAATCGCCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAg                                 <  1:1863774/93‑1 (MQ=255)
                                                         gcAATGGAGCCTTGCGGTGCTTCGGCAATCCCCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGc                                <  2:212350/93‑1 (MQ=255)
                                                         gcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGc                                <  2:4696854/93‑1 (MQ=255)
                                                         gcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGc                                >  1:3970542/1‑93 (MQ=255)
                                                         gcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGc                                <  1:3729462/93‑1 (MQ=255)
                                                          cAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCg                               <  2:3370981/93‑1 (MQ=255)
                                                           aaTGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCgg                              <  2:6028164/93‑1 (MQ=255)
                                                            aTGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGa                             <  2:3749345/93‑1 (MQ=255)
                                                             tGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGAc                            >  2:1962859/1‑93 (MQ=255)
                                                             tGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGAc                            <  2:1594043/93‑1 (MQ=255)
                                                             tGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCACTTGGTTGCCAGTTGCGCGGCAGAAGCGGAc                            >  1:5693217/1‑93 (MQ=255)
                                                                   cTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTgg                      >  1:3226556/1‑93 (MQ=255)
                                                                   cTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTgg                      <  2:2851067/93‑1 (MQ=255)
                                                                         gTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttca                <  1:5050511/93‑1 (MQ=255)
                                                                         gTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttca                <  2:3463649/93‑1 (MQ=255)
                                                                         gTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGaggtca                >  1:2110887/1‑93 (MQ=255)
                                                                          tGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcag               <  1:1074191/93‑1 (MQ=255)
                                                                          tGCTTCGGCAATCACCCGGCAGGCATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcag               <  1:3483129/93‑1 (MQ=255)
                                                                            cTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcagtt             <  2:4977125/93‑1 (MQ=255)
                                                                            cTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcagtt             <  2:2196965/93‑1 (MQ=255)
                                                                            cTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcagtt             >  2:4904567/1‑93 (MQ=255)
                                                                             ttCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcagttc            <  1:745290/93‑1 (MQ=255)
                                                                               cGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTcgc          <  1:1242828/93‑1 (MQ=255)
                                                                                ggCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTcgcg         >  2:851401/1‑93 (MQ=255)
                                                                                ggCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTcgcg         <  2:3881955/93‑1 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      <  2:3923460/93‑1 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      >  2:1247798/1‑93 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      >  1:99048/1‑93 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      >  1:5111456/1‑93 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      <  1:1939374/93‑1 (MQ=255)
                                                                                    aTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATc     <  2:2800849/93‑1 (MQ=255)
                                                                                    aTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATc     <  2:1217086/93‑1 (MQ=255)
                                                                                    aTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATc     >  1:4021549/1‑93 (MQ=255)
                                                                                    aTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATc     <  1:1667407/93‑1 (MQ=255)
                                                                                     tCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGTGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCa    <  2:2159592/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGG  >  NC_000913/4149212‑4149391

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGCAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACG  >  NC_000913/4149209‑4149395
                                                                                             |                                                                                             
AGCAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCC                                                                                               >  SRR4302532.4278411/1‑94 (MQ=60)
AGCAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCC                                                                                               >  SRR4302532.549580/1‑94 (MQ=60)
AGCAGCAGGGGGACAGCCAGTACGTCGGAAGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCGTGCGGGGCGTCGGCAATCACCG                                                                                               >  SRR4302532.549579/1‑94 (MQ=60)
  CAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGG                                                                                             >  SRR4302532.498457/1‑94 (MQ=60)
  CAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGG                                                                                             >  SRR4302532.2802278/1‑94 (MQ=60)
   AGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGC                                                                                            >  SRR4302532.4274382/1‑94 (MQ=60)
    GCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCA                                                                                           <  SRR4302532.3032440/94‑1 (MQ=60)
     CAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAG                                                                                          <  SRR4302532.1420521/94‑1 (MQ=60)
     CAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAG                                                                                          <  SRR4302532.181627/94‑1 (MQ=60)
     CAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGCGCTTCGGCAATCACCCGGCAG                                                                                          <  SRR4302532.3359255/94‑1 (MQ=60)
      AGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGG                                                                                         >  SRR4302532.3356127/1‑94 (MQ=60)
      AGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGG                                                                                         >  SRR4302532.4692781/1‑94 (MQ=60)
       GGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGT                                                                                        <  SRR4302532.3345344/94‑1 (MQ=60)
       GGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGT                                                                                        >  SRR4302532.1841340/1‑94 (MQ=60)
        GTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTA                                                                                       >  SRR4302532.5999044/1‑94 (MQ=60)
        GTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTA                                                                                       >  SRR4302532.5999045/1‑94 (MQ=60)
         TGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTAT                                                                                      <  SRR4302532.5417133/94‑1 (MQ=60)
          GGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATC                                                                                     >  SRR4302532.2511749/1‑94 (MQ=60)
           GACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCG                                                                                    <  SRR4302532.5134786/94‑1 (MQ=60)
            ACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGA                                                                                   <  SRR4302532.403632/94‑1 (MQ=60)
             CAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAG                                                                                  >  SRR4302532.1547246/1‑94 (MQ=60)
             CAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGACTTGCGGTGCTTCGGCAATCACCCGGCAGGTATGGAG                                                                                  >  SRR4302532.2434843/1‑94 (MQ=60)
              AGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGC                                                                                 <  SRR4302532.2124288/94‑1 (MQ=60)
              AGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGC                                                                                 <  SRR4302532.1072801/94‑1 (MQ=60)
               GCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCA                                                                                <  SRR4302532.3428501/94‑1 (MQ=60)
               GCCAGTACGTCGGACGGCGTTTTCGCCATCGAAATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCA                                                                                <  SRR4302532.6173179/94‑1 (MQ=60)
                CCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCAC                                                                               <  SRR4302532.4391744/94‑1 (MQ=60)
                CCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCAC                                                                               <  SRR4302532.646961/94‑1 (MQ=60)
                   GTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTC                                                                            >  SRR4302532.4995571/1‑94 (MQ=60)
                   GTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTC                                                                            <  SRR4302532.5137778/94‑1 (MQ=60)
                   GTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTC                                                                            >  SRR4302532.3465841/1‑94 (MQ=60)
                   GTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTC                                                                            >  SRR4302532.199293/1‑94 (MQ=60)
                   GTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTC                                                                            <  SRR4302532.1663337/94‑1 (MQ=60)
                    TACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCTCCCGGCAGGTATCGAGCACTTCG                                                                           <  SRR4302532.5462989/94‑1 (MQ=60)
                      CGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCG                                                                         <  SRR4302532.1971705/94‑1 (MQ=60)
                      CGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCCGCAGGTATCGAGCACTTCGCG                                                                         >  SRR4302532.1603778/1‑94 (MQ=60)
                       GTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGC                                                                        >  SRR4302532.6102014/1‑94 (MQ=60)
                         CGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGT                                                                      >  SRR4302532.4568415/1‑94 (MQ=60)
                          GGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTT                                                                     <  SRR4302532.3443861/94‑1 (MQ=60)
                             CGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCG                                                                  <  SRR4302532.6318822/94‑1 (MQ=60)
                              GGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGG                                                                 <  SRR4302532.2584356/94‑1 (MQ=60)
                                CGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCG                                                               <  SRR4302532.4927029/94‑1 (MQ=60)
                                 GTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGTGTTTCGGCGC                                                              >  SRR4302532.897939/1‑94 (MQ=60)
                                  TTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCT                                                             <  SRR4302532.4302989/94‑1 (MQ=60)
                                  TTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCT                                                             <  SRR4302532.1498824/94‑1 (MQ=60)
                                  TTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCT                                                             <  SRR4302532.5612911/94‑1 (MQ=60)
                                   TTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTT                                                            <  SRR4302532.4366410/94‑1 (MQ=60)
                                     TCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGG                                                          <  SRR4302532.1544262/94‑1 (MQ=60)
                                     TCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGG                                                          >  SRR4302532.3374334/1‑94 (MQ=60)
                                      CGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGT                                                         >  SRR4302532.4704647/1‑94 (MQ=60)
                                        CCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTG                                                       <  SRR4302532.4679446/94‑1 (MQ=60)
                                        CCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTG                                                       <  SRR4302532.3379267/94‑1 (MQ=60)
                                          ATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCC                                                     >  SRR4302532.993215/1‑94 (MQ=60)
                                          ATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCC                                                     >  SRR4302532.2550550/1‑94 (MQ=60)
                                           TCGAGATCCCGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCA                                                    <  SRR4302532.3420941/94‑1 (MQ=60)
                                            CGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAG                                                   >  SRR4302532.681089/1‑94 (MQ=60)
                                             GAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATTACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTCGTTGCCAGT                                                  >  SRR4302532.4019817/1‑94 (MQ=60)
                                             GAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGT                                                  <  SRR4302532.4504011/94‑1 (MQ=60)
                                               GATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTG                                                >  SRR4302532.2120505/1‑94 (MQ=60)
                                               GATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTAGCGCGTTTCGGCGCTTGGTTGCCAGTTG                                                >  SRR4302532.22830/1‑94 (MQ=60)
                                                ATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGC                                               >  SRR4302532.143428/1‑94 (MQ=60)
                                                ATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGC                                               >  SRR4302532.3210222/1‑94 (MQ=60)
                                                 TCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTTGCGCGTTTCGGCGCTTGGTTGCCAGTTGCG                                              >  SRR4302532.523399/1‑94 (MQ=60)
                                                  CACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGC                                             >  SRR4302532.3148643/1‑94 (MQ=60)
                                                    CGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGTCAGTTGCGCGG                                           >  SRR4302532.2920726/1‑94 (MQ=60)
                                                       AGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCTGCGCTTGGTTGCCAGTTGCGCGGCAG                                        >  SRR4302532.1536361/1‑94 (MQ=60)
                                                       AGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAG                                        <  SRR4302532.5380393/94‑1 (MQ=60)
                                                           GGCAATGGAGCCTTGCGGTGCTTCGGCAATCGCCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGC                                    <  SRR4302532.1925532/94‑1 (MQ=60)
                                                           GGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGC                                    >  SRR4302532.4104300/1‑94 (MQ=60)
                                                           GGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGC                                    >  SRR4302532.6121107/1‑94 (MQ=60)
                                                            GCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCG                                   <  SRR4302532.3854753/94‑1 (MQ=60)
                                                            GCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCG                                   <  SRR4302532.3854752/94‑1 (MQ=60)
                                                              TATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCCCTTCGCGCGTTTCGGCCCTTGGTTGCCAGTTGCGCGGCAGAAGCGGA                                 <  SRR4302532.4018369/94‑1 (MQ=60)
                                                              AATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCCGAAGCGGT                                 >  SRR4302532.1081899/1‑94 (MQ=60)
                                                              AATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGA                                 >  SRR4302532.605234/1‑94 (MQ=60)
                                                               ATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGAC                                >  SRR4302532.5009932/1‑94 (MQ=60)
                                                               ATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCACTTGGTTGCCAGTTGCGCGGCAGAAGCGGAC                                >  SRR4302532.5891850/1‑94 (MQ=60)
                                                               ATGGAGCCTTGCGGTGCTTCGGCAATCACCCGCCAGGTATCGCGCACTTCACGCGGTTCGCTGCTTGCTTTCCCGTTTCGCGGCAGAAGCCGAC                                >  SRR4302532.2907520/1‑94 (MQ=25)
                                                                 GGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCTCTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGT                              <  SRR4302532.1680524/94‑1 (MQ=60)
                                                                    GCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTG                           <  SRR4302532.5864497/94‑1 (MQ=60)
                                                                     CCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGG                          >  SRR4302532.3333651/1‑94 (MQ=60)
                                                                      CTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGG                         >  SRR4302532.3099204/1‑94 (MQ=60)
                                                                      CTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGA                         >  SRR4302532.3171459/1‑94 (MQ=60)
                                                                      CTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGA                         <  SRR4302532.3328864/94‑1 (MQ=60)
                                                                      CTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGA                         <  SRR4302532.4143857/94‑1 (MQ=60)
                                                                       TTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAG                        <  SRR4302532.3934981/94‑1 (MQ=60)
                                                                        TGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGT                       >  SRR4302532.822394/1‑94 (MQ=60)
                                                                        TGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGT                       <  SRR4302532.430430/94‑1 (MQ=60)
                                                                         GCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTT                      <  SRR4302532.1448109/94‑1 (MQ=60)
                                                                           GGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCG                    >  SRR4302532.2350838/1‑94 (MQ=60)
                                                                           GGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCA                    <  SRR4302532.4442338/94‑1 (MQ=60)
                                                                           GGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGGTCA                    >  SRR4302532.2181347/1‑94 (MQ=60)
                                                                            GTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAG                   <  SRR4302532.5223813/94‑1 (MQ=60)
                                                                            GTGCGTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAG                   <  SRR4302532.3029893/94‑1 (MQ=60)
                                                                             TGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGT                  >  SRR4302532.1205232/1‑94 (MQ=60)
                                                                             TGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGT                  <  SRR4302532.1110250/94‑1 (MQ=60)
                                                                             TGCTTCGGCAATCACCCGGCAGGCATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGT                  <  SRR4302532.3599795/94‑1 (MQ=60)
                                                                               CTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTC                >  SRR4302532.5519048/1‑94 (MQ=60)
                                                                               CTTCGGCAATCACCCGGCAGGTATCGAGCACGTCGCGCGTTTCGGCGCTTTGTTGCCAGTTGCGCGGCAGCACCGGCCGTTTGGAGTTCAGTTC                >  SRR4302532.6155770/1‑94 (MQ=60)
                                                                                TTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCG               <  SRR4302532.770790/94‑1 (MQ=60)
                                                                                TTCGGCAATAACCCGGCAGGTATCTAGCACTGCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCG               <  SRR4302532.2438066/94‑1 (MQ=60)
                                                                                  CGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCG             <  SRR4302532.1284415/94‑1 (MQ=60)
                                                                                     GAATCACCCGGCAGGTATCGAGCACTTCGCTCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAT          <  SRR4302532.3192780/94‑1 (MQ=60)
                                                                                     CAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAT          >  SRR4302532.102484/1‑94 (MQ=60)
                                                                                     CAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAT          >  SRR4302532.5287209/1‑94 (MQ=60)
                                                                                     CAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAC          >  SRR4302532.2025213/1‑94 (MQ=60)
                                                                                      AATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATC         <  SRR4302532.2003866/94‑1 (MQ=60)
                                                                                      AATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATC         >  SRR4302532.4157072/1‑94 (MQ=60)
                                                                                       ATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCA        <  SRR4302532.1723844/94‑1 (MQ=60)
                                                                                       ATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCA        <  SRR4302532.2913521/94‑1 (MQ=60)
                                                                                        TCCCCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAG       <  SRR4302532.2834678/94‑1 (MQ=60)
                                                                                        GCACCCGGCAGGTATCGAGCACGTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAG       <  SRR4302532.2255258/94‑1 (MQ=60)
                                                                                          ACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGATTTCAGTTCGCGGATCAGGA     <  SRR4302532.6227757/94‑1 (MQ=60)
                                                                                           CCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAA    <  SRR4302532.6005923/94‑1 (MQ=60)
                                                                                            CCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAAC   >  SRR4302532.2168834/1‑94 (MQ=60)
                                                                                            CCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAAC   <  SRR4302532.3738295/94‑1 (MQ=60)
                                                                                             CGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACG  <  SRR4302532.4748498/94‑1 (MQ=60)
                                                                                             |                                                                                             
AGCAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACG  >  NC_000913/4149209‑4149395

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