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 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,108,133 A→G G208G (GGA→GGG cdh → CDP‑diacylglycerol phosphotidylhydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,108,1330AG99.5% 671.1 / NA 209G208G (GGA→GGGcdhCDP‑diacylglycerol phosphotidylhydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (98/110);  minor base T (0/1);  total (98/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

GCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGACGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTC  >  NC_000913/4108044‑4108222
                                                                                         |                                                                                         
   cgTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTAc                                                                                     <  1:2745262/93‑1 (MQ=255)
   cgTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTAc                                                                                     >  2:1025951/1‑93 (MQ=255)
    gTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACg                                                                                    <  1:175719/93‑1 (MQ=255)
    gTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACg                                                                                    >  2:1929473/1‑93 (MQ=255)
    gTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACg                                                                                    <  1:2973832/93‑1 (MQ=255)
     tCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACgg                                                                                   <  1:3762646/93‑1 (MQ=255)
     tCGGGGAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTActg                                                                                   >  2:4489510/1‑91 (MQ=255)
       gggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGc                                                                                 <  2:4173373/93‑1 (MQ=255)
        ggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCt                                                                                <  2:4658394/93‑1 (MQ=255)
        ggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCt                                                                                <  2:6272694/93‑1 (MQ=255)
        ggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCt                                                                                >  1:20173/1‑93 (MQ=255)
        ggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGAAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCt                                                                                <  2:5919486/93‑1 (MQ=255)
         gTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTg                                                                               >  2:6256095/1‑93 (MQ=255)
         gTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTg                                                                               >  2:4796371/1‑93 (MQ=255)
         gTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTg                                                                               >  1:5932927/1‑93 (MQ=255)
          tAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTgg                                                                              <  1:1634322/93‑1 (MQ=255)
           aaCGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGc                                                                             <  1:1056885/93‑1 (MQ=255)
            aCGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCg                                                                            >  2:1302662/1‑93 (MQ=255)
             cGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGa                                                                           <  1:3966645/93‑1 (MQ=255)
             cGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGa                                                                           <  1:1220089/93‑1 (MQ=255)
               gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGCACATATGGGGCGCTACGGGCTGGCGGTg                                                                         >  2:1250042/1‑93 (MQ=255)
               gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATg                                                                         <  1:4756132/93‑1 (MQ=255)
               gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATg                                                                         >  2:5452283/1‑93 (MQ=255)
               gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATg                                                                         <  2:276778/93‑1 (MQ=255)
               gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATg                                                                         >  1:4995429/1‑93 (MQ=255)
                aaaGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATgg                                                                        <  1:4272552/93‑1 (MQ=255)
                   gCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTAc                                                                     <  1:2820186/93‑1 (MQ=255)
                   gCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTAc                                                                     <  2:1513572/93‑1 (MQ=255)
                     gAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACgg                                                                   <  1:1147937/93‑1 (MQ=255)
                     gAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACgg                                                                   >  2:3995782/1‑93 (MQ=255)
                       aCTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCa                                                                 <  2:5029680/93‑1 (MQ=255)
                       aCTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCag                                                                >  2:5355991/1‑93 (MQ=255)
                        cTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCag                                                                >  1:6035967/1‑93 (MQ=255)
                        cTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCag                                                                <  1:4732808/93‑1 (MQ=255)
                        cTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCag                                                                >  2:1627659/1‑93 (MQ=255)
                         tGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCaga                                                               <  2:2136078/93‑1 (MQ=255)
                         tGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCaga                                                               <  2:2553872/93‑1 (MQ=255)
                         tGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCaga                                                               >  1:55384/1‑93 (MQ=255)
                           gTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGt                                                             <  2:5320314/93‑1 (MQ=255)
                           gTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGt                                                             <  1:523727/93‑1 (MQ=255)
                             aCAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGa                                                           <  2:4078011/93‑1 (MQ=255)
                             aCAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGa                                                           >  1:4530795/1‑93 (MQ=255)
                             aCAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGa                                                           >  2:5693404/1‑93 (MQ=255)
                              cAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGAt                                                          <  2:2674345/93‑1 (MQ=255)
                               aaCGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATa                                                         <  1:223799/93‑1 (MQ=255)
                                aCGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACGTATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATaa                                                        <  1:4436060/93‑1 (MQ=255)
                                 cGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAAc                                                       >  1:5648042/1‑93 (MQ=255)
                                     gCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAt                                                   >  1:5461250/1‑93 (MQ=255)
                                     gCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAt                                                   >  2:4402999/1‑93 (MQ=255)
                                     gCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAt                                                   >  1:4823859/1‑93 (MQ=255)
                                       ccATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAttt                                                 <  1:1822118/93‑1 (MQ=255)
                                       ccATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAttt                                                 <  2:5852662/93‑1 (MQ=255)
                                       ccATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAttt                                                 <  1:4889167/93‑1 (MQ=255)
                                       ccATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAttt                                                 <  2:3943121/93‑1 (MQ=255)
                                        cATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTg                                                >  2:252052/1‑93 (MQ=255)
                                          tttATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTa                                              >  1:3281559/1‑93 (MQ=255)
                                          tttATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTa                                              >  2:3238430/1‑93 (MQ=255)
                                           ttATGATGCTGGCAGAAGAAGTACCTGAGTCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAt                                             <  1:1441946/93‑1 (MQ=255)
                                           ttATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAt                                             <  2:4239461/93‑1 (MQ=255)
                                              tgatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGc                                          <  1:5116618/93‑1 (MQ=255)
                                              tgatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGc                                          <  1:2250379/93‑1 (MQ=255)
                                               gatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCt                                         <  1:2940682/93‑1 (MQ=255)
                                                atgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTg                                        <  1:2544278/93‑1 (MQ=255)
                                                atgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTg                                        <  1:4004143/93‑1 (MQ=255)
                                                 tgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTgg                                       <  1:4323884/93‑1 (MQ=255)
                                                 tgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTgg                                       >  1:1358698/1‑93 (MQ=255)
                                                  gCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGc                                      >  2:4324943/1‑93 (MQ=255)
                                                  gCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGc                                      <  2:816521/93‑1 (MQ=255)
                                                    tGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGTTGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGa                                    <  1:580571/93‑1 (MQ=255)
                                                        agaagaAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGAcaca                                <  2:3786832/93‑1 (MQ=255)
                                                         ggagaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACaa                               >  2:3296860/3‑93 (MQ=255)
                                                         gaagaaGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTTATTTGTATTGCTGGCGACACaa                               >  1:1711265/1‑93 (MQ=255)
                                                         gaagaaGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACaa                               >  1:681795/1‑93 (MQ=255)
                                                          aagaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAAc                              <  2:1839530/93‑1 (MQ=255)
                                                          aagaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAAc                              <  2:2312776/93‑1 (MQ=255)
                                                              agTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGaaa                          <  1:3583733/93‑1 (MQ=255)
                                                              agTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGaaa                          >  1:2160873/1‑93 (MQ=255)
                                                              agTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGaaa                          >  2:1340752/1‑93 (MQ=255)
                                                              agTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGaaa                          <  2:1206824/93‑1 (MQ=255)
                                                               gTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAAt                         >  1:4039027/1‑93 (MQ=255)
                                                                tACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGGGGTGGTACGGCAGAGTGAGAAATCACTTGTATTTCTTGGTACACAACGaaaaa                        >  2:3143579/1‑91 (MQ=255)
                                                                tACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATc                        >  2:6217608/1‑93 (MQ=255)
                                                                  ccTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTa                      <  1:2787489/93‑1 (MQ=255)
                                                                   cTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTAc                     <  2:705815/93‑1 (MQ=255)
                                                                   cTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTAc                     >  1:2336550/1‑93 (MQ=255)
                                                                    tGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACt                    >  2:4220173/1‑93 (MQ=255)
                                                                    tGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACt                    <  1:1576073/93‑1 (MQ=255)
                                                                    tGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACt                    <  2:3127029/93‑1 (MQ=255)
                                                                      aGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGa                  <  1:5990646/93‑1 (MQ=255)
                                                                            gggAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTa            <  2:4781727/93‑1 (MQ=255)
                                                                              gAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAAc          >  2:3883594/1‑93 (MQ=255)
                                                                               aaCATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAAcc         <  1:5720278/93‑1 (MQ=255)
                                                                                  atatGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTg      >  1:3712996/1‑93 (MQ=255)
                                                                                  atatGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTg      >  1:703651/1‑93 (MQ=255)
                                                                                   tatGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGc     >  2:1764788/1‑93 (MQ=255)
                                                                                   tatGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGc     >  1:1308661/1‑93 (MQ=255)
                                                                                     tGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTCACGCTTAACCGTGCtt   <  1:3319808/93‑1 (MQ=255)
                                                                                      ggggCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTc  <  1:2077959/93‑1 (MQ=255)
                                                                                      ggggCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTc  <  2:4880930/93‑1 (MQ=255)
                                                                                      ggggCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTc  <  1:982281/93‑1 (MQ=255)
                                                                                         |                                                                                         
GCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGACGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTC  >  NC_000913/4108044‑4108222

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ 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
CCTGGCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGACGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCC  >  NC_000913/4108040‑4108226
                                                                                             |                                                                                             
CCTGGCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGG                                                                                               >  SRR4302533.6063634/1‑94 (MQ=60)
CCTGGCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGG                                                                                               <  SRR4302533.3857597/94‑1 (MQ=60)
CCTGGCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGG                                                                                               >  SRR4302533.5300951/1‑94 (MQ=60)
CCTGGCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGG                                                                                               <  SRR4302533.554096/94‑1 (MQ=60)
   GGCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGC                                                                                            >  SRR4302533.5263788/1‑94 (MQ=60)
       CGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGGACATATGGGGCGCTACG                                                                                        >  SRR4302533.5791873/1‑94 (MQ=60)
       CGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACG                                                                                        <  SRR4302533.2762091/94‑1 (MQ=60)
       CGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACG                                                                                        >  SRR4302533.4596383/1‑94 (MQ=60)
        GTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGG                                                                                       <  SRR4302533.176825/94‑1 (MQ=60)
        GTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGG                                                                                       <  SRR4302533.2992029/94‑1 (MQ=60)
        GTCGGGTAACGGAAACCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGG                                                                                       <  SRR4302533.3646372/94‑1 (MQ=60)
        GTCGGGCAACGGAAAGCGAACTGGTACAAAGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGG                                                                                       >  SRR4302533.5746576/1‑94 (MQ=60)
         TCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGG                                                                                      <  SRR4302533.3785711/94‑1 (MQ=60)
           GGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCT                                                                                    >  SRR4302533.22110/1‑94 (MQ=60)
           GGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCT                                                                                    <  SRR4302533.2042119/94‑1 (MQ=60)
           GGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCT                                                                                    >  SRR4302533.20310/1‑94 (MQ=60)
            GGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTG                                                                                   <  SRR4302533.3959394/94‑1 (MQ=60)
            GGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTG                                                                                   <  SRR4302533.3259711/94‑1 (MQ=60)
            GGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTG                                                                                   >  SRR4302533.5968794/1‑94 (MQ=60)
             GTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGG                                                                                  <  SRR4302533.1313126/94‑1 (MQ=60)
              TAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGC                                                                                 <  SRR4302533.1644221/94‑1 (MQ=60)
              TAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGC                                                                                 <  SRR4302533.4864679/94‑1 (MQ=60)
              TAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGC                                                                                 <  SRR4302533.193735/94‑1 (MQ=60)
               AACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCG                                                                                <  SRR4302533.1063100/94‑1 (MQ=60)
                 CGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGAT                                                                              <  SRR4302533.3990808/94‑1 (MQ=60)
                 CGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGAT                                                                              <  SRR4302533.1227412/94‑1 (MQ=60)
                 CGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGAG                                                                              >  SRR4302533.6092954/1‑94 (MQ=60)
                  GGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATG                                                                             >  SRR4302533.5025713/1‑94 (MQ=60)
                   GAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGG                                                                            <  SRR4302533.4784995/94‑1 (MQ=60)
                   GAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGG                                                                            >  SRR4302533.1401976/1‑94 (MQ=60)
                   GAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGG                                                                            <  SRR4302533.1456767/94‑1 (MQ=60)
                   GAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGG                                                                            >  SRR4302533.5555254/1‑94 (MQ=60)
                    AAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGT                                                                           <  SRR4302533.4298265/94‑1 (MQ=60)
                    AAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGT                                                                           >  SRR4302533.1001919/1‑94 (MQ=60)
                     AAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTA                                                                          <  SRR4302533.4823413/94‑1 (MQ=60)
                      AGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTAC                                                                         <  SRR4302533.3161104/94‑1 (MQ=60)
                       GCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACG                                                                        <  SRR4302533.1692303/94‑1 (MQ=60)
                       GCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACG                                                                        <  SRR4302533.2837542/94‑1 (MQ=60)
                         GAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGC                                                                      <  SRR4302533.3315887/94‑1 (MQ=60)
                         GAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGC                                                                      <  SRR4302533.1154791/94‑1 (MQ=60)
                          AACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCA                                                                     >  SRR4302533.1068110/1‑94 (MQ=60)
                           ACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAG                                                                    >  SRR4302533.6072545/1‑94 (MQ=60)
                            CTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGA                                                                   <  SRR4302533.4761499/94‑1 (MQ=60)
                            CTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGA                                                                   >  SRR4302533.55756/1‑94 (MQ=60)
                             TGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGTCGCTACGGGCTGGCGATGGTACGGCAGAG                                                                  <  SRR4302533.2601608/94‑1 (MQ=60)
                              GGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGG                                                                 >  SRR4302533.6031941/1‑94 (MQ=60)
                               GTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTG                                                                <  SRR4302533.526835/94‑1 (MQ=60)
                                TACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGA                                                               >  SRR4302533.4558230/1‑94 (MQ=60)
                                  CAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATA                                                             <  SRR4302533.418019/94‑1 (MQ=60)
                                  CAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATA                                                             >  SRR4302533.3879531/1‑94 (MQ=60)
                                   AACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGCTGGTACGGCAGAGTGATAA                                                            >  SRR4302533.5889799/1‑94 (MQ=60)
                                   AACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAA                                                            >  SRR4302533.4051856/1‑94 (MQ=60)
                                   AACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAA                                                            <  SRR4302533.225196/94‑1 (MQ=60)
                                   AACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAA                                                            >  SRR4302533.2063673/1‑94 (MQ=60)
                                    ACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACGTATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAAC                                                           <  SRR4302533.4462877/94‑1 (MQ=60)
                                    ACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAAC                                                           >  SRR4302533.5682185/1‑94 (MQ=60)
                                      GCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTC                                                         <  SRR4302533.1479120/94‑1 (MQ=60)
                                      GCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGAGCTACGGGCTGGCGCTGGTACGGCAGAGTGATAACTC                                                         >  SRR4302533.2934775/1‑94 (MQ=60)
                                        AGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAT                                                       >  SRR4302533.5494152/1‑94 (MQ=60)
                                        AGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAT                                                       >  SRR4302533.4853080/1‑94 (MQ=60)
                                        AGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAT                                                       >  SRR4302533.564534/1‑94 (MQ=60)
                                          CCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTT                                                     <  SRR4302533.1319097/94‑1 (MQ=60)
                                          CCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTT                                                     <  SRR4302533.3569994/94‑1 (MQ=60)
                                           CCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTG                                                    <  SRR4302533.4918793/94‑1 (MQ=60)
                                           CCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTG                                                    <  SRR4302533.1833249/94‑1 (MQ=60)
                                           CCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTG                                                    <  SRR4302533.1738721/94‑1 (MQ=60)
                                             ATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTA                                                  >  SRR4302533.277414/1‑94 (MQ=60)
                                             ATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTA                                                  >  SRR4302533.3301819/1‑94 (MQ=60)
                                              TTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAT                                                 >  SRR4302533.1832612/1‑94 (MQ=60)
                                              TTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAT                                                 <  SRR4302533.5086923/94‑1 (MQ=60)
                                              TTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAT                                                 <  SRR4302533.377747/94‑1 (MQ=60)
                                               TTATGATGCTGGCAGAAGAAGTACCTGAGTCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATT                                                <  SRR4302533.1450646/94‑1 (MQ=60)
                                                 ATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGC                                              <  SRR4302533.2458987/94‑1 (MQ=60)
                                                 ATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGC                                              >  SRR4302533.4522980/1‑94 (MQ=60)
                                                 ATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGC                                              >  SRR4302533.692141/1‑94 (MQ=60)
                                                  TGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCT                                             <  SRR4302533.2264288/94‑1 (MQ=60)
                                                  TGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCT                                             <  SRR4302533.5147539/94‑1 (MQ=60)
                                                   GATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTG                                            <  SRR4302533.2958643/94‑1 (MQ=60)
                                                    ATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGG                                           <  SRR4302533.4028521/94‑1 (MQ=60)
                                                    ATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGG                                           <  SRR4302533.4967825/94‑1 (MQ=60)
                                                    ATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGG                                           >  SRR4302533.1366865/1‑94 (MQ=60)
                                                    ATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGG                                           <  SRR4302533.2559935/94‑1 (MQ=60)
                                                    ATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGG                                           >  SRR4302533.2833412/1‑94 (MQ=60)
                                                    ATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGG                                           >  SRR4302533.1158953/1‑94 (MQ=60)
                                                     TGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGC                                          <  SRR4302533.4349921/94‑1 (MQ=60)
                                                       cgggaAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGA                                        <  SRR4302533.5279684/89‑1 (MQ=60)
                                                        TGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGTTGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGAC                                       <  SRR4302533.583987/94‑1 (MQ=60)
                                                         GGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACA                                      <  SRR4302533.2436915/94‑1 (MQ=60)
                                                            AGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTTATTTGTATTGCTGGCGACACAA                                   >  SRR4302533.1721634/1‑94 (MQ=60)
                                                            AGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAA                                   >  SRR4302533.685827/1‑94 (MQ=60)
                                                             GAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAAC                                  <  SRR4302533.5994267/94‑1 (MQ=60)
                                                               AGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACTCAACGA                                >  SRR4302533.2253662/1‑94 (MQ=60)
                                                                 AAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAA                              >  SRR4302533.2174173/1‑94 (MQ=60)
                                                                  AGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAAT                             <  SRR4302533.2272536/94‑1 (MQ=60)
                                                                  AGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAAT                             >  SRR4302533.4063604/1‑94 (MQ=60)
                                                                  AGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAAT                             <  SRR4302533.3605805/94‑1 (MQ=60)
                                                                  AGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGATGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAAT                             >  SRR4302533.2152443/1‑94 (MQ=60)
                                                                   GTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATC                            <  SRR4302533.4018011/94‑1 (MQ=60)
                                                                   GTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATC                            <  SRR4302533.499741/94‑1 (MQ=60)
                                                                    TACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCT                           <  SRR4302533.5684879/94‑1 (MQ=60)
                                                                      CCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTAC                         >  SRR4302533.2351012/1‑94 (MQ=60)
                                                                      CCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTAC                         <  SRR4302533.2804617/94‑1 (MQ=60)
                                                                       CTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACT                        <  SRR4302533.4326999/94‑1 (MQ=60)
                                                                       CTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACT                        >  SRR4302533.6206744/1‑94 (MQ=60)
                                                                        TGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTG                       >  SRR4302533.4508915/1‑94 (MQ=60)
                                                                        TGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTG                       <  SRR4302533.1585561/94‑1 (MQ=60)
                                                                          AGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGAC                     <  SRR4302533.6026933/94‑1 (MQ=60)
                                                                            GCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGC                   <  SRR4302533.2973064/94‑1 (MQ=60)
                                                                               CGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTA                >  SRR4302533.4769213/1‑94 (MQ=60)
                                                                                GGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACCCAACGAAATCTACTGACGCTTAA               >  SRR4302533.4482313/1‑94 (MQ=60)
                                                                                 GGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAAC              >  SRR4302533.959794/1‑94 (MQ=60)
                                                                                  GAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACC             <  SRR4302533.5324913/94‑1 (MQ=60)
                                                                                   AACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCG            <  SRR4302533.5754816/94‑1 (MQ=60)
                                                                                     CATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTG          >  SRR4302533.707811/1‑94 (MQ=60)
                                                                                     CATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTG          >  SRR4302533.3735796/1‑94 (MQ=60)
                                                                                     CATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTG          >  SRR4302533.4397631/1‑94 (MQ=60)
                                                                                      ATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGC         >  SRR4302533.1316517/1‑94 (MQ=60)
                                                                                      ATATGGGGCGCTACGGGCCGGCGATGGTACGGCAGAGAGATCACTCCTTTGTATTGCTGGCCACACAACGAAATCTACTGACGCTTAACCGTGA         >  SRR4302533.2111824/1‑94 (MQ=60)
                                                                                        ATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTT       >  SRR4302533.1612867/1‑94 (MQ=60)
                                                                                         TGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTCACGCTTAACCGTGCTTC      <  SRR4302533.3340314/94‑1 (MQ=60)
                                                                                          GGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCA     <  SRR4302533.988037/94‑1 (MQ=60)
                                                                                          GGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCA     <  SRR4302533.2090749/94‑1 (MQ=60)
                                                                                          GGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCA     >  SRR4302533.3176780/1‑94 (MQ=60)
                                                                                          GGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCA     >  SRR4302533.4649563/1‑94 (MQ=60)
                                                                                             GCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCC  <  SRR4302533.3532974/94‑1 (MQ=60)
                                                                                             |                                                                                             
CCTGGCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGACGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCC  >  NC_000913/4108040‑4108226

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: