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,411 +C coding (494/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4101.C100.0% 257.5 / NA 74A165P (GCT→CCT) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (27/47);  total (27/47)

TCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGC  >  NC_000913/4008346‑4008492
                                                                        |                                                                                           
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                                                       >  2:5403718/1‑93 (MQ=255)
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                                                       <  2:2774904/93‑1 (MQ=255)
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                                                       <  1:4486244/93‑1 (MQ=255)
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                                                       >  2:869479/1‑93 (MQ=255)
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                                                       <  2:3787930/93‑1 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                                                                     <  2:464874/93‑3 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                                                                     >  2:6138832/1‑91 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                                                                     >  1:6016849/1‑91 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                                                                     >  1:6016825/1‑91 (MQ=255)
   ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttgg                                                                    <  1:3506191/93‑4 (MQ=255)
   ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttgg                                                                    >  2:3320361/1‑90 (MQ=255)
    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttggt                                                                   <  2:4040145/93‑5 (MQ=255)
    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttggt                                                                   <  2:1252666/93‑5 (MQ=255)
    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttggt                                                                   <  2:3483846/93‑5 (MQ=255)
    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttggg                                                                   >  2:2507347/1‑89 (MQ=255)
          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACggg                                                             <  1:4940357/93‑1 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                                                            >  1:2083031/1‑93 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                                                            <  1:6405431/93‑1 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                                                            <  2:1394688/93‑1 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                                                            <  2:5796766/93‑1 (MQ=255)
            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt                                                           <  2:3310173/93‑1 (MQ=255)
            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt                                                           <  1:303919/93‑1 (MQ=255)
              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGATGGGATGa                                                         >  1:4832777/1‑93 (MQ=255)
              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGa                                                         <  2:5573378/93‑1 (MQ=255)
               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                                                        <  1:4711400/93‑1 (MQ=255)
               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                                                        <  2:3173286/93‑1 (MQ=255)
               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                                                        <  2:3417995/93‑1 (MQ=255)
               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                                                        <  2:3446146/93‑1 (MQ=255)
                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                                                       <  1:4031251/93‑1 (MQ=255)
                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                                                       <  1:3439681/93‑1 (MQ=255)
                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                                                       >  1:6362609/1‑93 (MQ=255)
                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacga                                                     <  1:1316320/93‑1 (MQ=255)
                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacga                                                     <  1:4625671/93‑1 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac                                                    >  2:1044481/1‑93 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac                                                    <  1:1552969/93‑1 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac                                                    <  2:763229/93‑1 (MQ=255)
                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg                                                   <  2:4952602/93‑1 (MQ=255)
                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg                                                   >  1:1426432/1‑93 (MQ=255)
                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg                                                   >  2:4992638/1‑93 (MQ=255)
                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg                                                   >  2:5946770/1‑93 (MQ=255)
                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg                                                   >  1:6033668/1‑93 (MQ=255)
                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGc                                                  >  1:1552563/1‑93 (MQ=255)
                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGc                                                  <  2:5168728/93‑1 (MQ=255)
                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGc                                                  <  1:426660/93‑1 (MQ=255)
                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGc                                                  <  2:4958176/93‑1 (MQ=255)
                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGcc                                                 >  2:2256472/1‑93 (MQ=255)
                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa                                               <  1:6062415/89‑1 (MQ=255)
                         tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaa                                              >  2:5083323/4‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaag                                            <  2:2846134/92‑2 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt                                            >  1:6035252/2‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt                                            >  2:4617817/2‑93 (MQ=255)
                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATg                                           <  1:5923994/93‑1 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc                                          <  2:1736980/93‑1 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc                                          <  1:3789541/93‑1 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc                                          >  2:5226036/1‑93 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt                                         <  1:6098058/93‑1 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt                                         >  1:2452004/1‑93 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt                                         >  1:82540/1‑93 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt                                         <  1:1120950/93‑1 (MQ=255)
                                                            actgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg         <  1:5532678/89‑1 (MQ=255)
                                                            actgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg         <  2:3150136/89‑1 (MQ=255)
                                                             ctgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc        >  1:1890430/4‑93 (MQ=255)
                                                              tgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCa       >  2:2729629/3‑93 (MQ=255)
                                                              tgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCa       <  2:2266755/91‑1 (MQ=255)
                                                              tgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGTAAGGCTGCATTATGGGCa       >  1:3112540/3‑93 (MQ=255)
                                                              cgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCa       <  2:2410320/91‑1 (MQ=255)
                                                               gggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAg      <  2:2017840/92‑1 (MQ=255)
                                                               gggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAg      <  2:1036539/92‑1 (MQ=255)
                                                                ggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt     <  2:4972239/93‑1 (MQ=255)
                                                                ggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt     <  2:96549/93‑1 (MQ=255)
                                                                 gCTACAGCCTGAAGGATTGTATTGGG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg    <  2:691933/93‑1 (MQ=255)
                                                                  cTACAGCCTGAAGGATTGTATTGCGTTTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc   <  1:3519216/93‑2 (MQ=255)
                                                                  cTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc   <  2:4043434/93‑2 (MQ=255)
                                                                  cTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc   <  2:2123975/93‑2 (MQ=255)
                                                                   tACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg  >  2:3052832/1‑91 (MQ=255)
                                                                        |                                                                                           
TCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGC  >  NC_000913/4008346‑4008492

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ 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
TCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGC  >  NC_000913/4008371‑4008499
                                               |                                                                                                
caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAT                                                  >  SRR4302537.6076720/3‑94 (MQ=18)
  AAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGC                                                <  SRR4302537.5964645/94‑1 (MQ=21)
   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCT                                               >  SRR4302537.2469257/1‑94 (MQ=21)
   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCT                                               <  SRR4302537.3816140/94‑1 (MQ=21)
   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCT                                               >  SRR4302537.83144/1‑94 (MQ=21)
    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTG                                              <  SRR4302537.6139961/94‑1 (MQ=23)
    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTG                                              >  SRR4302537.3086211/1‑94 (MQ=23)
    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTG                                              <  SRR4302537.1128719/94‑1 (MQ=23)
     GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGT                                             >  SRR4302537.6155501/1‑94 (MQ=25)
        catgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGA                                          <  SRR4302537.4436875/64‑1 (MQ=42)
        catgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGA                                          <  SRR4302537.307524/64‑1 (MQ=42)
         atgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGAT                                         >  SRR4302537.5413712/30‑94 (MQ=45)
         atgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGAT                                         <  SRR4302537.1350829/65‑1 (MQ=45)
            cgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGG                                      >  SRR4302537.1785606/27‑94 (MQ=48)
            cgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC                                      >  SRR4302537.4279114/27‑94 (MQ=51)
             gctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG                                    >  SRR4302537.6097263/26‑94 (MQ=53)
             gctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG                                    <  SRR4302537.2809454/69‑1 (MQ=53)
              ctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGG                                   >  SRR4302537.1414020/25‑94 (MQ=55)
              ctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGG                                   <  SRR4302537.61936/70‑1 (MQ=55)
              ctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGG                                   >  SRR4302537.6098262/25‑94 (MQ=55)
               tggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG                                  <  SRR4302537.6075178/71‑1 (MQ=57)
               tggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG                                  <  SRR4302537.6491540/71‑1 (MQ=57)
               tggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG                                  <  SRR4302537.6491541/71‑1 (MQ=57)
                 gaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑GA                              >  SRR4302537.2918606/22‑94 (MQ=56)
                 gaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑GA                              <  SRR4302537.2289760/73‑1 (MQ=56)
                  aagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑GAA                             <  SRR4302537.3465938/74‑1 (MQ=55)
                   agcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAA                             >  SRR4302537.6220179/20‑94 (MQ=54)
                    gcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAG                            <  SRR4302537.3515312/76‑1 (MQ=53)
                    gcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAG                            >  SRR4302537.219749/19‑94 (MQ=53)
                    gcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAG                            <  SRR4302537.706654/76‑1 (MQ=53)
                     cggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaagg                            <  SRR4302537.3353664/77‑7 (MQ=57)
                      ggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggc                           >  SRR4302537.3442692/17‑87 (MQ=57)
                      ggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggc                           <  SRR4302537.2329728/78‑8 (MQ=57)
                      ggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggc                           >  SRR4302537.1843874/17‑87 (MQ=57)
                       gtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggct                          >  SRR4302537.2380942/16‑86 (MQ=57)
                       gtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggct                          <  SRR4302537.775492/79‑9 (MQ=57)
                       gtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggct                          <  SRR4302537.5497392/79‑9 (MQ=57)
                       gtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggct                          >  SRR4302537.5735591/16‑86 (MQ=57)
                        tggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggctg                         <  SRR4302537.995395/80‑10 (MQ=57)
                         ggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggctgc                        <  SRR4302537.5215171/81‑11 (MQ=57)
                           cgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAT                      >  SRR4302537.2087440/12‑94 (MQ=49)
                            gaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATT                     >  SRR4302537.1498215/11‑94 (MQ=49)
                            gaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATT                     >  SRR4302537.5536242/11‑94 (MQ=49)
                            gaagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATT                     >  SRR4302537.923093/11‑94 (MQ=49)
                             aagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACTCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTA                    >  SRR4302537.1611067/10‑94 (MQ=37)
                             aagaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTA                    <  SRR4302537.874927/85‑1 (MQ=48)
                              agaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTAT                   >  SRR4302537.3877400/9‑94 (MQ=48)
                               gaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATG                  >  SRR4302537.5945483/8‑94 (MQ=49)
                               gaaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATG                  <  SRR4302537.6080847/87‑1 (MQ=49)
                                aaactgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGG                 <  SRR4302537.2985492/88‑1 (MQ=51)
                                 AACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGG                <  SRR4302537.2846596/94‑1 (MQ=50)
                                 AACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGG                <  SRR4302537.5909300/94‑1 (MQ=50)
                                  ACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGC               <  SRR4302537.5570442/94‑1 (MQ=53)
                                  ACTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGC               >  SRR4302537.1903568/1‑94 (MQ=53)
                                   ctgGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGTAAGGCTGCATTATGGGCA              >  SRR4302537.3134280/4‑94 (MQ=57)
                                        CTACAGCCTGAAGGATTGTATTGCGTTTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcg         <  SRR4302537.3544116/94‑3 (MQ=60)
                                          ACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcgca       <  SRR4302537.2416802/94‑5 (MQ=60)
                                          ACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcgca       >  SRR4302537.2174283/1‑90 (MQ=60)
                                            agcCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC     <  SRR4302537.6354651/91‑1 (MQ=60)
                                            agcCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC     >  SRR4302537.6178180/4‑94 (MQ=60)
                                            agcCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC     <  SRR4302537.736139/91‑1 (MQ=60)
                                             CCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAG   <  SRR4302537.2539279/94‑1 (MQ=60)
                                               CTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC  >  SRR4302537.3313670/1‑94 (MQ=60)
                                               |                                                                                                
TCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGC  >  NC_000913/4008371‑4008499

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