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 882 254.4 12957588 98.0% 12698436 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,393 C→G A159G (GCT→GGT)  b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3930CG100.0% 244.3 / NA 71A159G (GCT→GGT) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (24/47);  total (24/47)

CGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑C‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008310‑4008457
                                                                                      |                                                                         
cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                <  1:1280000/93‑1 (MQ=255)
cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                <  2:916174/93‑1 (MQ=255)
cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                <  2:602061/93‑1 (MQ=255)
cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                <  2:491565/93‑1 (MQ=255)
cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                <  2:3722069/93‑1 (MQ=255)
  aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCC‑AAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                              >  1:3484754/1‑91 (MQ=255)
  aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                              <  1:5387331/93‑3 (MQ=255)
  aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                              <  2:4351240/93‑3 (MQ=255)
  aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                              <  1:4787452/93‑3 (MQ=255)
  aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                              <  2:6281707/93‑3 (MQ=255)
   tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                             <  1:1513701/93‑4 (MQ=255)
   tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                             >  2:6443111/1‑90 (MQ=255)
    cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                            <  2:3015233/93‑5 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                <  1:4486244/93‑1 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                >  2:5403718/1‑93 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                <  2:3787930/93‑1 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                <  2:2774904/93‑1 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                >  2:869479/1‑93 (MQ=255)
                                      tGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                              <  2:464874/93‑3 (MQ=255)
                                      tGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                              >  2:6138832/1‑91 (MQ=255)
                                      tGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                              >  1:6016825/1‑91 (MQ=255)
                                      tGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                              >  1:6016849/1‑91 (MQ=255)
                                       ggAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                             <  1:3506191/93‑4 (MQ=255)
                                       ggAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                             >  2:3320361/1‑90 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  2:3483846/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  2:4040145/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  2:1252666/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggg                            >  2:2507347/1‑89 (MQ=255)
                                              aTGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                      <  1:4940357/93‑1 (MQ=255)
                                               tGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                     <  1:6405431/93‑1 (MQ=255)
                                               tGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                     <  2:1394688/93‑1 (MQ=255)
                                               tGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                     >  1:2083031/1‑93 (MQ=255)
                                               tGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                     <  2:5796766/93‑1 (MQ=255)
                                                ggcgggcgGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                    <  2:3310173/93‑1 (MQ=255)
                                                ggcgggcgGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                    <  1:303919/93‑1 (MQ=255)
                                                  cgggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGATGGGATGa                  >  1:4832777/1‑93 (MQ=255)
                                                  cgggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                  <  2:5573378/93‑1 (MQ=255)
                                                   gggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                 <  2:3446146/93‑1 (MQ=255)
                                                   gggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                 <  2:3417995/93‑1 (MQ=255)
                                                   gggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                 <  1:4711400/93‑1 (MQ=255)
                                                   gggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                 <  2:3173286/93‑1 (MQ=255)
                                                    ggcggcGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                <  1:4031251/93‑1 (MQ=255)
                                                    ggcggcGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                <  1:3439681/93‑1 (MQ=255)
                                                    ggcggcGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                >  1:6362609/1‑93 (MQ=255)
                                                      cggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga              <  1:4625671/93‑1 (MQ=255)
                                                      cggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga              <  1:1316320/93‑1 (MQ=255)
                                                       ggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac             <  2:763229/93‑1 (MQ=255)
                                                       ggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac             <  1:1552969/93‑1 (MQ=255)
                                                       ggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac             >  2:1044481/1‑93 (MQ=255)
                                                        gcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg            >  1:6033668/1‑93 (MQ=255)
                                                        gcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg            >  2:5946770/1‑93 (MQ=255)
                                                        gcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg            >  1:1426432/1‑93 (MQ=255)
                                                        gcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg            <  2:4952602/93‑1 (MQ=255)
                                                        gcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg            >  2:4992638/1‑93 (MQ=255)
                                                         cgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc           >  1:1552563/1‑93 (MQ=255)
                                                         cgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc           <  1:426660/93‑1 (MQ=255)
                                                         cgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc           <  2:4958176/93‑1 (MQ=255)
                                                         cgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc           <  2:5168728/93‑1 (MQ=255)
                                                          gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          >  2:2256472/1‑93 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGa        <  1:6062415/89‑1 (MQ=255)
                                                              tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa       >  2:5083323/4‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaag     <  2:2846134/92‑2 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt     >  1:6035252/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt     >  2:4617817/2‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg    <  1:5923994/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc   <  1:3789541/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc   >  2:5226036/1‑93 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc   <  2:1736980/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  <  1:6098058/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  >  1:82540/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  >  1:2452004/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  <  1:1120950/93‑1 (MQ=255)
                                                                                      |                                                                         
CGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑C‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008310‑4008457

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
GGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGT  >  NC_000913/4008308‑4008459
                                                                                       |                                                                           
GGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑A                                                                   <  SRR4302537.2532390/94‑1 (MQ=60)
                                                                caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT       >  SRR4302537.6076720/3‑94 (MQ=18)
                                                                  AAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGC     <  SRR4302537.5964645/94‑1 (MQ=21)
                                                                   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT    >  SRR4302537.2469257/1‑94 (MQ=21)
                                                                   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT    <  SRR4302537.3816140/94‑1 (MQ=21)
                                                                   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT    >  SRR4302537.83144/1‑94 (MQ=21)
                                                                    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG   <  SRR4302537.1128719/94‑1 (MQ=23)
                                                                    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG   >  SRR4302537.3086211/1‑94 (MQ=23)
                                                                    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG   <  SRR4302537.6139961/94‑1 (MQ=23)
                                                                     GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT  >  SRR4302537.6155501/1‑94 (MQ=25)
                                                                                       |                                                                           
GGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGT  >  NC_000913/4008308‑4008459

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑