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 1873 183.1 9376148 98.2% 9207377 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% 561.7 / NA 155Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (69/86);  total (69/86)
*NC_0009133,222,9450AG99.4% 484.8 / ‑4.5 155Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (69/85);  minor base T (0/1);  total (69/86)
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.41e-01
*NC_0009133,222,9460AC100.0% 490.1 / NA 155Q891S (CAA→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (69/86);  total (69/86)
*NC_0009133,222,9470GA99.4% 542.5 / NA 155A892T (GCC→ACC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (68/86);  minor base C (1/0);  total (69/86)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222863‑3223027
                                                                                 ||||                                                                                
gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                          >  1:106484/1‑93 (MQ=255)
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                         >  2:1193883/1‑93 (MQ=255)
  tCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCCCCGTCGATGc                                                                        >  2:2274211/1‑93 (MQ=255)
  tCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGc                                                                        <  2:2598767/93‑1 (MQ=255)
   cAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGcc                                                                       >  2:4031627/1‑93 (MQ=255)
   cAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGcc                                                                       <  2:171021/93‑1 (MQ=255)
   cAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGcc                                                                       <  2:4551444/93‑1 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGACAACTACGCCAACAGCCAGCAGGCTAACATCATCGATCTCTGGCGCAGCACCGTCGATGCCa                                                                      <  2:4687916/93‑1 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                      <  2:2818397/93‑1 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                      <  1:4108472/93‑1 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                     >  1:596564/1‑93 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                     <  2:103929/93‑1 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                     <  1:2675394/93‑1 (MQ=255)
       tGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                   >  2:1829077/1‑93 (MQ=255)
       tGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                   >  2:193016/1‑93 (MQ=255)
        ggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                  <  1:1724725/93‑1 (MQ=255)
        ggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                  <  2:2418937/93‑1 (MQ=255)
        ggCGTATTACGGTCGTGGACAGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                  <  1:1794791/93‑1 (MQ=255)
         gCGTATTCCGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                                                 <  2:28650/93‑1 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                               <  2:2118105/93‑1 (MQ=255)
            tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                              >  1:1849287/1‑93 (MQ=255)
            tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                              >  2:4482718/1‑93 (MQ=255)
             aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgaga                                                             >  2:451289/1‑93 (MQ=255)
             aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgaga                                                             >  2:451288/1‑93 (MQ=255)
             aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgaga                                                             >  2:3736075/1‑93 (MQ=255)
              ttACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                            <  2:4330055/93‑1 (MQ=255)
               tACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                                           <  1:2798304/93‑1 (MQ=255)
                 cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                         >  2:3231660/1‑93 (MQ=255)
                 cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                         <  2:3127317/93‑1 (MQ=255)
                  ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                        <  1:3713986/93‑1 (MQ=255)
                  ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                        >  2:2043927/1‑93 (MQ=255)
                   gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                       >  2:2357930/1‑93 (MQ=255)
                    tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                      <  1:1173469/93‑1 (MQ=255)
                     cGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                                     <  2:1107400/93‑1 (MQ=255)
                     cGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                                     <  1:1855206/93‑1 (MQ=255)
                      gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                    <  1:3190637/93‑1 (MQ=255)
                        ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCATCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                  <  1:1621430/93‑1 (MQ=255)
                        ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                  <  1:588135/93‑1 (MQ=255)
                        ggACCGGCCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                  <  2:833891/93‑1 (MQ=255)
                         gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                 >  2:2174858/1‑93 (MQ=255)
                              ggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGcgg                                            <  2:2445821/93‑3 (MQ=255)
                                cGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaa                                          <  1:351662/93‑1 (MQ=255)
                                 gAAAACTCCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                         <  2:3424128/93‑1 (MQ=255)
                                 gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                         >  1:796956/1‑93 (MQ=255)
                                 gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                         <  2:1451393/93‑1 (MQ=255)
                                 gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                         <  1:226170/93‑1 (MQ=255)
                                  aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                                        >  1:3790785/1‑93 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       >  2:2092541/1‑93 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       >  1:4267467/1‑93 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       <  1:2746624/93‑1 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       <  2:4210950/93‑1 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                      <  2:3935649/93‑1 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                      <  2:4233741/93‑1 (MQ=255)
                                      cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                    <  2:64402/93‑1 (MQ=255)
                                      cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                    >  1:4105675/1‑93 (MQ=255)
                                       caCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                   <  2:3335580/92‑1 (MQ=255)
                                         cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCTAGAACTATCCCTTCCCGCAGAACAACGGTaa                                 >  1:4486138/1‑93 (MQ=255)
                                         cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                 <  2:3309945/93‑1 (MQ=255)
                                          gCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                                <  1:401912/93‑1 (MQ=255)
                                          cccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                                <  1:2424751/92‑1 (MQ=255)
                                           ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                               <  1:18419/93‑1 (MQ=255)
                                            cGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCg                              >  1:972134/1‑93 (MQ=255)
                                             gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGt                             >  1:1469297/1‑93 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGTTAACCGTc                            <  2:4579253/93‑1 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                            >  2:3998594/1‑93 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                            <  1:50935/93‑1 (MQ=255)
                                                agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAg                          >  1:1958733/1‑93 (MQ=255)
                                                  ccagcAGTCTAACTTCATCGATAGCTGGCGCAGCCCCGTCGATGCGAGGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa                        <  2:2754523/93‑1 (MQ=255)
                                                  ccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa                        <  1:3471932/93‑1 (MQ=255)
                                                   cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt                       <  1:3934834/93‑1 (MQ=255)
                                                   cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt                       >  1:3116744/1‑93 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                      <  2:4096012/93‑1 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                      <  1:350724/93‑1 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                      >  1:3876648/1‑93 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                      <  2:2180655/93‑1 (MQ=255)
                                                     gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                     >  2:1428849/1‑93 (MQ=255)
                                                     gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                     >  1:1129658/1‑93 (MQ=255)
                                                      cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                    <  2:3886422/93‑1 (MQ=255)
                                                         gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                 <  2:490023/93‑1 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg               >  1:3840325/1‑93 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg              <  2:3610169/93‑1 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg              <  1:500073/93‑1 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa             >  2:3480225/1‑93 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa             >  1:4011762/1‑93 (MQ=255)
                                                               atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACg           <  2:673707/93‑1 (MQ=255)
                                                               atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACg           >  2:4337827/1‑93 (MQ=255)
                                                                 catcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGc         <  2:952283/93‑1 (MQ=255)
                                                                 catcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGc         <  2:3721092/93‑1 (MQ=255)
                                                                  atcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCa        >  2:1882939/1‑93 (MQ=255)
                                                                   tcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAc       <  1:2921648/93‑1 (MQ=255)
                                                                   tcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAc       <  1:2343022/93‑1 (MQ=255)
                                                                   tcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAc       <  1:2236208/93‑1 (MQ=255)
                                                                   tcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAc       <  2:928567/93‑1 (MQ=255)
                                                                      gtatCTGGCGCATCCCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGa    <  1:4586049/92‑1 (MQ=255)
                                                                      atatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGa    >  1:3266891/1‑93 (MQ=255)
                                                                      atatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGa    <  2:520508/93‑1 (MQ=255)
                                                                      atatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGa    <  1:2062096/93‑1 (MQ=255)
                                                                      atatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGa    <  1:1096581/93‑1 (MQ=255)
                                                                       tatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAc   <  1:2150346/93‑1 (MQ=255)
                                                                        atCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACt  <  2:1556196/93‑1 (MQ=255)
                                                                                 ||||                                                                                
GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222863‑3223027

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ 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
AACGGCGAATACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑CAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222851‑3223027
                                                                                               |                                                                                   
AACGGCGAATACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑A                                                                                     >  SRR4302535.4587281/1‑94 (MQ=60)
AACGGCGAATACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑A                                                                                     >  SRR4302535.849663/1‑94 (MQ=60)
             ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGC                                                                        <  SRR4302535.396632/94‑1 (MQ=60)
               CAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCA                                                                      <  SRR4302535.3528696/94‑1 (MQ=60)
                AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     <  SRR4302535.4136599/94‑1 (MQ=60)
                AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     >  SRR4302535.600744/1‑94 (MQ=60)
                AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCAT                                                                     >  SRR4302535.1351395/1‑94 (MQ=60)
                 GGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATG                                                                    <  SRR4302535.2694536/94‑1 (MQ=60)
                    GGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTC                                                                 <  SRR4302535.1737199/94‑1 (MQ=60)
                    GGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTC                                                                 <  SRR4302535.1855614/94‑1 (MQ=60)
                    GGCGTATTACGGTCGTGGACAGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTC                                                                 <  SRR4302535.1807765/94‑1 (MQ=60)
                       GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAG                                                              >  SRR4302535.1862614/1‑94 (MQ=60)
                       GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAG                                                              <  SRR4302535.3736793/94‑1 (MQ=60)
                         ATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAA                                                            >  SRR4302535.3831209/1‑94 (MQ=60)
                           TACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACT                                                          <  SRR4302535.2818217/94‑1 (MQ=60)
                             CGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTAT                                                        >  SRR4302535.4229347/1‑94 (MQ=60)
                              GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       <  SRR4302535.3739717/94‑1 (MQ=60)
                                TCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCC                                                     <  SRR4302535.1181846/94‑1 (MQ=60)
                                 CGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCT                                                    <  SRR4302535.1868572/94‑1 (MQ=60)
                                 CGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCT                                                    <  SRR4302535.522637/94‑1 (MQ=60)
                                  GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   >  SRR4302535.1895030/1‑94 (MQ=60)
                                  GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   <  SRR4302535.3212703/94‑1 (MQ=60)
                                    GGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCATCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCC                                                 <  SRR4302535.1633131/94‑1 (MQ=60)
                                    GGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCC                                                 <  SRR4302535.592257/94‑1 (MQ=60)
                                      ACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCG                                               <  SRR4302535.977228/94‑1 (MQ=60)
                                            CGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAAC                                         >  SRR4302535.802552/1‑94 (MQ=60)
                                            CGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAAC                                         <  SRR4302535.354087/94‑1 (MQ=60)
                                             GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        <  SRR4302535.227723/94‑1 (MQ=60)
                                             GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        >  SRR4302535.3816924/1‑94 (MQ=60)
                                              AAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAA                                       >  SRR4302535.4296687/1‑94 (MQ=60)
                                               AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      <  SRR4302535.2766203/94‑1 (MQ=60)
                                               AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      >  SRR4302535.2903694/1‑94 (MQ=60)
                                                 ACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGG                                    >  SRR4302535.4133786/1‑94 (MQ=60)
                                                 ACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGG                                    <  SRR4302535.3959462/94‑1 (MQ=60)
                                                    ACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCTAGAACTATCCCTTCCCGCAGAACAACGGTAA                                 >  SRR4302535.4516933/1‑94 (MQ=60)
                                                      GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               <  SRR4302535.404730/94‑1 (MQ=60)
                                                      CCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               <  SRR4302535.2442222/94‑1 (MQ=60)
                                                       CCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCG                              <  SRR4302535.18557/94‑1 (MQ=60)
                                                       CCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCG                              >  SRR4302535.979031/1‑94 (MQ=60)
                                                        CGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGT                             >  SRR4302535.1479822/1‑94 (MQ=60)
                                                          ACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCA                           <  SRR4302535.51309/94‑1 (MQ=60)
                                                           CAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAG                          >  SRR4302535.1972832/1‑94 (MQ=60)
                                                           CAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAG                          >  SRR4302535.1803931/1‑94 (MQ=60)
                                                             GCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCA                        >  SRR4302535.4073421/1‑94 (MQ=60)
                                                              CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       >  SRR4302535.3138442/1‑94 (MQ=60)
                                                              CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       <  SRR4302535.3496056/94‑1 (MQ=60)
                                                              CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       >  SRR4302535.4510394/1‑94 (MQ=60)
                                                               CAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG                      >  SRR4302535.3903274/1‑94 (MQ=60)
                                                               CAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG                      <  SRR4302535.3961683/94‑1 (MQ=60)
                                                                AGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGT                     <  SRR4302535.353141/94‑1 (MQ=60)
                                                                AGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGT                     >  SRR4302535.1137722/1‑94 (MQ=60)
                                                                 GCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTC                    <  SRR4302535.2748535/94‑1 (MQ=60)
                                                                  CAGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCC                   <  SRR4302535.1123128/94‑1 (MQ=60)
                                                                   AGGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCG                  >  SRR4302535.2633040/1‑94 (MQ=60)
                                                                    GGCTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGC                 >  SRR4302535.2281369/1‑94 (MQ=60)
                                                                      CTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               >  SRR4302535.3866789/1‑94 (MQ=60)
                                                                      CTAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               <  SRR4302535.1685348/94‑1 (MQ=60)
                                                                       TAACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGG              >  SRR4302535.474387/1‑94 (MQ=60)
                                                                        AACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGA             <  SRR4302535.503577/94‑1 (MQ=60)
                                                                        AACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGA             <  SRR4302535.2117461/94‑1 (MQ=60)
                                                                        AACATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGA             >  SRR4302535.4039207/1‑94 (MQ=60)
                                                                           ATCATCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGG          >  SRR4302535.3203935/1‑94 (MQ=60)
                                                                               TCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACT      <  SRR4302535.2359789/94‑1 (MQ=60)
                                                                               TCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACT      <  SRR4302535.2942258/94‑1 (MQ=60)
                                                                               TCGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACT      <  SRR4302535.2252116/94‑1 (MQ=60)
                                                                                CGATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTG     >  SRR4302535.335082/1‑94 (MQ=60)
                                                                                 GATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGA    >  SRR4302535.3289587/1‑94 (MQ=60)
                                                                                 GATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGA    >  SRR4302535.743692/1‑94 (MQ=60)
                                                                                  ATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC   >  SRR4302535.779809/1‑94 (MQ=60)
                                                                                  ATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC   <  SRR4302535.2076956/94‑1 (MQ=60)
                                                                                  ATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC   >  SRR4302535.141185/1‑94 (MQ=60)
                                                                                  ATATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC   <  SRR4302535.1104387/94‑1 (MQ=60)
                                                                                   TATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  <  SRR4302535.2165941/94‑1 (MQ=60)
                                                                                   TATCTGGCGCAGCA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  <  SRR4302535.2165810/94‑1 (MQ=60)
                                                                                               |                                                                                   
AACGGCGAATACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑CAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222851‑3223027

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