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 593 264.2 13311096 97.8% 13018251 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% 821.1 / NA 212Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (104/108);  total (104/108)
*NC_0009133,222,9450AG100.0% 768.2 / NA 212Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (104/108);  total (104/108)
*NC_0009133,222,9460AC100.0% 772.9 / NA 212Q891S (CAA→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (104/108);  total (104/108)
*NC_0009133,222,9470GA99.1% 800.0 / NA 212A892T (GCC→ACC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (104/106);  minor base C (0/2);  total (104/108)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC  >  NC_000913/3222863‑3223026
                                                                                 ||||                                                                               
gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                         <  2:4503572/93‑1 (MQ=255)
gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                         >  2:2318611/1‑93 (MQ=255)
gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                         <  1:837315/93‑1 (MQ=255)
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                        >  2:2479770/1‑93 (MQ=255)
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                        >  1:922717/1‑93 (MQ=255)
  tCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGc                                                                       <  2:58721/93‑1 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                    <  2:5375220/93‑1 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                    <  2:5331079/93‑1 (MQ=255)
      gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                   <  1:5241885/93‑1 (MQ=255)
       tGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                  >  1:6586243/1‑93 (MQ=255)
        ggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                 >  1:1140879/1‑93 (MQ=255)
         gCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                                                <  2:4315344/93‑1 (MQ=255)
         gCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                                                >  2:5820880/1‑93 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                              >  2:4415266/1‑93 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                              >  1:5681957/1‑93 (MQ=255)
            tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                             >  1:3843364/1‑93 (MQ=255)
             aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgaga                                                            <  1:1255107/93‑1 (MQ=255)
              ttACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                           >  2:5520019/1‑93 (MQ=255)
                 cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                        <  1:5865271/93‑1 (MQ=255)
                 cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                        <  2:3884893/93‑1 (MQ=255)
                 cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                        >  2:5315354/1‑93 (MQ=255)
                  ggTCGTGTACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCCCCGTCGATGCCATGTTCGAGAACTAt                                                       <  1:2194694/93‑1 (MQ=255)
                  ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                       >  1:869574/1‑93 (MQ=255)
                  ggTCGTGGACCGAGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                       <  2:1038068/93‑1 (MQ=255)
                  ggTCGTCGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                       <  1:4713948/93‑1 (MQ=255)
                   gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                      >  2:427497/1‑93 (MQ=255)
                   gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                      <  1:4894194/93‑1 (MQ=255)
                     cGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                                    <  2:394994/93‑1 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                  <  1:6515820/93‑1 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                  <  2:3149118/93‑1 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                  >  1:407925/1‑93 (MQ=255)
                         gcccGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                <  1:4908791/91‑1 (MQ=255)
                         gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                <  1:4985750/93‑1 (MQ=255)
                         gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                <  2:3199660/93‑1 (MQ=255)
                         gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                <  2:1122052/93‑1 (MQ=255)
                         gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                <  2:4930899/93‑1 (MQ=255)
                           ccGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                                              >  1:1513045/1‑93 (MQ=255)
                           ccGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                                              >  1:3922033/1‑93 (MQ=255)
                            cGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                                             <  2:1818452/93‑1 (MQ=255)
                             gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                                            <  1:718163/93‑1 (MQ=255)
                              ggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAg                                           >  1:2824239/1‑93 (MQ=255)
                                cGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaa                                         <  2:5470614/93‑1 (MQ=255)
                                  aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                                       >  1:4738232/1‑93 (MQ=255)
                                  aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                                       <  2:2075789/93‑1 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                      >  2:2980471/1‑93 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                     <  1:6274779/93‑1 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                     >  1:5622075/1‑93 (MQ=255)
                                     aCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                                    >  2:4144047/1‑93 (MQ=255)
                                       tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                  >  1:2414427/1‑93 (MQ=255)
                                       tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCATTCCCGCAGAACAACgga                                  >  1:1432314/1‑92 (MQ=255)
                                        aCGCTGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                                 <  1:6379367/93‑1 (MQ=255)
                                        aCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                                 >  2:2377416/1‑93 (MQ=255)
                                        aCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                                 <  2:6117546/93‑1 (MQ=255)
                                          gCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                               >  2:3738837/1‑93 (MQ=255)
                                          gCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                               <  2:6256847/93‑1 (MQ=255)
                                           ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                              >  1:2688474/1‑93 (MQ=255)
                                           ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                              >  1:3508510/1‑93 (MQ=255)
                                           ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                              >  2:1210944/1‑93 (MQ=255)
                                            cGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCg                             >  2:5197560/1‑93 (MQ=255)
                                             gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGt                            <  1:5334000/93‑1 (MQ=255)
                                             gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGt                            >  1:1825327/1‑93 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                           >  2:1542942/1‑93 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                           >  1:4752750/1‑93 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGGAACCGTc                           >  1:5813290/1‑93 (MQ=255)
                                               cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa                          <  2:2241614/93‑1 (MQ=255)
                                               cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa                          >  2:1394583/1‑93 (MQ=255)
                                               cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa                          >  1:4291619/1‑93 (MQ=255)
                                                agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAg                         <  1:2964047/93‑1 (MQ=255)
                                                 gccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGc                        >  2:2550901/1‑93 (MQ=255)
                                                 gccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGc                        >  2:1323963/1‑93 (MQ=255)
                                                  ccagcAGGCTAACATCATCGGTATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa                       >  1:5409371/1‑93 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                     >  2:3077155/1‑93 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                     <  2:5210387/93‑1 (MQ=255)
                                                     gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                    >  1:5064339/1‑93 (MQ=255)
                                                      cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                   >  2:1782334/1‑93 (MQ=255)
                                                       agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTcc                  >  2:3402798/1‑93 (MQ=255)
                                                       agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTcc                  <  2:3836201/93‑1 (MQ=255)
                                                        ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCg                 >  2:6359827/1‑93 (MQ=255)
                                                         gCTACCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                <  1:3478906/93‑1 (MQ=255)
                                                         gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                >  2:1235273/1‑93 (MQ=255)
                                                          cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt               >  2:6512794/1‑93 (MQ=255)
                                                          cTAACATCATCGATATCTGGCCCAGCATCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt               <  2:3529291/93‑1 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTTTCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg              >  1:2575001/1‑93 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg              <  2:2865833/93‑1 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCACGCAGAACAACGGTAACCGTCAGCATGTCCGCTg              >  2:3479939/1‑93 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg             <  1:3731027/93‑1 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg             >  1:2113753/1‑93 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg             <  2:2156940/93‑1 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa            <  1:398995/93‑1 (MQ=255)
                                                                tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACgg         >  2:4396305/1‑93 (MQ=255)
                                                                   tcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAc      >  1:4613164/1‑93 (MQ=255)
                                                                   tcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAc      >  2:4989160/1‑93 (MQ=255)
                                                                    cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACt     <  2:1712307/93‑1 (MQ=255)
                                                                    cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACt     <  1:2244291/93‑1 (MQ=255)
                                                                     gATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTg    >  2:597282/1‑93 (MQ=255)
                                                                      atatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGa   <  1:3180590/93‑1 (MQ=255)
                                                                      atatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGa   <  2:4650845/93‑1 (MQ=255)
                                                                       tatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAc  >  1:2532028/1‑93 (MQ=255)
                                                                       tatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAc  <  1:1893006/93‑1 (MQ=255)
                                                                       tatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAc  >  2:749458/1‑93 (MQ=255)
                                                                                 ||||                                                                               
GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC  >  NC_000913/3222863‑3223026

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑CAAG‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222866‑3223027
                                                                                |                                                                                     
CAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCA                                                                      >  SRR4302547.1219141/1‑94 (MQ=60)
CAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCA                                                                      >  SRR4302547.1562755/1‑94 (MQ=60)
   GTGGCGTATTACGGTCGTGGCCCGGGAGAAACCTACGCAAACAGCCCGCAGTCTAACATCATCGATATCTGGCGC‑‑‑‑AGCCCCGTCGAAGCCATGT                                                                   <  SRR4302547.401969/94‑1 (MQ=34)
   GTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGT                                                                   <  SRR4302547.5349784/94‑1 (MQ=60)
   GTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGT                                                                   >  SRR4302547.6723370/1‑94 (MQ=60)
   GTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGT                                                                   >  SRR4302547.2225366/1‑94 (MQ=60)
    TGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTT                                                                  >  SRR4302547.1163916/1‑94 (MQ=60)
    TGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTT                                                                  <  SRR4302547.632528/94‑1 (MQ=60)
     GGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTC                                                                 <  SRR4302547.857256/94‑1 (MQ=60)
       CGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGA                                                               >  SRR4302547.5799272/1‑94 (MQ=60)
        GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAG                                                              >  SRR4302547.3921625/1‑94 (MQ=60)
          ATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAA                                                            <  SRR4302547.1280417/94‑1 (MQ=60)
           TTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAAC                                                           <  SRR4302547.4847190/94‑1 (MQ=60)
             ACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTA                                                         >  SRR4302547.3561942/1‑94 (MQ=60)
              CGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTAT                                                        <  SRR4302547.5986663/94‑1 (MQ=60)
              CGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTAT                                                        >  SRR4302547.6488063/1‑94 (MQ=60)
              CGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTAT                                                        >  SRR4302547.886228/1‑94 (MQ=60)
               GGTCGTGTACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑‑‑AGCCCCGTCGATGCCATGTTCGAGAACTATC                                                       <  SRR4302547.2238765/94‑1 (MQ=60)
               GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       <  SRR4302547.480236/94‑1 (MQ=60)
               GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       >  SRR4302547.2553542/1‑94 (MQ=60)
               GGTCGTCGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       <  SRR4302547.4811536/94‑1 (MQ=60)
                GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCC                                                      <  SRR4302547.4995555/94‑1 (MQ=60)
                GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCC                                                      <  SRR4302547.6560472/94‑1 (MQ=60)
                  CGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCT                                                    <  SRR4302547.4347493/94‑1 (MQ=60)
                   GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   >  SRR4302547.3318136/1‑94 (MQ=60)
                   GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   >  SRR4302547.415492/1‑94 (MQ=60)
                   GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   >  SRR4302547.3963940/1‑94 (MQ=60)
                    TGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTC                                                  <  SRR4302547.2689328/94‑1 (MQ=60)
                    TGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTC                                                  <  SRR4302547.6651472/94‑1 (MQ=60)
                      GCCCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCC                                                <  SRR4302547.5010448/94‑1 (MQ=60)
                      GACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCC                                                <  SRR4302547.5088837/94‑1 (MQ=60)
                      GACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCC                                                <  SRR4302547.1877855/94‑1 (MQ=60)
                       ACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCG                                               <  SRR4302547.3238942/94‑1 (MQ=60)
                       ACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCG                                               >  SRR4302547.4001781/1‑94 (MQ=60)
                       ACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCG                                               >  SRR4302547.1543196/1‑94 (MQ=60)
                         CGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑‑‑AGCCCCGTCGATGCCATGTTCGAGAACTTTCCCTTCCCGCA                                             <  SRR4302547.5182080/94‑1 (MQ=60)
                          GGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAG                                            >  SRR4302547.2880729/1‑94 (MQ=60)
                          GGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAG                                            <  SRR4302547.731487/94‑1 (MQ=60)
                          GGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAG                                            >  SRR4302547.6024079/1‑94 (MQ=60)
                              GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        <  SRR4302547.1506023/94‑1 (MQ=60)
                              GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        >  SRR4302547.4836435/1‑94 (MQ=60)
                              GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        >  SRR4302547.5123762/1‑94 (MQ=60)
                              GAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACA                                        >  SRR4302547.5975620/1‑94 (MQ=60)
                                AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      <  SRR4302547.2242119/94‑1 (MQ=60)
                                AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      >  SRR4302547.5738099/1‑94 (MQ=60)
                                AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      >  SRR4302547.4104736/1‑94 (MQ=60)
                                 AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     <  SRR4302547.5797044/94‑1 (MQ=60)
                                 AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     <  SRR4302547.2589614/94‑1 (MQ=60)
                                 AACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACG                                     <  SRR4302547.6405131/94‑1 (MQ=60)
                                  ACTACGCCGACAGCCGGCAGGCTAACATCATCGATATCTGGCGC‑‑AACA‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCCCAGAACAAGGG                                    >  SRR4302547.861159/1‑94 (MQ=60)
                                  ACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGG                                    >  SRR4302547.4758855/1‑94 (MQ=60)
                                   CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   >  SRR4302547.2463314/1‑94 (MQ=60)
                                   CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   <  SRR4302547.2072609/94‑1 (MQ=60)
                                   CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   >  SRR4302547.3743155/1‑94 (MQ=60)
                                   CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCATTCCCGCAGAACAACGGA                                   >  SRR4302547.1460949/1‑94 (MQ=60)
                                     ACGCTGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAA                                 <  SRR4302547.6512494/94‑1 (MQ=60)
                                     ACGCCGACAGCCAGCAGGCTAACATCATCGATATCTCGCGCAGCA‑‑‑‑CCGTCGATGGCATGCTCGAGAACTATCCCTTCCCGCAGAACAACGGTAA                                 >  SRR4302547.135919/1‑94 (MQ=60)
                                      CGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAC                                >  SRR4302547.3012249/1‑94 (MQ=60)
                                      CGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAC                                >  SRR4302547.5840378/1‑94 (MQ=60)
                                       GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               >  SRR4302547.2742440/1‑94 (MQ=60)
                                       GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               <  SRR4302547.5098847/94‑1 (MQ=60)
                                       GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               <  SRR4302547.5537513/94‑1 (MQ=60)
                                       GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               >  SRR4302547.3580120/1‑94 (MQ=60)
                                         CGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGT                             >  SRR4302547.1861410/1‑94 (MQ=60)
                                          GACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTC                            <  SRR4302547.5443915/94‑1 (MQ=60)
                                          GACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTC                            >  SRR4302547.4851249/1‑94 (MQ=60)
                                          GACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGGAACCGTC                            >  SRR4302547.5933573/1‑94 (MQ=60)
                                           ACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCA                           <  SRR4302547.5510729/94‑1 (MQ=60)
                                           ACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCA                           >  SRR4302547.4379830/1‑94 (MQ=60)
                                           ACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCA                           <  SRR4302547.1368059/94‑1 (MQ=60)
                                             AGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGC                         <  SRR4302547.3023039/94‑1 (MQ=60)
                                              GCCAGCAGGCTAACATCATCGGTATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCA                        >  SRR4302547.5521214/1‑94 (MQ=60)
                                              GCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCA                        <  SRR4302547.5050155/94‑1 (MQ=60)
                                              GCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCA                        >  SRR4302547.5005164/1‑94 (MQ=60)
                                                CAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG                      <  SRR4302547.3460456/94‑1 (MQ=60)
                                                CAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG                      >  SRR4302547.4461924/1‑94 (MQ=60)
                                                 AGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGT                     >  SRR4302547.5168845/1‑94 (MQ=60)
                                                  GCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTC                    <  SRR4302547.4044169/94‑1 (MQ=60)
                                                   CAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCAATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCC                   >  SRR4302547.674598/1‑94 (MQ=60)
                                                    AGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCG                  <  SRR4302547.920212/94‑1 (MQ=60)
                                                      GCTACCATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCT                <  SRR4302547.3549869/94‑1 (MQ=60)
                                                       CTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTTTCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               >  SRR4302547.2626977/1‑94 (MQ=60)
                                                        TAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGG              >  SRR4302547.2155865/1‑94 (MQ=60)
                                                         AACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGA             <  SRR4302547.3807203/94‑1 (MQ=60)
                                                          ACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGAC            <  SRR4302547.406419/94‑1 (MQ=60)
                                                          ACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGAC            >  SRR4302547.667063/1‑94 (MQ=60)
                                                          ACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGAC            >  SRR4302547.490586/1‑94 (MQ=60)
                                                          ACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAAATATCGCTTCCAGCAGAACAACGGTAACCGTCAGCATGTTCGCTGGAC            >  SRR4302547.6285138/1‑94 (MQ=60)
                                                           CATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACG           >  SRR4302547.2449003/1‑94 (MQ=60)
                                                               ATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAC       >  SRR4302547.4708659/1‑94 (MQ=60)
                                                                TCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACT      <  SRR4302547.2815323/94‑1 (MQ=60)
                                                                 CGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTG     <  SRR4302547.2289567/94‑1 (MQ=60)
                                                                 CGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTG     <  SRR4302547.5386737/94‑1 (MQ=60)
                                                                  GATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGA    <  SRR4302547.1003361/94‑1 (MQ=60)
                                                                  GATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCATTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGA    >  SRR4302547.3856274/1‑94 (MQ=60)
                                                                   ATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC   <  SRR4302547.3243979/94‑1 (MQ=60)
                                                                   ATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC   >  SRR4302547.2583201/1‑94 (MQ=60)
                                                                   ATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC   >  SRR4302547.6265199/1‑94 (MQ=60)
                                                                    TATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  <  SRR4302547.833932/94‑1 (MQ=60)
                                                                    TATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  <  SRR4302547.1930415/94‑1 (MQ=60)
                                                                                |                                                                                     
CAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑CAAG‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222866‑3223027

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