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,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% 297.6 / NA 92A165P (GCT→CCT) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (40/52);  total (40/52)

TCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AA‑‑‑‑‑CTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGC  >  NC_000913/4008346‑4008492
                                                                             |                                                                                          
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                                                      >  2:5367497/1‑93 (MQ=255)
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                                                      <  1:3067422/93‑1 (MQ=255)
 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                                                                     <  1:3423220/93‑2 (MQ=255)
 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                                                                     >  2:2690426/1‑92 (MQ=255)
 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                                                                     >  2:2607369/1‑92 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                                                                    >  1:5995372/1‑91 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                                                                    >  2:4708075/1‑91 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                                                                    >  1:5951006/1‑91 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                                                                    >  2:1156237/1‑91 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                                                                    <  1:1051706/93‑3 (MQ=255)
   ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttgg                                                                   >  2:630739/1‑90 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                                                           <  2:5970655/93‑1 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                                                           <  2:349093/93‑1 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                                                           <  1:100493/93‑1 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                                                           >  2:6045798/1‑93 (MQ=255)
            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt                                                          <  2:2148939/93‑1 (MQ=255)
            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AAGAAACCGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt                                                          <  2:1686766/93‑1 (MQ=255)
             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg                                                         <  2:2236559/93‑1 (MQ=255)
             gcgggcggCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg                                                         <  2:790947/93‑1 (MQ=255)
              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGa                                                        <  1:5478423/93‑1 (MQ=255)
               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                                                       <  1:2635131/93‑1 (MQ=255)
               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                                                       <  2:3445771/93‑1 (MQ=255)
                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                                                      >  1:3993000/1‑93 (MQ=255)
                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                                                      <  2:2772791/93‑1 (MQ=255)
                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                                                      <  2:331868/93‑1 (MQ=255)
                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                                                      <  1:4957056/93‑1 (MQ=255)
                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg                                                     <  1:2957114/93‑1 (MQ=255)
                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg                                                     >  1:4814006/1‑93 (MQ=255)
                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacga                                                    <  1:2546223/93‑1 (MQ=255)
                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacga                                                    <  2:2166732/93‑1 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTT‑‑GTGACGGGATGAacgac                                                   >  2:5773184/1‑93 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac                                                   >  2:3119511/1‑93 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac                                                   <  2:4279944/93‑1 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac                                                   <  2:4498313/93‑1 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac                                                   >  2:5281397/1‑93 (MQ=255)
                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg                                                  >  2:4300379/1‑93 (MQ=255)
                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGc                                                 >  2:5938791/1‑93 (MQ=255)
                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGcc                                                <  2:4808827/93‑1 (MQ=255)
                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGcc                                                <  1:4418546/93‑1 (MQ=255)
                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa                                              <  2:2580522/89‑1 (MQ=255)
                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa                                              >  1:6270373/5‑93 (MQ=255)
                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa                                              >  2:5932837/5‑93 (MQ=255)
                         tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaa                                             >  2:476989/4‑93 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa                                            <  2:723500/91‑1 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa                                            <  2:808001/91‑1 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa                                            <  1:4553544/91‑1 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa                                            >  1:4825534/3‑93 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa                                            <  1:5766444/91‑1 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa                                            >  1:4825533/3‑93 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa                                            <  1:3434171/91‑1 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa                                            <  1:995345/91‑1 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt                                           >  2:4130459/2‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt                                           >  2:3678743/2‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt                                           >  1:5825289/2‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt                                           <  2:2946514/92‑1 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt                                           >  2:1526869/2‑93 (MQ=255)
                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGCAATg                                          >  2:2085646/1‑93 (MQ=255)
                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATg                                          >  1:290367/1‑93 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc                                         <  1:2935716/93‑1 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc                                         >  2:3848605/1‑93 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc                                         <  1:6192030/93‑1 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt                                        <  1:3735094/93‑1 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt                                        <  1:1197056/93‑1 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt                                        >  2:6045723/1‑93 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt                                        <  2:193059/93‑1 (MQ=255)
                               ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg                                       <  1:5385678/93‑1 (MQ=255)
                               ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg                                       <  1:241840/93‑1 (MQ=255)
                               ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg                                       <  2:2269624/93‑1 (MQ=255)
                               ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg                                       <  1:4092752/93‑1 (MQ=255)
                                                                 actgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATggg         <  1:1993023/89‑1 (MQ=255)
                                                                 actgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATggg         <  1:2656113/89‑1 (MQ=255)
                                                                 actgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATggg         >  1:2688598/5‑93 (MQ=255)
                                                                 actgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATggg         >  2:5517684/5‑93 (MQ=255)
                                                                 actgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGTCTGCATTATggg         >  1:3118034/5‑93 (MQ=255)
                                                                  ctgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        <  2:1450881/90‑1 (MQ=255)
                                                                  ctgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        >  2:5902450/4‑93 (MQ=255)
                                                                  ctgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        <  1:844949/90‑1 (MQ=255)
                                                                  ctgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        >  2:2214314/4‑93 (MQ=255)
                                                                  ctgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        <  1:4035461/90‑1 (MQ=255)
                                                                  ctgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        >  1:4618390/4‑93 (MQ=255)
                                                                   tgggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCa       >  1:1750273/3‑93 (MQ=255)
                                                                    gggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAg      <  2:5597564/92‑1 (MQ=255)
                                                                    gggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAg      >  2:2095688/2‑93 (MQ=255)
                                                                     ggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt     <  2:236929/93‑1 (MQ=255)
                                                                     ggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt     >  1:1525946/1‑93 (MQ=255)
                                                                     ggCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt     <  2:3176408/93‑1 (MQ=255)
                                                                       cTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc   >  1:720482/1‑92 (MQ=255)
                                                                       cTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc   <  2:2350087/93‑2 (MQ=255)
                                                                       cTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc   >  1:6162543/1‑92 (MQ=255)
                                                                       cTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc   <  2:3528898/93‑2 (MQ=255)
                                                                       cTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc   <  2:6200703/93‑2 (MQ=255)
                                                                       cTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc   <  1:5539427/93‑2 (MQ=255)
                                                                        tACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg  >  2:2417089/1‑91 (MQ=255)
                                                                        tACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg  >  1:5748970/1‑91 (MQ=255)
                                                                             |                                                                                          
TCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AA‑‑‑‑‑CTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGC  >  NC_000913/4008346‑4008492

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
TTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGC  >  NC_000913/4008370‑4008499
                                                |                                                                                               
tcaAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAA                                                    >  SRR4302533.4854776/4‑94 (MQ=16)
tcaAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAA                                                    >  SRR4302533.4854775/4‑94 (MQ=16)
 caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                   <  SRR4302533.1001159/92‑1 (MQ=18)
 caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                   <  SRR4302533.5801257/92‑1 (MQ=18)
 caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                   >  SRR4302533.5860484/3‑94 (MQ=18)
 caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                   <  SRR4302533.4581160/92‑1 (MQ=18)
 caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                   <  SRR4302533.3455385/92‑1 (MQ=18)
  AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG                                                  >  SRR4302533.292135/1‑94 (MQ=19)
    AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT                                                <  SRR4302533.2953665/94‑1 (MQ=21)
    AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT                                                <  SRR4302533.6229189/94‑1 (MQ=21)
     aggccatgcgcgggaagcggttgcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG                                               <  SRR4302533.433927/60‑1 (MQ=34)
     AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG                                               >  SRR4302533.2757655/1‑94 (MQ=23)
     AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG                                               <  SRR4302533.3758008/94‑1 (MQ=23)
     AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG                                               <  SRR4302533.1204246/94‑1 (MQ=23)
     AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG                                               >  SRR4302533.49854/1‑94 (MQ=23)
      GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT                                              <  SRR4302533.5418036/94‑1 (MQ=25)
      GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT                                              <  SRR4302533.243338/94‑1 (MQ=25)
      GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT                                              <  SRR4302533.4117665/94‑1 (MQ=25)
       GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC                                             <  SRR4302533.20499/94‑1 (MQ=26)
       GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC                                             <  SRR4302533.3654174/94‑1 (MQ=26)
       GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC                                             <  SRR4302533.5623384/94‑1 (MQ=26)
       GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC                                             >  SRR4302533.39867/1‑94 (MQ=26)
       GCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTC                                             <  SRR4302533.1174613/94‑1 (MQ=26)
         catgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGA                                           >  SRR4302533.3672547/31‑94 (MQ=42)
         catgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGA                                           >  SRR4302533.538534/31‑94 (MQ=42)
          atgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGAT                                          >  SRR4302533.835372/30‑94 (MQ=45)
          atgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGAT                                          <  SRR4302533.159236/65‑1 (MQ=45)
           tgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATG                                         <  SRR4302533.227141/66‑1 (MQ=47)
            gcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGG                                        <  SRR4302533.5769847/67‑1 (MQ=49)
             cgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC                                       <  SRR4302533.3219274/68‑1 (MQ=51)
              gctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG                                     >  SRR4302533.2823135/26‑94 (MQ=53)
              gctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG                                     <  SRR4302533.1249754/69‑1 (MQ=53)
               ctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGG                                    >  SRR4302533.3414928/25‑94 (MQ=55)
                tggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG                                   >  SRR4302533.4347199/24‑94 (MQ=57)
                tggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG                                   <  SRR4302533.725472/71‑1 (MQ=57)
                 ggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑G                                <  SRR4302533.1492162/72‑1 (MQ=56)
                  gaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑GA                               <  SRR4302533.1488020/73‑1 (MQ=56)
                  gaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑GA                               <  SRR4302533.4639773/73‑1 (MQ=56)
                  gaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑GA                               >  SRR4302533.570631/22‑94 (MQ=56)
                   aagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑GAA                              >  SRR4302533.3519055/21‑94 (MQ=55)
                     gcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAG                             <  SRR4302533.2615802/76‑1 (MQ=53)
                      cggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaagg                             <  SRR4302533.4427539/77‑7 (MQ=57)
                       ggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggc                            >  SRR4302533.1248457/17‑87 (MQ=57)
                       ggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggc                            <  SRR4302533.6042832/78‑8 (MQ=57)
                        gtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggct                           >  SRR4302533.5714838/16‑86 (MQ=57)
                        gtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggct                           <  SRR4302533.1951459/79‑9 (MQ=57)
                         tggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggctg                          >  SRR4302533.43879/15‑85 (MQ=57)
                         tggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggctg                          <  SRR4302533.355138/80‑10 (MQ=57)
                            cgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAT                      >  SRR4302533.4625322/12‑94 (MQ=49)
                            cgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAT                      >  SRR4302533.2915384/12‑94 (MQ=49)
                             gaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATT                     >  SRR4302533.3589461/11‑94 (MQ=49)
                               agaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTAT                   >  SRR4302533.3330090/9‑94 (MQ=48)
                               agaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTAT                   <  SRR4302533.2666254/86‑1 (MQ=48)
                                gaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATG                  <  SRR4302533.1052833/87‑1 (MQ=49)
                                 aaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGG                 <  SRR4302533.4657144/88‑1 (MQ=51)
                                 aaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGG                 <  SRR4302533.2348022/88‑1 (MQ=51)
                                  AACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGTCTGCATTATGGG                >  SRR4302533.3137206/1‑94 (MQ=48)
                                  AACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGG                >  SRR4302533.2705082/1‑94 (MQ=50)
                                  AACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGG                <  SRR4302533.1644244/94‑1 (MQ=50)
                                   ACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGC               >  SRR4302533.4646378/1‑94 (MQ=53)
                                   ACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGC               <  SRR4302533.2672427/94‑1 (MQ=53)
                                   ACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGC               <  SRR4302533.2005247/94‑1 (MQ=53)
                                    ctgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCA              <  SRR4302533.4060026/91‑1 (MQ=57)
                                    ctgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCA              <  SRR4302533.849893/91‑1 (MQ=57)
                                    ctgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCA              >  SRR4302533.1760930/4‑94 (MQ=57)
                                     GGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGT            >  SRR4302533.1535130/1‑94 (MQ=60)
                                        GCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGC          >  SRR4302533.6199597/1‑94 (MQ=60)
                                        GCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGC          >  SRR4302533.724731/1‑94 (MQ=60)
                                         CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcg         <  SRR4302533.5572877/94‑3 (MQ=60)
                                         CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcg         >  SRR4302533.5783651/1‑92 (MQ=60)
                                           ACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcgca       <  SRR4302533.6277876/94‑5 (MQ=60)
                                             agcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC     >  SRR4302533.1025366/4‑94 (MQ=60)
                                             agcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC     >  SRR4302533.5693748/4‑94 (MQ=60)
                                             agcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC     <  SRR4302533.5237275/91‑1 (MQ=60)
                                             agcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC     >  SRR4302533.1562424/4‑94 (MQ=60)
                                              gcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCA    <  SRR4302533.3114813/92‑1 (MQ=60)
                                              gcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCA    >  SRR4302533.97895/3‑94 (MQ=60)
                                              CCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAG   >  SRR4302533.2061001/1‑94 (MQ=60)
                                                CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAATGCTGCATTATGGGCAGTGCGCACCAGC  >  SRR4302533.2649175/1‑94 (MQ=60)
                                                CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC  <  SRR4302533.5893172/94‑1 (MQ=60)
                                                CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC  >  SRR4302533.5281497/1‑94 (MQ=60)
                                                CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC  >  SRR4302533.6140799/1‑94 (MQ=60)
                                                CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC  <  SRR4302533.5240768/94‑1 (MQ=60)
                                                CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC  >  SRR4302533.1691103/1‑94 (MQ=60)
                                                |                                                                                               
TTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGC  >  NC_000913/4008370‑4008499

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