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 I2 R1 882 254.4 12957588 98.0% 12698436 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,222,944 4 bp→AGCA coding (2671‑2674/3093 nt) b3076 → evolved beta‑D‑galactosidase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,222,9440CA100.0% 705.7 / NA 198Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (93/105);  total (93/105)
*NC_0009133,222,9450AG100.0% 680.6 / NA 198Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (93/105);  total (93/105)
*NC_0009133,222,9460AC99.5% 674.4 / NA 198Q891S (CAA→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (92/105);  total (93/105)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,222,9470GA99.5% 685.7 / NA 198A892T (GCC→ACC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (93/104);  total (93/105)
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.95e-01

Only 100 of 202 total aligned reads displayed.

TACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222860‑3223027
                                                                                    ||||                                                                                
tACGATTAGATCGCGTATTATGGTCGTGGCCCGGGCGAAAACTACGCCGCAAGCAAGCAGGCTAACCTCATCGATATCTGGCGCAGCGCCGTc                                                                             <  2:3210236/93‑1 (MQ=255)
   gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                          >  2:4685836/1‑93 (MQ=255)
   gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                          >  1:1918728/1‑93 (MQ=255)
    aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                         >  1:2373192/1‑93 (MQ=255)
     tCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGc                                                                        >  1:1303677/1‑93 (MQ=255)
      cAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGcc                                                                       <  2:2701889/93‑1 (MQ=255)
       aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTACCATCATCGATATCTGGCGCAGCACCGCCGATGCCa                                                                      >  1:888963/1‑93 (MQ=255)
       aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGcct                                                                      >  2:813860/1‑92 (MQ=255)
       aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                      <  1:1825780/93‑1 (MQ=255)
        ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                     <  1:2583898/93‑1 (MQ=255)
        ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGAACCGTCGATGCCAt                                                                     >  1:4222167/1‑93 (MQ=255)
        ggTGGCGTATGACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTACCATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                     <  2:220477/93‑1 (MQ=255)
         gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                    <  1:3676854/93‑1 (MQ=255)
         gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                    >  2:1566613/1‑93 (MQ=255)
          tGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGGCATGt                                                                   <  1:5868322/93‑1 (MQ=255)
          tGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                   >  1:1413634/1‑93 (MQ=255)
           ggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                  >  2:5617204/1‑93 (MQ=255)
               tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                              >  2:3556712/1‑93 (MQ=255)
               tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                              >  2:287988/1‑93 (MQ=255)
               tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                              >  2:1021880/1‑93 (MQ=255)
                aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgaga                                                             >  2:2127411/1‑93 (MQ=255)
                aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgaga                                                             <  1:1523797/93‑1 (MQ=255)
                 ttACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                            <  2:1819650/93‑1 (MQ=255)
                 ttACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                            >  2:6051568/1‑93 (MQ=255)
                    cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                         <  2:5488630/93‑1 (MQ=255)
                    cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                         <  1:1935740/93‑1 (MQ=255)
                     ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                        <  1:3396707/93‑1 (MQ=255)
                     ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                        >  2:5953273/1‑93 (MQ=255)
                     ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                        >  1:1273527/1‑93 (MQ=255)
                      gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                       >  1:3096791/1‑93 (MQ=255)
                      gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                       <  2:6186209/93‑1 (MQ=255)
                       tCGTGGACCGGGCGAACACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                      <  2:2798339/93‑1 (MQ=255)
                       gcGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                      <  1:4797247/92‑1 (MQ=255)
                       gcGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                      <  2:599386/92‑1 (MQ=255)
                        cGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                                     <  2:5881986/93‑1 (MQ=255)
                          tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                   >  1:5185875/1‑93 (MQ=255)
                            gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                 >  1:3721487/1‑93 (MQ=255)
                             aCCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTccc                                                <  1:4168907/93‑1 (MQ=255)
                               cGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                                              >  2:3067482/1‑93 (MQ=255)
                               cGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                                              >  2:1882186/1‑93 (MQ=255)
                               cGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                                              <  1:2235003/93‑1 (MQ=255)
                                gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                                             <  2:1267211/93‑1 (MQ=255)
                                 ggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAg                                            >  1:6205632/1‑93 (MQ=255)
                                 ggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAg                                            >  2:3907487/1‑93 (MQ=255)
                                   cGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaa                                          <  1:2665510/93‑1 (MQ=255)
                                    gAAAACTACGCGGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                         <  2:3139748/93‑1 (MQ=255)
                                    gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                         <  1:1595757/93‑1 (MQ=255)
                                    gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                         >  2:1989194/1‑93 (MQ=255)
                                     aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                                        >  2:2683910/1‑93 (MQ=255)
                                      aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       <  1:3546277/93‑1 (MQ=255)
                                      aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       <  1:4889585/93‑1 (MQ=255)
                                       aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                      >  2:6234228/1‑93 (MQ=255)
                                       aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                      <  1:1733131/93‑1 (MQ=255)
                                       aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                      <  2:4258853/93‑1 (MQ=255)
                                       aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                      <  1:4006735/93‑1 (MQ=255)
                                        aCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                                     <  2:4496306/93‑1 (MQ=255)
                                         cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                    <  1:3630819/93‑1 (MQ=255)
                                         cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                    <  1:605856/93‑1 (MQ=255)
                                          tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                   >  2:1373639/1‑93 (MQ=255)
                                          tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                   <  2:2057265/93‑1 (MQ=255)
                                           aCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                                  >  1:4977287/1‑93 (MQ=255)
                                            cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                 <  1:1551854/93‑1 (MQ=255)
                                            cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                 <  1:3049877/93‑1 (MQ=255)
                                            cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                 <  1:1101633/93‑1 (MQ=255)
                                            cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                 <  1:609008/93‑1 (MQ=255)
                                              ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                               <  1:5627651/93‑1 (MQ=255)
                                               cGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCg                              <  2:4965246/93‑1 (MQ=255)
                                               cGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATACCTTCCCGCAGAACAACGGTAACCg                              >  2:5125784/1‑93 (MQ=255)
                                                gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACACCGGTAACCGt                             <  2:4237709/93‑1 (MQ=255)
                                                 aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCTCTTCCCGCAGAACAACGGTAACTGTc                            >  1:3574234/1‑93 (MQ=255)
                                                 aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                            <  2:2486559/93‑1 (MQ=255)
                                                 aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                            >  1:2116962/1‑93 (MQ=255)
                                                   agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGTGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAg                          <  2:5292275/93‑1 (MQ=255)
                                                   agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAg                          <  1:2367753/93‑1 (MQ=255)
                                                   agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAg                          <  1:1453446/93‑1 (MQ=255)
                                                   agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAg                          >  2:6162276/1‑93 (MQ=255)
                                                     ccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa                        <  1:735398/93‑1 (MQ=255)
                                                      cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt                       >  1:4540541/1‑93 (MQ=255)
                                                      cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt                       <  2:5236686/93‑1 (MQ=255)
                                                       agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                      <  2:4437050/93‑1 (MQ=255)
                                                       agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                      <  1:492806/93‑1 (MQ=255)
                                                        gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTTCCGCAGAACAAAGGTAACAGTCAGCATga                     >  1:5091669/1‑92 (MQ=255)
                                                        gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                     >  2:5450640/1‑93 (MQ=255)
                                                        gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                     <  2:4369813/93‑1 (MQ=255)
                                                        gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                     >  2:4579513/1‑93 (MQ=255)
                                                            gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                 >  1:832956/1‑93 (MQ=255)
                                                             cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGACCTATCCCTCCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                <  2:2277978/93‑1 (MQ=255)
                                                               aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg              <  1:6471338/93‑1 (MQ=255)
                                                                aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa             <  2:3763798/93‑1 (MQ=255)
                                                                  atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGTACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACg           >  2:6476229/1‑93 (MQ=255)
                                                                   tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACgg          <  2:4179343/93‑1 (MQ=255)
                                                                    catcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGc         <  1:2925682/93‑1 (MQ=255)
                                                                     atcGATATCTGGCGCAGCACCGTCGATGCCATGTTAGAGAACTATCCCTTCCCGCAGAACAACGGCAACCGTCAGCATGTCCGCTGGACGGCa        >  2:2342596/1‑93 (MQ=255)
                                                                      tcGATATCTGGCGCAGCACCGTCGATCCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAc       <  2:4137647/93‑1 (MQ=255)
                                                                         atatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGa    <  2:3402103/93‑1 (MQ=255)
                                                                         atatCTGGCGCAGCACCATCGATGCCATGTTCGAGAACTATCCCTCCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGa    <  1:6323239/93‑1 (MQ=255)
                                                                          tatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAc   <  2:896296/93‑1 (MQ=255)
                                                                          tatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAc   >  2:1720636/1‑93 (MQ=255)
                                                                          tatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGCACGGCACTGAc   >  1:5137686/1‑93 (MQ=255)
                                                                           atCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACt  <  1:4685836/93‑1 (MQ=255)
                                                                                    ||||                                                                                
TACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222860‑3223027

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑CAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222864‑3223027
                                                                                  |                                                                                   
ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGC                                                                        >  SRR4302537.5053421/1‑94 (MQ=60)
ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGC                                                                        <  SRR4302537.3646257/94‑1 (MQ=60)
ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGC                                                                        >  SRR4302537.1312742/1‑94 (MQ=60)
 TCAGGTGGCGTATTTCGGTCGTGGACCGGGCGAAACCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCC                                                                       <  SRR4302537.867802/94‑1 (MQ=60)
   AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     <  SRR4302537.1838504/94‑1 (MQ=60)
   AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑AGAACCGTCGATGCCAT                                                                     >  SRR4302537.4251577/1‑94 (MQ=60)
    GGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATG                                                                    <  SRR4302537.2954135/94‑1 (MQ=60)
    GGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATG                                                                    <  SRR4302537.1886122/94‑1 (MQ=60)
    GGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATG                                                                    <  SRR4302537.2601982/94‑1 (MQ=60)
     GTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGT                                                                   >  SRR4302537.1423473/1‑94 (MQ=60)
     GTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGT                                                                   >  SRR4302537.1552249/1‑94 (MQ=60)
     GTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGT                                                                   <  SRR4302537.3702721/94‑1 (MQ=60)
      TGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGGCATGTT                                                                  <  SRR4302537.5908582/94‑1 (MQ=60)
      TGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTT                                                                  <  SRR4302537.3920538/94‑1 (MQ=60)
       GGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTC                                                                 <  SRR4302537.3745076/94‑1 (MQ=60)
         CGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGA                                                               <  SRR4302537.1306054/94‑1 (MQ=60)
         CGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGA                                                               >  SRR4302537.2277411/1‑94 (MQ=60)
          GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAG                                                              >  SRR4302537.3641355/1‑94 (MQ=60)
          GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAG                                                              <  SRR4302537.5008218/94‑1 (MQ=60)
           TATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGA                                                             <  SRR4302537.517207/94‑1 (MQ=60)
            ATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAA                                                            <  SRR4302537.1534440/94‑1 (MQ=60)
              TACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACT                                                          >  SRR4302537.2187389/1‑94 (MQ=60)
               ACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTA                                                         <  SRR4302537.3313810/94‑1 (MQ=60)
                CGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTAT                                                        <  SRR4302537.1949202/94‑1 (MQ=60)
                CGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTAT                                                        <  SRR4302537.762952/94‑1 (MQ=60)
                CGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTAT                                                        >  SRR4302537.1282381/1‑94 (MQ=60)
                 GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       <  SRR4302537.3420766/94‑1 (MQ=60)
                 GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       <  SRR4302537.1954249/94‑1 (MQ=60)
                 GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       >  SRR4302537.3118402/1‑94 (MQ=60)
                  GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCC                                                      >  SRR4302537.1585058/1‑94 (MQ=60)
                  GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCC                                                      <  SRR4302537.5783232/94‑1 (MQ=60)
                  GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCC                                                      <  SRR4302537.2821430/94‑1 (MQ=60)
                   GCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCC                                                     <  SRR4302537.4830350/94‑1 (MQ=60)
                    CGTGGACCGGGCGAAACCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCT                                                    <  SRR4302537.4138830/94‑1 (MQ=60)
                    CGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCT                                                    >  SRR4302537.3570439/1‑94 (MQ=60)
                     GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   <  SRR4302537.669217/94‑1 (MQ=60)
                     GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   >  SRR4302537.5221545/1‑94 (MQ=60)
                      TGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTC                                                  <  SRR4302537.1942517/94‑1 (MQ=60)
                       GGACCGGGCGAAAACTACGCCGACAGCCGGCAGGCTAACATAATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCC                                                 >  SRR4302537.5187223/1‑94 (MQ=60)
                       GGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCC                                                 >  SRR4302537.5187224/1‑94 (MQ=60)
                       GGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCC                                                 >  SRR4302537.3747633/1‑94 (MQ=60)
                        GACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCC                                                <  SRR4302537.550070/94‑1 (MQ=60)
                         ACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCG                                               <  SRR4302537.4198026/94‑1 (MQ=60)
                         ACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCG                                               <  SRR4302537.1728499/94‑1 (MQ=60)
                           CGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCA                                             <  SRR4302537.2250768/94‑1 (MQ=60)
                           CGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCA                                             >  SRR4302537.2399916/1‑94 (MQ=60)
                            GGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAG                                            <  SRR4302537.4734585/94‑1 (MQ=60)
                            GGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAG                                            >  SRR4302537.6248396/1‑94 (MQ=60)
                               CGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAAC                                         <  SRR4302537.2684131/94‑1 (MQ=60)
                                GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        <  SRR4302537.1606907/94‑1 (MQ=60)
                                GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        >  SRR4302537.5945503/1‑94 (MQ=60)
                                GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        >  SRR4302537.3093626/1‑94 (MQ=60)
                                  AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      <  SRR4302537.4923351/94‑1 (MQ=60)
                                  AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      <  SRR4302537.3571329/94‑1 (MQ=60)
                                   AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     >  SRR4302537.4838066/1‑94 (MQ=60)
                                   AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     <  SRR4302537.1745223/94‑1 (MQ=60)
                                   AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     <  SRR4302537.4034777/94‑1 (MQ=60)
                                   AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     <  SRR4302537.422928/94‑1 (MQ=60)
                                     CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   <  SRR4302537.3656404/94‑1 (MQ=60)
                                     CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   <  SRR4302537.1456808/94‑1 (MQ=60)
                                     CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   <  SRR4302537.609996/94‑1 (MQ=60)
                                      TACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTA                                  >  SRR4302537.5011592/1‑94 (MQ=60)
                                        CGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAC                                <  SRR4302537.613170/94‑1 (MQ=60)
                                        CGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAC                                <  SRR4302537.1139231/94‑1 (MQ=60)
                                        CGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAC                                <  SRR4302537.1109246/94‑1 (MQ=60)
                                        CGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAC                                <  SRR4302537.1562680/94‑1 (MQ=60)
                                        CGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAC                                <  SRR4302537.3071128/94‑1 (MQ=60)
                                          CCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCG                              >  SRR4302537.4457037/1‑94 (MQ=60)
                                          CCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCG                              <  SRR4302537.5666150/94‑1 (MQ=60)
                                           CGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGT                             >  SRR4302537.2602534/1‑94 (MQ=60)
                                            GACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCTCTTCCCGCAGAACAACGGTAACTGTC                            >  SRR4302537.3599444/1‑94 (MQ=60)
                                            GACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTC                            >  SRR4302537.2131873/1‑94 (MQ=60)
                                               AGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGC                         <  SRR4302537.1463579/94‑1 (MQ=60)
                                               AGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGC                         <  SRR4302537.2384442/94‑1 (MQ=60)
                                                 CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       >  SRR4302537.4571969/1‑94 (MQ=60)
                                                 CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       <  SRR4302537.740415/94‑1 (MQ=60)
                                                 CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       >  SRR4302537.4571970/1‑94 (MQ=60)
                                                   AGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTTCCGCAGAACAAAGGTAACAGTCAGCATGA                     >  SRR4302537.5126727/1‑94 (MQ=60)
                                                   AGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGT                     <  SRR4302537.496202/94‑1 (MQ=60)
                                                    GCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTC                    >  SRR4302537.6337110/1‑94 (MQ=60)
                                                     CAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCC                   <  SRR4302537.2385332/94‑1 (MQ=60)
                                                     CAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCC                   <  SRR4302537.3661404/94‑1 (MQ=60)
                                                       GGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGC                 >  SRR4302537.838654/1‑94 (MQ=60)
                                                        GCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCT                <  SRR4302537.1828132/94‑1 (MQ=60)
                                                        GCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCT                <  SRR4302537.6104992/94‑1 (MQ=60)
                                                         CTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               >  SRR4302537.1501136/1‑94 (MQ=60)
                                                          TAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGC              >  SRR4302537.1343187/1‑94 (MQ=60)
                                                           AACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGA             <  SRR4302537.6515217/94‑1 (MQ=60)
                                                              ATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGG          >  SRR4302537.6172218/1‑94 (MQ=60)
                                                                CATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCA        <  SRR4302537.2946229/94‑1 (MQ=60)
                                                                 ATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAC       >  SRR4302537.6384696/1‑94 (MQ=60)
                                                                 ATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAC       >  SRR4302537.4953021/1‑94 (MQ=60)
                                                                   CGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTG     >  SRR4302537.2694793/1‑94 (MQ=60)
                                                                    GATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGA    >  SRR4302537.2045118/1‑94 (MQ=60)
                                                                     ATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGCACGGCACTGAC   >  SRR4302537.5173089/1‑94 (MQ=60)
                                                                      TATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  <  SRR4302537.953831/94‑1 (MQ=60)
                                                                                  |                                                                                   
ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑CAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222864‑3223027

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