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 1061 294.1 15065038 98.6% 14854127 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,9440CA99.6% 926.5 / NA 242Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (114/127);  total (115/127)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,222,9450AG100.0% 818.9 / NA 242Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (115/127);  total (115/127)
*NC_0009133,222,9460AC100.0% 815.1 / NA 242Q891S (CAA→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (115/127);  total (115/127)
*NC_0009133,222,9470GA100.0% 889.4 / NA 242A892T (GCC→ACC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (115/127);  total (115/127)

Only 100 of 242 total aligned reads displayed.

GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACTA  >  NC_000913/3222863‑3223028
                                                                                 ||||                                                                                 
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                          >  2:5436519/1‑93 (MQ=255)
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                          >  1:6944985/1‑93 (MQ=255)
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                          >  1:6015069/1‑93 (MQ=255)
  tCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGc                                                                         <  1:4595635/93‑1 (MQ=255)
    nngTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                       <  1:1907262/91‑1 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                       >  1:5976126/1‑93 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                       >  2:5894351/1‑93 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                       <  2:2345292/93‑1 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                       <  1:4486832/93‑1 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                      <  2:3789946/93‑1 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                      >  1:5761379/1‑93 (MQ=255)
       tGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                    >  1:7109621/1‑93 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                                >  1:6355235/1‑93 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                                <  1:4045320/93‑1 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                                >  2:4446764/1‑93 (MQ=255)
            tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                               >  2:6121852/1‑93 (MQ=255)
            tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                               >  1:5683573/1‑93 (MQ=255)
            tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                               >  2:3893012/1‑93 (MQ=255)
             aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgaga                                                              <  2:1952828/93‑1 (MQ=255)
              ttACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                             <  2:249027/93‑1 (MQ=255)
              ttACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                             <  2:4057069/93‑1 (MQ=255)
               tACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                                            <  1:421581/93‑1 (MQ=255)
               tACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                                            <  2:1237175/93‑1 (MQ=255)
               tACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                                            <  2:5065284/93‑1 (MQ=255)
                 cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                          <  2:3013209/93‑1 (MQ=255)
                 cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                          <  1:3718646/93‑1 (MQ=255)
                 cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                          <  2:5998668/93‑1 (MQ=255)
                  ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                         <  2:4309939/93‑1 (MQ=255)
                   gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                        <  2:608629/93‑1 (MQ=255)
                   gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                        <  2:4983429/93‑1 (MQ=255)
                    tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                       >  1:5468104/1‑93 (MQ=255)
                    tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATACCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                       >  1:5666034/1‑93 (MQ=255)
                     cGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                                      <  2:6656137/93‑1 (MQ=255)
                     cGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                                      <  2:6396759/93‑1 (MQ=255)
                      gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                     <  2:283141/93‑1 (MQ=255)
                      gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                     <  1:2384243/93‑1 (MQ=255)
                      gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                     <  1:3835071/93‑1 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                    >  2:620597/1‑93 (MQ=255)
                          aCCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTccc                                                 <  2:5406558/93‑1 (MQ=255)
                             gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                                              >  1:1168475/1‑93 (MQ=255)
                               gCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGa                                            >  2:2199081/1‑93 (MQ=255)
                                cGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaa                                           <  2:3464599/93‑1 (MQ=255)
                                 gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                          <  1:5213440/93‑1 (MQ=255)
                                 gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                          >  2:5737160/1‑93 (MQ=255)
                                 gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                          <  2:1943343/93‑1 (MQ=255)
                                 gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                          >  1:4391556/1‑93 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                        >  1:2847093/1‑93 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                       <  1:2016468/93‑1 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                       <  1:5243904/93‑1 (MQ=255)
                                     aCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                                      <  2:6238988/93‑1 (MQ=255)
                                      cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                     >  1:1309975/1‑93 (MQ=255)
                                       tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                    >  1:4849313/1‑93 (MQ=255)
                                       tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                    <  2:2234812/93‑1 (MQ=255)
                                        aCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                                   <  1:333944/93‑1 (MQ=255)
                                        aCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                                   <  2:1608705/93‑1 (MQ=255)
                                         cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                  <  2:2605164/93‑1 (MQ=255)
                                         cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                  <  2:2742334/93‑1 (MQ=255)
                                          gCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                                 <  1:3914721/93‑1 (MQ=255)
                                           ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                                >  1:937364/1‑93 (MQ=255)
                                           ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                                <  2:7193361/93‑1 (MQ=255)
                                           ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                                <  2:7481521/93‑1 (MQ=255)
                                             gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGt                              >  1:1774836/1‑93 (MQ=255)
                                             gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGt                              >  2:5810021/1‑93 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                             >  2:880678/1‑93 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                             >  2:1058663/1‑93 (MQ=255)
                                               cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa                            >  2:7128189/1‑93 (MQ=255)
                                               cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa                            <  2:6814854/93‑1 (MQ=255)
                                               cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa                            >  2:6433255/1‑93 (MQ=255)
                                                  ccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa                         <  2:6977065/93‑1 (MQ=255)
                                                   cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt                        <  1:547663/93‑1 (MQ=255)
                                                   cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt                        >  1:1880421/1‑93 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                       >  2:183540/1‑93 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                       >  2:4667101/1‑93 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                       >  1:3891662/1‑93 (MQ=255)
                                                     gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                      <  1:5158575/93‑1 (MQ=255)
                                                       agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTcc                    >  1:4326969/1‑93 (MQ=255)
                                                       agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTcc                    <  1:7197364/93‑1 (MQ=255)
                                                         gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                  <  1:5019136/93‑1 (MQ=255)
                                                         gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                  >  2:4706944/1‑93 (MQ=255)
                                                         gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                  <  2:1563159/93‑1 (MQ=255)
                                                          cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                 <  1:6912612/93‑1 (MQ=255)
                                                          cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                 <  1:686193/93‑1 (MQ=255)
                                                          cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                 >  2:5972143/1‑93 (MQ=255)
                                                          cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                 <  1:1539155/93‑1 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg                <  2:3346041/93‑1 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg                >  1:6771604/1‑93 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg                >  1:871003/1‑93 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg               >  1:1439512/1‑93 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa              <  2:5114575/93‑1 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa              >  1:3543682/1‑93 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa              <  2:7283525/93‑1 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa              <  2:799133/93‑1 (MQ=255)
                                                              catcatCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGAc             >  2:1334617/1‑93 (MQ=255)
                                                               atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACg            <  2:2956114/93‑1 (MQ=255)
                                                               atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACg            >  2:7377614/1‑93 (MQ=255)
                                                               atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACg            <  1:1001025/93‑1 (MQ=255)
                                                                  atcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCa         >  2:4039383/1‑93 (MQ=255)
                                                                  atcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCa         >  2:3936849/1‑93 (MQ=255)
                                                                   tcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAc        >  1:2584596/1‑93 (MQ=255)
                                                                        atCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACt   <  2:3510374/93‑1 (MQ=255)
                                                                                 ||||                                                                                 
GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACTA  >  NC_000913/3222863‑3223028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ 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                                                                        <  SRR4302542.7495830/94‑1 (MQ=60)
 TCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCC                                                                       <  SRR4302542.4640573/94‑1 (MQ=60)
  CAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCA                                                                      >  SRR4302542.6033961/1‑94 (MQ=60)
   NNGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     <  SRR4302542.1923743/94‑1 (MQ=60)
   AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     <  SRR4302542.1486371/94‑1 (MQ=60)
   AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     <  SRR4302542.1977788/94‑1 (MQ=60)
   AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     <  SRR4302542.4530867/94‑1 (MQ=60)
   AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     <  SRR4302542.7227560/94‑1 (MQ=60)
   AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     >  SRR4302542.5817306/1‑94 (MQ=60)
     GTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGT                                                                   >  SRR4302542.7179407/1‑94 (MQ=60)
      TGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTT                                                                  <  SRR4302542.6994973/94‑1 (MQ=60)
        GCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCG                                                                <  SRR4302542.552184/94‑1 (MQ=60)
         CGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGA                                                               >  SRR4302542.6417833/1‑94 (MQ=60)
          GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAG                                                              >  SRR4302542.1398746/1‑94 (MQ=60)
          GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAG                                                              >  SRR4302542.920946/1‑94 (MQ=60)
          GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAG                                                              >  SRR4302542.5738761/1‑94 (MQ=60)
          GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAG                                                              <  SRR4302542.4085231/94‑1 (MQ=60)
          GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAG                                                              <  SRR4302542.5977813/94‑1 (MQ=60)
           TATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGA                                                             >  SRR4302542.6959918/1‑94 (MQ=60)
            ATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAA                                                            <  SRR4302542.3821874/94‑1 (MQ=60)
            ATTACGGTCGCGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAA                                                            >  SRR4302542.4688927/1‑94 (MQ=60)
             TTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAAC                                                           >  SRR4302542.3940912/1‑94 (MQ=60)
             TTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAAC                                                           >  SRR4302542.5395766/1‑94 (MQ=60)
              TACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACT                                                          <  SRR4302542.4930551/94‑1 (MQ=60)
              TACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACT                                                          <  SRR4302542.425200/94‑1 (MQ=60)
               ACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTA                                                         >  SRR4302542.3391231/1‑94 (MQ=60)
                CGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTAT                                                        <  SRR4302542.3755436/94‑1 (MQ=60)
                 GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       <  SRR4302542.6418783/94‑1 (MQ=60)
                 GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       <  SRR4302542.2068004/94‑1 (MQ=60)
                  GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCC                                                      >  SRR4302542.5521221/1‑94 (MQ=60)
                  GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATACCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCC                                                      >  SRR4302542.5721070/1‑94 (MQ=60)
                   TCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCC                                                     <  SRR4302542.3777037/94‑1 (MQ=60)
                   TCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCC                                                     <  SRR4302542.1091298/94‑1 (MQ=60)
                    CGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCT                                                    >  SRR4302542.4375841/1‑94 (MQ=60)
                    CGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCT                                                    <  SRR4302542.278635/94‑1 (MQ=60)
                    CGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCT                                                    >  SRR4302542.6872152/1‑94 (MQ=60)
                     GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   <  SRR4302542.2406277/94‑1 (MQ=60)
                     GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   <  SRR4302542.6005799/94‑1 (MQ=60)
                     GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   >  SRR4302542.1914997/1‑94 (MQ=60)
                     GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTA                                                   <  SRR4302542.3872999/94‑1 (MQ=60)
                      TGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTC                                                  <  SRR4302542.4307240/94‑1 (MQ=60)
                        GACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCC                                                <  SRR4302542.3290673/94‑1 (MQ=60)
                           CGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCA                                             >  SRR4302542.1178566/1‑94 (MQ=60)
                            GGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAG                                            >  SRR4302542.3653666/1‑94 (MQ=60)
                               CGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAAC                                         >  SRR4302542.4434644/1‑94 (MQ=60)
                                GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        <  SRR4302542.5264140/94‑1 (MQ=60)
                                GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        >  SRR4302542.789151/1‑94 (MQ=60)
                                 AAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAA                                       >  SRR4302542.2874461/1‑94 (MQ=60)
                                  AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      >  SRR4302542.4871803/1‑94 (MQ=60)
                                  AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      <  SRR4302542.2491934/94‑1 (MQ=60)
                                   AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     <  SRR4302542.5294877/94‑1 (MQ=60)
                                   AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     >  SRR4302542.3285290/1‑94 (MQ=60)
                                   AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     <  SRR4302542.2033885/94‑1 (MQ=60)
                                    ACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGG                                    >  SRR4302542.1321263/1‑94 (MQ=60)
                                    ACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGG                                    <  SRR4302542.4094291/94‑1 (MQ=60)
                                     CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   >  SRR4302542.2500985/1‑94 (MQ=60)
                                     CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   >  SRR4302542.4896534/1‑94 (MQ=60)
                                       ACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAA                                 <  SRR4302542.404148/94‑1 (MQ=60)
                                       ACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAA                                 <  SRR4302542.336843/94‑1 (MQ=60)
                                        CGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAC                                >  SRR4302542.5292377/1‑94 (MQ=60)
                                         GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               <  SRR4302542.448839/94‑1 (MQ=60)
                                         GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               >  SRR4302542.945439/1‑94 (MQ=60)
                                         GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               <  SRR4302542.3953350/94‑1 (MQ=60)
                                         GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               <  SRR4302542.5210645/94‑1 (MQ=60)
                                          CCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCG                              <  SRR4302542.5511647/94‑1 (MQ=60)
                                           CGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGT                             <  SRR4302542.136852/94‑1 (MQ=60)
                                           CGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGT                             >  SRR4302542.1790219/1‑94 (MQ=60)
                                           CGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGT                             >  SRR4302542.5777839/1‑94 (MQ=60)
                                             ACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCA                           <  SRR4302542.6933509/94‑1 (MQ=60)
                                             ACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCA                           <  SRR4302542.4313582/94‑1 (MQ=60)
                                                 CCAGCAGGCTAACATCGTCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       <  SRR4302542.1648574/94‑1 (MQ=60)
                                                 CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       <  SRR4302542.1006723/94‑1 (MQ=60)
                                                 CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       <  SRR4302542.6063665/94‑1 (MQ=60)
                                                 CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       >  SRR4302542.1896699/1‑94 (MQ=60)
                                                 CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       >  SRR4302542.3932584/1‑94 (MQ=60)
                                                  CAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG                      <  SRR4302542.552389/94‑1 (MQ=60)
                                                  CAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG                      <  SRR4302542.388279/94‑1 (MQ=60)
                                                  CAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG                      >  SRR4302542.3930082/1‑94 (MQ=60)
                                                    GCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTC                    <  SRR4302542.5208706/94‑1 (MQ=60)
                                                     CAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCC                   >  SRR4302542.5733275/1‑94 (MQ=60)
                                                     CAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCC                   >  SRR4302542.4369506/1‑94 (MQ=60)
                                                     CAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCC                   >  SRR4302542.5245867/1‑94 (MQ=60)
                                                      AGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCG                  <  SRR4302542.7267992/94‑1 (MQ=60)
                                                      AGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCG                  >  SRR4302542.1941805/1‑94 (MQ=60)
                                                        GCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCT                <  SRR4302542.5067878/94‑1 (MQ=60)
                                                        GCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCT                >  SRR4302542.614627/1‑94 (MQ=60)
                                                         CTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               >  SRR4302542.878512/1‑94 (MQ=60)
                                                         CTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               >  SRR4302542.6838190/1‑94 (MQ=60)
                                                         CTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               <  SRR4302542.1552446/94‑1 (MQ=60)
                                                         CTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               <  SRR4302542.692119/94‑1 (MQ=60)
                                                         CTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               <  SRR4302542.6980474/94‑1 (MQ=60)
                                                          TAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGG              <  SRR4302542.5607272/94‑1 (MQ=60)
                                                          TAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGG              >  SRR4302542.1451946/1‑94 (MQ=60)
                                                           AACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGA             >  SRR4302542.3576786/1‑94 (MQ=60)
                                                             CATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACG           <  SRR4302542.564363/94‑1 (MQ=60)
                                                             CATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACG           <  SRR4302542.7153769/94‑1 (MQ=60)
                                                              ATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGG          <  SRR4302542.1009639/94‑1 (MQ=60)
                                                                 ATCTATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAC       >  SRR4302542.444992/1‑94 (MQ=60)
                                                                 ATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAC       >  SRR4302542.5102080/1‑94 (MQ=60)
                                                                 ATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAC       >  SRR4302542.2609631/1‑94 (MQ=60)
                                                                  TCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACT      <  SRR4302542.7021510/94‑1 (MQ=60)
                                                                  TCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACT      <  SRR4302542.1627531/94‑1 (MQ=60)
                                                                   CGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTG     <  SRR4302542.1791761/94‑1 (MQ=60)
                                                                     ATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC   >  SRR4302542.4432474/1‑94 (MQ=60)
                                                                     ATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC   >  SRR4302542.3630203/1‑94 (MQ=60)
                                                                      TATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGGCGGCACTGACT  <  SRR4302542.6213336/94‑1 (MQ=60)
                                                                      TATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  <  SRR4302542.4474623/94‑1 (MQ=60)
                                                                      TATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  SRR4302542.112706/1‑94 (MQ=60)
                                                                                  |                                                                                   
ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑CAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222864‑3223027

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