Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 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% 755.8 / NA 204Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (97/107);  total (97/107)
*NC_0009133,222,9450AG100.0% 724.2 / NA 204Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (97/107);  total (97/107)
*NC_0009133,222,9460AC100.0% 718.9 / NA 204Q891S (CAA→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (97/107);  total (97/107)
*NC_0009133,222,9470GA99.5% 727.9 / NA 204A892T (GCC→ACC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (97/106);  minor base C (0/1);  total (97/107)
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.98e-01

Only 100 of 207 total aligned reads displayed.

GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC  >  NC_000913/3222863‑3223026
                                                                                 ||||                                                                               
gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                         >  1:4646550/1‑93 (MQ=255)
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                        >  2:2850655/1‑93 (MQ=255)
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                        >  1:2605970/1‑93 (MQ=255)
  tCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGc                                                                       <  2:334495/93‑1 (MQ=255)
  tCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGc                                                                       >  1:5334173/1‑93 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                     >  1:4974640/1‑93 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                    >  1:308324/1‑93 (MQ=255)
      gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                   <  2:2076901/93‑1 (MQ=255)
      gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                   <  1:6047166/93‑1 (MQ=255)
      gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                   <  1:5449049/93‑1 (MQ=255)
      gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                   <  2:4014012/93‑1 (MQ=255)
       tGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                  >  2:1184744/1‑93 (MQ=255)
       tGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                  <  1:2884394/93‑1 (MQ=255)
        ggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                 <  1:3248001/93‑1 (MQ=255)
        ggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                 <  1:2081854/93‑1 (MQ=255)
        ggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                 >  2:2161698/1‑93 (MQ=255)
        ggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                 >  2:6043272/1‑93 (MQ=255)
         gCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                                                >  1:5798964/1‑93 (MQ=255)
         gCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                                                >  1:3984378/1‑93 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                              <  2:1714320/93‑1 (MQ=255)
            tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                             >  2:3046769/1‑93 (MQ=255)
            tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                             >  2:2470344/1‑93 (MQ=255)
              ttACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                           <  2:272656/93‑1 (MQ=255)
               tACGGTCGTGGACCGGGCGAAACCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                                          <  1:4586406/93‑1 (MQ=255)
                   gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                      >  1:5164859/1‑93 (MQ=255)
                    tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                     >  1:3012627/1‑93 (MQ=255)
                    tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                     <  2:5138608/93‑1 (MQ=255)
                      gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                   <  1:4692295/93‑1 (MQ=255)
                      gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                   >  2:5835709/1‑93 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                  <  2:3692754/93‑1 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                  <  1:5923377/93‑1 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                  >  1:6017354/1‑93 (MQ=255)
                        ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                 <  1:857974/93‑1 (MQ=255)
                        ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                 <  1:110157/93‑1 (MQ=255)
                         gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                <  2:5228577/93‑1 (MQ=255)
                         gACCGGGCGAAAACAACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                <  1:2708866/93‑1 (MQ=255)
                          aCCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTccc                                               >  1:4143366/1‑93 (MQ=255)
                            cGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                                             >  1:1998186/1‑93 (MQ=255)
                            cGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                                             <  2:4792218/93‑1 (MQ=255)
                            cGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGCCGATGCCATGTTCGTGAACTATACCTTCCCGc                                             >  1:3376659/1‑93 (MQ=255)
                             gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                                            <  1:1026372/93‑1 (MQ=255)
                              ggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAg                                           >  1:44040/1‑93 (MQ=255)
                                 gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                        >  2:2266070/1‑93 (MQ=255)
                                  aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                                       <  1:1393049/93‑1 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                      <  1:4324861/93‑1 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                     <  2:3644707/93‑1 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                     <  2:4167183/93‑1 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                     >  1:5404744/1‑93 (MQ=255)
                                     aCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                                    >  2:613898/1‑93 (MQ=255)
                                     aCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                                    <  1:5880157/93‑1 (MQ=255)
                                      cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                   >  1:946770/1‑93 (MQ=255)
                                       tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                  >  1:6115483/1‑93 (MQ=255)
                                       tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                  <  2:4232421/93‑1 (MQ=255)
                                       tACCCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                  <  1:2371839/93‑1 (MQ=255)
                                         cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                <  2:5899822/93‑1 (MQ=255)
                                         cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                <  2:1589720/93‑1 (MQ=255)
                                          cccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                               <  1:1936076/92‑1 (MQ=255)
                                           ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                              >  1:2193138/1‑93 (MQ=255)
                                           ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                              <  1:4218606/93‑1 (MQ=255)
                                            cGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCg                             <  1:2840564/93‑1 (MQ=255)
                                             gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGt                            >  2:3485818/1‑93 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                           <  2:5697665/93‑1 (MQ=255)
                                              aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                           <  1:2487478/93‑1 (MQ=255)
                                                 gccagcGGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCCGAACAACGGTAACCGTCAGc                        >  1:917623/1‑93 (MQ=255)
                                                 gccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGc                        >  1:3539201/1‑93 (MQ=255)
                                                   cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt                      >  2:4252966/1‑93 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                     >  2:2410901/1‑93 (MQ=255)
                                                    agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCCTCCCGCAGAACAACGGTAACCGCCAGCATg                     >  1:3629844/1‑93 (MQ=255)
                                                      cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                   <  2:4611739/93‑1 (MQ=255)
                                                       agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTcc                  >  2:2606649/1‑93 (MQ=255)
                                                        ggttAACAACATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCg                 <  2:2451475/90‑1 (MQ=255)
                                                        ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCCGCGTgcgcg                 >  2:563222/1‑89 (MQ=255)
                                                        ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCg                 >  2:5392878/1‑93 (MQ=255)
                                                        ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCg                 >  2:129570/1‑93 (MQ=255)
                                                        ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCg                 >  2:1738732/1‑93 (MQ=255)
                                                        ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCg                 >  1:1328246/1‑93 (MQ=255)
                                                         gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                >  1:1545391/1‑93 (MQ=255)
                                                         gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                >  2:3299174/1‑93 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg              >  1:2304366/1‑93 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg              >  1:4543457/1‑93 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg              <  2:5439590/93‑1 (MQ=255)
                                                           tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCCTCAGCATGTCCGCTg              >  1:4447527/1‑93 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg             >  1:373191/1‑93 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg             >  2:3897139/1‑93 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg             >  1:4731914/1‑93 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa            <  2:3545493/93‑1 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa            <  2:4863870/93‑1 (MQ=255)
                                                             aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa            <  2:1493644/93‑1 (MQ=255)
                                                              catcatCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGAc           <  1:383291/93‑1 (MQ=255)
                                                              catcatCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGAc           <  1:1760985/93‑1 (MQ=255)
                                                               atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACg          <  1:1105698/93‑1 (MQ=255)
                                                               atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACg          <  2:3409387/93‑1 (MQ=255)
                                                                tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACgg         >  2:484992/1‑93 (MQ=255)
                                                                tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACgg         <  2:1014689/93‑1 (MQ=255)
                                                                ccaccGACATCTGGCGCAGCACCGTCGATCGCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACgg         <  2:399722/89‑1 (MQ=255)
                                                                 catcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGc        >  2:1968348/1‑93 (MQ=255)
                                                                  ctcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCa       <  1:1667941/92‑1 (MQ=255)
                                                                   tcGATATCTGGCGCAGCACCGTCGATGCCATGTACGAGAACTATCCCTTCCCGCAGAACAACGGTAACAGTCAGCATGTCCGCTGGACGGCAc      <  2:5941917/93‑1 (MQ=255)
                                                                    cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACt     <  1:2951586/93‑1 (MQ=255)
                                                                       tatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAc  <  1:2978303/93‑1 (MQ=255)
                                                                                 ||||                                                                               
GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAC  >  NC_000913/3222863‑3223026

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AACGGCGAATACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑CAAG‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222851‑3223027
                                                                                               |                                                                                     
AACGGCGAATACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑A                                                                                       <  SRR4302532.2056587/94‑1 (MQ=60)
             ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGC                                                                        >  SRR4302532.5519167/1‑94 (MQ=60)
             ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGC                                                                        >  SRR4302532.5577145/1‑94 (MQ=60)
             ATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGC                                                                        >  SRR4302532.4557715/1‑94 (MQ=60)
              TCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCC                                                                       >  SRR4302532.39629/1‑94 (MQ=60)
               CAGGTGTCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCA                                                                      <  SRR4302532.1074731/94‑1 (MQ=60)
               CAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCA                                                                      >  SRR4302532.5145163/1‑94 (MQ=60)
                AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGTCGCAGCA‑‑‑‑CCGTCGATGCCAT                                                                     <  SRR4302532.4813828/94‑1 (MQ=60)
                AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCAT                                                                     >  SRR4302532.319011/1‑94 (MQ=60)
                AGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCAT                                                                     >  SRR4302532.2635912/1‑94 (MQ=60)
                 GGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATG                                                                    >  SRR4302532.3001694/1‑94 (MQ=60)
                 GGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATG                                                                    <  SRR4302532.4863421/94‑1 (MQ=60)
                 GGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATG                                                                    >  SRR4302532.5297100/1‑94 (MQ=60)
                  GTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGT                                                                   <  SRR4302532.5637910/94‑1 (MQ=60)
                  GTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGT                                                                   <  SRR4302532.6255976/94‑1 (MQ=60)
                   TGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTT                                                                  <  SRR4302532.2980786/94‑1 (MQ=60)
                    GGCGTATTGCGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTC                                                                 <  SRR4302532.476554/94‑1 (MQ=60)
                    GGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTC                                                                 >  SRR4302532.6000510/1‑94 (MQ=60)
                    GGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTC                                                                 <  SRR4302532.3355837/94‑1 (MQ=60)
                    GGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTC                                                                 >  SRR4302532.4118597/1‑94 (MQ=60)
                    GGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTC                                                                 <  SRR4302532.2151247/94‑1 (MQ=60)
                    GGCGTATTACGGTCGTGGACCGGGCGAAAACCAAGCCGACAGCCCGCAAGCTAACATCATCTATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTC                                                                 >  SRR4302532.1684594/1‑94 (MQ=57)
                      CGTATTACGGTCGTGGACCGGGCGAAACCTACACCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGA                                                               <  SRR4302532.2577644/94‑1 (MQ=60)
                       GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAG                                                              <  SRR4302532.3175027/94‑1 (MQ=60)
                       GTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAG                                                              >  SRR4302532.49662/1‑94 (MQ=60)
                        TATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGA                                                             >  SRR4302532.4166788/1‑94 (MQ=60)
                        TATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGA                                                             >  SRR4302532.5633242/1‑94 (MQ=60)
                        TATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGA                                                             >  SRR4302532.683191/1‑94 (MQ=60)
                         ATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAA                                                            >  SRR4302532.2155053/1‑94 (MQ=60)
                         ATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAA                                                            >  SRR4302532.3666979/1‑94 (MQ=60)
                          TTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAAC                                                           >  SRR4302532.3092144/1‑94 (MQ=60)
                           TACGGTCGTGGACCGGGCGAAACCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACT                                                          <  SRR4302532.4742787/94‑1 (MQ=60)
                           TACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACT                                                          <  SRR4302532.4050434/94‑1 (MQ=60)
                             CGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTAT                                                        <  SRR4302532.1158135/94‑1 (MQ=60)
                              GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       >  SRR4302532.5342660/1‑94 (MQ=60)
                              GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       <  SRR4302532.6195184/94‑1 (MQ=60)
                              GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATC                                                       >  SRR4302532.3068867/1‑94 (MQ=60)
                               GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCC                                                      >  SRR4302532.3112074/1‑94 (MQ=60)
                                tcggtgacaGGGCGACAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGCCGATGCCATGTTCGAGAACTATCCC                                                     <  SRR4302532.31822/85‑1 (MQ=60)
                                TCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCC                                                     >  SRR4302532.555657/1‑94 (MQ=60)
                                TCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCC                                                     <  SRR4302532.6026715/94‑1 (MQ=60)
                                 CGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCT                                                    <  SRR4302532.6294885/94‑1 (MQ=60)
                                  GTGGACCGGGCGAAAACTACGCCTACAGGCAGCAGGCTAACATCATCCATTTCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   <  SRR4302532.4421208/94‑1 (MQ=60)
                                  GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   <  SRR4302532.963129/94‑1 (MQ=60)
                                  GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   <  SRR4302532.1378730/94‑1 (MQ=60)
                                  GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   <  SRR4302532.4852332/94‑1 (MQ=60)
                                  GTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTT                                                   >  SRR4302532.6225149/1‑94 (MQ=60)
                                   TGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTC                                                  >  SRR4302532.5434428/1‑94 (MQ=60)
                                   TGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTC                                                  <  SRR4302532.6128486/94‑1 (MQ=60)
                                    GGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCC                                                 <  SRR4302532.113961/94‑1 (MQ=60)
                                    GGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCC                                                 <  SRR4302532.113962/94‑1 (MQ=60)
                                    GGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCC                                                 <  SRR4302532.4734130/94‑1 (MQ=60)
                                    GGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCC                                                 <  SRR4302532.887454/94‑1 (MQ=60)
                                     GACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCC                                                >  SRR4302532.4282368/1‑94 (MQ=60)
                                     GACCGGGCGAAAACAACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCC                                                <  SRR4302532.2799963/94‑1 (MQ=60)
                                       CCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGC                                              >  SRR4302532.2064762/1‑94 (MQ=60)
                                       CCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGCCGATGCCATGTTCGTGAACTATACCTTCCCGC                                              >  SRR4302532.3489584/1‑94 (MQ=60)
                                         GGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAG                                            >  SRR4302532.45602/1‑94 (MQ=60)
                                         GGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAG                                            <  SRR4302532.1061097/94‑1 (MQ=60)
                                            CGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAAC                                         >  SRR4302532.3236019/1‑94 (MQ=60)
                                            CGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAAC                                         <  SRR4302532.2084622/94‑1 (MQ=60)
                                              AAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAA                                       <  SRR4302532.1439756/94‑1 (MQ=60)
                                               AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      >  SRR4302532.5592002/1‑94 (MQ=60)
                                               AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      <  SRR4302532.2762234/94‑1 (MQ=60)
                                               AAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAAC                                      <  SRR4302532.4471000/94‑1 (MQ=60)
                                                 ACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGG                                    <  SRR4302532.6083845/94‑1 (MQ=60)
                                                 ACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGG                                    >  SRR4302532.979167/1‑94 (MQ=60)
                                                  CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   >  SRR4302532.6326535/1‑94 (MQ=60)
                                                  CTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGT                                   <  SRR4302532.2854711/94‑1 (MQ=60)
                                                   TACCCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTA                                  <  SRR4302532.2450910/94‑1 (MQ=60)
                                                      GCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               >  SRR4302532.2266209/1‑94 (MQ=60)
                                                      CCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACC                               <  SRR4302532.2000504/94‑1 (MQ=60)
                                                       CCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCG                              >  SRR4302532.4421640/1‑94 (MQ=60)
                                                       CCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCG                              >  SRR4302532.2431676/1‑94 (MQ=60)
                                                       CCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCG                              <  SRR4302532.4360201/94‑1 (MQ=60)
                                                        CGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGT                             <  SRR4302532.2935508/94‑1 (MQ=60)
                                                          ACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCA                           <  SRR4302532.1753812/94‑1 (MQ=60)
                                                          ACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCA                           <  SRR4302532.2570358/94‑1 (MQ=60)
                                                           CAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAG                          >  SRR4302532.4537244/1‑94 (MQ=60)
                                                            AGCCAGCGGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCCGAACAACGGTAACCGTCAGC                         >  SRR4302532.949021/1‑94 (MQ=60)
                                                            AGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGC                         <  SRR4302532.2459178/94‑1 (MQ=60)
                                                            AGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGC                         >  SRR4302532.3657605/1‑94 (MQ=60)
                                                            AGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGC                         >  SRR4302532.3484969/1‑94 (MQ=60)
                                                              CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       >  SRR4302532.5973742/1‑94 (MQ=60)
                                                              CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       <  SRR4302532.3548661/94‑1 (MQ=60)
                                                              CCAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAT                       >  SRR4302532.1930725/1‑94 (MQ=60)
                                                               CAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG                      >  SRR4302532.2961880/1‑94 (MQ=60)
                                                               CAGCAGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCCTCCCGCAGAACAACGGTAACCGCCAGCATG                      >  SRR4302532.3751596/1‑94 (MQ=60)
                                                                   AGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCG                  <  SRR4302532.2319687/94‑1 (MQ=60)
                                                                   AGGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCG                  >  SRR4302532.1372637/1‑94 (MQ=60)
                                                                   AGGCTAACATCATAGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGACAACTCTCCCTTCCCGCAGAACAACGGCAACCGTCAGGATGTCCG                  >  SRR4302532.5589995/1‑94 (MQ=57)
                                                                    GGCTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGC                 >  SRR4302532.1597611/1‑94 (MQ=60)
                                                                      CTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               >  SRR4302532.4698544/1‑94 (MQ=60)
                                                                      CTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTG               >  SRR4302532.2380979/1‑94 (MQ=60)
                                                                      CTAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCCTCAGCATGTCCGCTG               >  SRR4302532.4599105/1‑94 (MQ=60)
                                                                       TAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGG              >  SRR4302532.386145/1‑94 (MQ=60)
                                                                       TAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGG              >  SRR4302532.3760745/1‑94 (MQ=60)
                                                                       TAACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGG              >  SRR4302532.4893424/1‑94 (MQ=60)
                                                                         ACATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGAC            <  SRR4302532.6088938/94‑1 (MQ=60)
                                                                          CATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACG           <  SRR4302532.396683/94‑1 (MQ=60)
                                                                          CATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACG           <  SRR4302532.1820060/94‑1 (MQ=60)
                                                                          CATCATCGATATCTGGCGC‑‑‑‑AGCCCCGTCGATGCCATGTTCTAGAACTATCCCTTCCTGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACG           <  SRR4302532.6226311/94‑1 (MQ=60)
                                                                           ATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGcggg          >  SRR4302532.1607385/1‑90 (MQ=60)
                                                                           ATCATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGG          <  SRR4302532.1142662/94‑1 (MQ=60)
                                                                             CATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCA        >  SRR4302532.4333208/1‑94 (MQ=60)
                                                                              CTCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAC       <  SRR4302532.1724389/94‑1 (MQ=60)
                                                                              ATCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAC       >  SRR4302532.4667136/1‑94 (MQ=60)
                                                                               TCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACT      <  SRR4302532.1571229/94‑1 (MQ=60)
                                                                               TCGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACT      <  SRR4302532.944224/94‑1 (MQ=60)
                                                                                CGATATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTG     <  SRR4302532.3049434/94‑1 (MQ=60)
                                                                                   TATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  <  SRR4302532.3076840/94‑1 (MQ=60)
                                                                                   TATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  SRR4302532.1772637/1‑94 (MQ=60)
                                                                                   TATCTGGCGCAGCA‑‑‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  SRR4302532.3878602/1‑94 (MQ=60)
                                                                                               |                                                                                     
AACGGCGAATACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGC‑‑CAAG‑‑CCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222851‑3223027

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