Sample Resequencing Stats

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

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,554,088 +C coding (2/624 nt) yhgKJ ← RNA phosphage cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,554,0881.C99.6% 797.5 / ‑3.5 229coding (2/624 nt)yhgKJRNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (120/108);  total (120/109)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 259 total aligned reads displayed.

CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTT  >  NC_000913/3553996‑3554173
                                                                                             |                                                                                      
 ggCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgc                                                                                        <  2:6040678/93‑2 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                       <  2:3527553/93‑3 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                       <  2:1765783/93‑3 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                       <  1:2025265/93‑3 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                       <  2:5381457/93‑3 (MQ=255)
   gaGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                      <  2:4916079/92‑4 (MQ=255)
       gTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcgg                                                                                  <  2:3263300/93‑1 (MQ=255)
        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                 <  1:122036/93‑1 (MQ=255)
        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                 >  2:2303384/1‑93 (MQ=255)
         aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGgcggcg                                                                                >  1:104675/1‑93 (MQ=255)
           aaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGcc                                                                              >  1:1075909/1‑93 (MQ=255)
           aaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGcc                                                                              >  1:391725/1‑93 (MQ=255)
             aCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCgc                                                                            >  1:1940846/1‑92 (MQ=255)
              ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                                           <  1:3181859/93‑1 (MQ=255)
               cGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                                          >  1:4692426/1‑93 (MQ=255)
                gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCCGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAccg                                                                         >  2:2002296/1‑91 (MQ=255)
                gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACgg                                                                         <  1:5659741/93‑1 (MQ=255)
                 tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                                        <  2:4560786/93‑1 (MQ=255)
                  ggCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGt                                                                       <  1:6018572/93‑1 (MQ=255)
                  ggCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGt                                                                       <  2:5508680/93‑1 (MQ=255)
                   gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGtt                                                                      >  1:2417162/1‑93 (MQ=255)
                     gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTAt                                                                    >  2:2477630/1‑93 (MQ=255)
                      tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                   >  1:1846628/1‑93 (MQ=255)
                       aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCg                                                                  <  1:1462447/93‑1 (MQ=255)
                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                 <  2:2729890/93‑1 (MQ=255)
                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                 <  1:462468/93‑1 (MQ=255)
                          aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtg                                                               >  2:5103459/1‑92 (MQ=255)
                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                              >  1:5877899/1‑91 (MQ=255)
                            cGTGGTTTGCTGATTAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgcc                                                             >  2:2806885/1‑90 (MQ=255)
                              tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                         >  2:1264731/1‑93 (MQ=255)
                              tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                         >  1:1129754/1‑93 (MQ=255)
                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAACTCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                        <  1:6116065/93‑2 (MQ=255)
                                  ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccgct                                                     >  1:4881201/1‑89 (MQ=255)
                                    gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCa                                                    >  1:3595083/1‑93 (MQ=255)
                                     cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                   >  1:273616/1‑93 (MQ=255)
                                     cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                   <  2:140838/93‑1 (MQ=255)
                                      tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                  <  1:3614906/93‑1 (MQ=255)
                                         tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGGCGC‑CGCTCACTTCa                                               >  2:72051/1‑93 (MQ=255)
                                         tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCa                                               >  2:2501438/1‑93 (MQ=255)
                                         tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATACGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCa                                               >  1:6136700/1‑93 (MQ=255)
                                          gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                              >  1:5307856/1‑93 (MQ=255)
                                          gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                              >  1:4330893/1‑93 (MQ=255)
                                          gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGAGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                              >  2:1520486/1‑93 (MQ=255)
                                           aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAac                                             >  2:4769927/1‑93 (MQ=255)
                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                            >  1:5152306/1‑93 (MQ=255)
                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                            <  2:2872012/93‑1 (MQ=255)
                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                            >  1:4044213/1‑93 (MQ=255)
                                             ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAacac                                           >  2:4200314/1‑93 (MQ=255)
                                              tCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcg                                          >  1:2621689/1‑93 (MQ=255)
                                               ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgc                                         >  2:5132071/1‑93 (MQ=255)
                                               ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgc                                         >  2:2386483/1‑93 (MQ=255)
                                                 tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGa                                       >  2:5888844/1‑93 (MQ=255)
                                                 tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGa                                       >  1:5062295/1‑93 (MQ=255)
                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                      <  2:908857/93‑1 (MQ=255)
                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                      <  1:1559320/93‑1 (MQ=255)
                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                      >  1:3272259/1‑93 (MQ=255)
                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                      >  1:4607025/1‑93 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                     >  1:5414023/1‑93 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                     <  1:2528060/93‑1 (MQ=255)
                                                    tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                    <  1:3819659/93‑1 (MQ=255)
                                                     ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                   <  1:557569/93‑1 (MQ=255)
                                                      tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTc                                  <  2:3652803/93‑1 (MQ=255)
                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                <  1:218442/93‑1 (MQ=255)
                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                >  1:4117577/1‑93 (MQ=255)
                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                               >  2:3959039/1‑93 (MQ=255)
                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                               <  1:5613411/93‑1 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                              >  2:3126761/1‑93 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                              <  1:4747339/93‑1 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACCCGAAGGTCCATc                              >  1:3791410/1‑93 (MQ=255)
                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                             >  2:1935657/1‑93 (MQ=255)
                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                             >  1:975540/1‑93 (MQ=255)
                                                            gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                            >  1:3109289/1‑93 (MQ=255)
                                                             cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                           >  2:2525851/1‑93 (MQ=255)
                                                             cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                           <  1:2057942/93‑1 (MQ=255)
                                                               ggcgCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACGTCAACACGCGAAGGTCCATCGGCaa                         <  2:4243033/90‑1 (MQ=255)
                                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGGCCATCGGCaa                         >  1:5780353/1‑93 (MQ=255)
                                                                cGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                        <  2:369143/93‑1 (MQ=255)
                                                                cGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                        <  2:3985804/93‑1 (MQ=255)
                                                                 ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                       >  2:4208693/1‑93 (MQ=255)
                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGCAGGTCCATCGGCAAAcc                      >  1:3325581/1‑93 (MQ=255)
                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCCAAcc                      >  2:2364643/1‑93 (MQ=255)
                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                      >  2:4307065/1‑93 (MQ=255)
                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                      <  2:3856358/93‑1 (MQ=255)
                                                                   cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                     <  1:5695375/93‑1 (MQ=255)
                                                                     ccAGCCCCCGGCGGAGAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                   <  2:619616/93‑1 (MQ=255)
                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                  >  1:1806606/1‑93 (MQ=255)
                                                                       agaccccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                 <  1:367236/89‑1 (MQ=255)
                                                                        gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                >  2:1099353/1‑93 (MQ=255)
                                                                         cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc               >  2:5321643/1‑93 (MQ=255)
                                                                          ccccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg              <  2:3572140/92‑1 (MQ=255)
                                                                          aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg              >  2:5538665/1‑93 (MQ=255)
                                                                            ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg            >  1:879338/1‑93 (MQ=255)
                                                                            ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg            >  2:1023053/1‑93 (MQ=255)
                                                                             cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc           >  2:1067786/1‑93 (MQ=255)
                                                                             cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc           <  2:4634965/93‑1 (MQ=255)
                                                                              ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca          <  1:5822580/93‑1 (MQ=255)
                                                                              ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca          <  1:4446806/93‑1 (MQ=255)
                                                                                 caaTAAAATCCGCAGGCGGCGCCGAAGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca       <  1:5201813/91‑1 (MQ=255)
                                                                                  gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc      <  2:839045/93‑1 (MQ=255)
                                                                                      aaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt  >  1:2376141/1‑93 (MQ=255)
                                                                                             |                                                                                      
CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTT  >  NC_000913/3553996‑3554173

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑G‑AGGCGGCGCCGACGGGTTATCGGT‑‑CGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAC  >  NC_000913/3553995‑3554181
                                                                                              |                                                                                                 
CCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑G                                                                                                   >  SRR4302532.500068/1‑94 (MQ=60)
  GGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑Gca                                                                                                 <  SRR4302532.395875/94‑3 (MQ=60)
   GCAGGGTAAAACCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑G‑AAG                                                                                               <  SRR4302532.189337/94‑1 (MQ=60)
   GCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑Gcag                                                                                                <  SRR4302532.2092785/94‑4 (MQ=60)
     AGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGC                                                                                              <  SRR4302532.5781933/94‑1 (MQ=60)
      GGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCG                                                                                             <  SRR4302532.3105339/94‑1 (MQ=60)
      GGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCG                                                                                             <  SRR4302532.4875756/94‑1 (MQ=60)
         TAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCG                                                                                          <  SRR4302532.5468937/94‑1 (MQ=60)
         TAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCG                                                                                          <  SRR4302532.126187/94‑1 (MQ=60)
         TAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCG                                                                                          >  SRR4302532.108288/1‑94 (MQ=60)
         TAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCG                                                                                          >  SRR4302532.108287/1‑94 (MQ=60)
           AAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCC                                                                                        >  SRR4302532.1112016/1‑94 (MQ=60)
           AAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCC                                                                                        >  SRR4302532.405429/1‑94 (MQ=60)
             AACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGC                                                                                      >  SRR4302532.2005384/1‑94 (MQ=60)
               CCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACG                                                                                    <  SRR4302532.3287525/94‑1 (MQ=60)
               CCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACG                                                                                    >  SRR4302532.4852466/1‑94 (MQ=60)
               CCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACG                                                                                    <  SRR4302532.393632/94‑1 (MQ=60)
                CGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGG                                                                                   <  SRR4302532.4269110/94‑1 (MQ=60)
                 GTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGG                                                                                  <  SRR4302532.5856917/94‑1 (MQ=60)
                 GTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGG                                                                                  <  SRR4302532.1129464/94‑1 (MQ=60)
                  TGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGT                                                                                 <  SRR4302532.4576462/94‑1 (MQ=60)
                  TGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGT                                                                                 <  SRR4302532.4609547/94‑1 (MQ=60)
                   GGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTT                                                                                <  SRR4302532.6226414/94‑1 (MQ=60)
                   GGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTT                                                                                >  SRR4302532.2497540/1‑94 (MQ=60)
                   GGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGT                                                                                >  SRR4302532.489440/1‑94 (MQ=60)
                     CGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTAT                                                                              >  SRR4302532.3456067/1‑94 (MQ=60)
                     CGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGTAT                                                                              >  SRR4302532.4223399/1‑94 (MQ=60)
                     CGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCAAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACcgggtaa                                                                              >  SRR4302532.1004991/1‑87 (MQ=60)
                      GTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATC                                                                             >  SRR4302532.1907932/1‑94 (MQ=60)
                      GTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGTATC                                                                             >  SRR4302532.2046984/1‑94 (MQ=60)
                        AACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGG                                                                           <  SRR4302532.1511725/94‑1 (MQ=60)
                         ACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT                                                                          <  SRR4302532.478559/94‑1 (MQ=60)
                           AGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGTgc                                                                        >  SRR4302532.6081514/1‑92 (MQ=60)
                           AGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGTgc                                                                        >  SRR4302532.1523809/1‑92 (MQ=60)
                             CGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCACCACCCGACGGATAAAATCC‑GCCGCCGGCGCCGACGGGTTATCGtagccg                                                                      >  SRR4302532.2179184/1‑88 (MQ=60)
                              GTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCG‑C                                                                 >  SRR4302532.1167454/1‑94 (MQ=60)
                               TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGTGCCGCC                                                                    >  SRR4302532.5945758/1‑94 (MQ=60)
                               TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGTGCCGCC                                                                    <  SRR4302532.4341421/94‑1 (MQ=60)
                                GGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAACTCC‑GCAGGCGGCGCCGACGGGTTATCGGTGCCGCCG                                                                   <  SRR4302532.6327140/94‑1 (MQ=60)
                                GGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGTGCCGCCG                                                                   <  SRR4302532.3380335/94‑1 (MQ=60)
                                  TTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGTgccgccgct                                                                 >  SRR4302532.5049044/1‑85 (MQ=60)
                                   TTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGTgccgccgctc                                                                >  SRR4302532.2431302/1‑84 (MQ=60)
                                   TTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGTgccgccgctc                                                                >  SRR4302532.2431301/1‑84 (MQ=60)
                                    TGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCA                                                            >  SRR4302532.3715307/1‑94 (MQ=60)
                                     GCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCAC                                                           <  SRR4302532.4186234/94‑1 (MQ=60)
                                     GCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCAC                                                           >  SRR4302532.283092/1‑94 (MQ=60)
                                      CTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACT                                                          <  SRR4302532.3724081/94‑1 (MQ=60)
                                      CTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCAGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACT                                                          >  SRR4302532.4807665/1‑94 (MQ=60)
                                       TGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGTgctgccgctcactt                                                            >  SRR4302532.2337534/1‑80 (MQ=60)
                                       TGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTT                                                         <  SRR4302532.3735806/94‑1 (MQ=60)
                                       TGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTT                                                         >  SRR4302532.6000525/1‑94 (MQ=60)
                                       TGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGACGACGGGTTGTCGGT‑‑‑GCCGCCGCTCACTT                                                         >  SRR4302532.2218508/1‑94 (MQ=60)
                                        GATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTC                                                        <  SRR4302532.286394/94‑1 (MQ=60)
                                         ATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCA                                                       >  SRR4302532.2536094/1‑94 (MQ=60)
                                         ATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATAC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCA                                                       >  SRR4302532.6348394/1‑94 (MQ=60)
                                          TGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAA                                                      >  SRR4302532.5491375/1‑94 (MQ=60)
                                          TGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAA                                                      >  SRR4302532.4477207/1‑94 (MQ=60)
                                           GAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAAC                                                     >  SRR4302532.3786782/1‑94 (MQ=60)
                                            AATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACA                                                    >  SRR4302532.1130481/1‑94 (MQ=60)
                                            AATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACA                                                    >  SRR4302532.5329556/1‑94 (MQ=60)
                                            AATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACA                                                    >  SRR4302532.4180387/1‑94 (MQ=60)
                                             ATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACAC                                                   <  SRR4302532.5994311/94‑1 (MQ=60)
                                             ATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACAC                                                   >  SRR4302532.380866/1‑94 (MQ=60)
                                              TTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACG                                                  >  SRR4302532.2709467/1‑94 (MQ=60)
                                              TTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACG                                                  >  SRR4302532.4081427/1‑94 (MQ=60)
                                               TCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGC                                                 >  SRR4302532.1415778/1‑94 (MQ=60)
                                                 CTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGA                                               >  SRR4302532.1747960/1‑94 (MQ=60)
                                                 CTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGA                                               >  SRR4302532.5236094/1‑94 (MQ=60)
                                                  TATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAA                                              >  SRR4302532.4764109/1‑94 (MQ=60)
                                                  TATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAA                                              >  SRR4302532.3381084/1‑94 (MQ=60)
                                                  TATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAA                                              >  SRR4302532.3336897/1‑94 (MQ=60)
                                                   ATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAG                                             <  SRR4302532.1612057/94‑1 (MQ=60)
                                                   ATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAG                                             <  SRR4302532.1306136/94‑1 (MQ=60)
                                                   ATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAG                                             >  SRR4302532.5601616/1‑94 (MQ=60)
                                                    TTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGG                                            <  SRR4302532.2612416/94‑1 (MQ=60)
                                                    TTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGG                                            <  SRR4302532.3240674/94‑1 (MQ=60)
                                                     TTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGT                                           <  SRR4302532.3948611/94‑1 (MQ=60)
                                                     TTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGT                                           <  SRR4302532.2483741/94‑1 (MQ=60)
                                                      TTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTC                                          <  SRR4302532.210134/94‑1 (MQ=60)
                                                      TTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTC                                          <  SRR4302532.577012/94‑1 (MQ=60)
                                                      TTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTC                                          >  SRR4302532.2918402/1‑94 (MQ=60)
                                                        CGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCA                                        >  SRR4302532.4255777/1‑94 (MQ=60)
                                                        CGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCA                                        >  SRR4302532.5634994/1‑94 (MQ=60)
                                                        CGCCACCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCA                                        <  SRR4302532.2597435/94‑1 (MQ=60)
                                                         GCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAT                                       <  SRR4302532.4774000/94‑1 (MQ=60)
                                                         GCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAT                                       <  SRR4302532.225898/94‑1 (MQ=60)
                                                          CCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATC                                      <  SRR4302532.5209373/94‑1 (MQ=60)
                                                          CCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATC                                      <  SRR4302532.5808196/94‑1 (MQ=60)
                                                          CCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACCCGAAGGTCCATC                                      >  SRR4302532.3919333/1‑94 (MQ=60)
                                                           CAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCG                                     <  SRR4302532.6303167/94‑1 (MQ=60)
                                                           CAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCG                                     <  SRR4302532.4909568/94‑1 (MQ=60)
                                                           CAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCG                                     >  SRR4302532.1008842/1‑94 (MQ=60)
                                                            AGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGG                                    >  SRR4302532.3212424/1‑94 (MQ=60)
                                                              CAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCA                                  <  SRR4302532.2126483/94‑1 (MQ=60)
                                                              CAGCGGCCCCAGCCCCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCA                                  <  SRR4302532.1255469/94‑1 (MQ=60)
                                                               AGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAA                                 >  SRR4302532.5310616/1‑94 (MQ=60)
                                                               AGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGGCCATCGGCAA                                 >  SRR4302532.5981412/1‑94 (MQ=60)
                                                                GCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAA                                >  SRR4302532.6327782/1‑94 (MQ=60)
                                                                  GGCTCCAGCCCCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACC                              <  SRR4302532.2718328/94‑1 (MQ=60)
                                                                  GGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGCAGGTCCATCGGCAAACC                              >  SRR4302532.3436480/1‑94 (MQ=60)
                                                                   GCTCCAGCACCCGGCGGATAAAATCT‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCA                             >  SRR4302532.3894418/1‑94 (MQ=60)
                                                                    CTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCAC                            <  SRR4302532.5894063/94‑1 (MQ=60)
                                                                    CTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCAC                            <  SRR4302532.4890966/94‑1 (MQ=60)
                                                                     cccacCACCCGGCGGATAAAATCC‑GCAGGCGGGGCCGACGGGTTAGCGGT‑‑‑GCCGCCGCTAACTTCAACACGCGAAGGTCCATCGGCAAACCACA                           <  SRR4302532.329794/89‑1 (MQ=60)
                                                                     TCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                           >  SRR4302532.4761851/1‑94 (MQ=60)
                                                                     TCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                           >  SRR4302532.2667519/1‑94 (MQ=60)
                                                                     TCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                           >  SRR4302532.2770200/1‑94 (MQ=60)
                                                                     TCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                           <  SRR4302532.6130334/94‑1 (MQ=60)
                                                                     TCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                           <  SRR4302532.3664470/94‑1 (MQ=60)
                                                                     TCCAGCACCCGGCGGATAAAAGCC‑GCAGGCGGCGCCGACGGATTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                           <  SRR4302532.679527/94‑1 (MQ=60)
                                                                      CCAGCACCCGGCGGATAAAATCCTC‑AGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAG                          <  SRR4302532.1833771/94‑1 (MQ=60)
                                                                      CCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAG                          >  SRR4302532.1866899/1‑94 (MQ=60)
                                                                      CCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAG                          >  SRR4302532.5919235/1‑94 (MQ=60)
                                                                       CAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGCCGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGC                         <  SRR4302532.5866500/94‑1 (MQ=60)
                                                                       CAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGC                         <  SRR4302532.784513/94‑1 (MQ=60)
                                                                       CAGCACCCGGCGAATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGC                         <  SRR4302532.3753405/94‑1 (MQ=60)
                                                                        aggcCCCGGCGGATAAAATCG‑GCAGCCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCG                        <  SRR4302532.2059126/90‑1 (MQ=60)
                                                                        agacCCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCG                        <  SRR4302532.379997/90‑1 (MQ=60)
                                                                        AGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCG                        <  SRR4302532.5350210/94‑1 (MQ=60)
                                                                          CACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCG                      >  SRR4302532.5491309/1‑94 (MQ=60)
                                                                          CACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCG                      >  SRR4302532.6063291/1‑94 (MQ=60)
                                                                           CCCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGG                     <  SRR4302532.5853777/94‑1 (MQ=60)
                                                                           ACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGG                     <  SRR4302532.6341392/94‑1 (MQ=60)
                                                                            CCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGG                    >  SRR4302532.909408/1‑94 (MQ=60)
                                                                             CCGGCGGATAAAATCC‑GCAGTCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGC                   <  SRR4302532.2128745/94‑1 (MQ=60)
                                                                               GGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAG                 <  SRR4302532.4598361/94‑1 (MQ=60)
                                                                               GGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAG                 >  SRR4302532.5733771/1‑94 (MQ=60)
                                                                               GGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAG                 <  SRR4302532.6024780/94‑1 (MQ=60)
                                                                               GGAGGAGAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑CCCGC‑CGCTCACGGCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAG                 <  SRR4302532.3966696/94‑1 (MQ=60)
                                                                                GCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGC                >  SRR4302532.5843361/1‑94 (MQ=60)
                                                                                 CGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCC               >  SRR4302532.4825452/1‑94 (MQ=60)
                                                                                  caATAAAATCC‑GCAGGCGGCGCCGAAGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAC              <  SRR4302532.5381178/92‑1 (MQ=60)
                                                                                  GGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAC              >  SRR4302532.5569367/1‑94 (MQ=60)
                                                                                  GGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAC              >  SRR4302532.186805/1‑94 (MQ=60)
                                                                                  GGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAC              <  SRR4302532.5381177/94‑1 (MQ=60)
                                                                                   GATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACC             <  SRR4302532.6088966/94‑1 (MQ=60)
                                                                                   GATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACC             <  SRR4302532.520318/94‑1 (MQ=60)
                                                                                    ATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCG            <  SRR4302532.3390266/94‑1 (MQ=60)
                                                                                      AAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTT          >  SRR4302532.2455340/1‑94 (MQ=60)
                                                                                      AAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTT          >  SRR4302532.2361573/1‑94 (MQ=60)
                                                                                       AAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT         <  SRR4302532.1980937/94‑1 (MQ=60)
                                                                                              C‑AGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAC  >  SRR4302532.2829016/1‑94 (MQ=60)
                                                                                              C‑AGGCGGCGCCGACGGGTTATCGGT‑‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAC  >  SRR4302532.4351133/1‑94 (MQ=60)
                                                                                              |                                                                                                 
CCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑G‑AGGCGGCGCCGACGGGTTATCGGT‑‑CGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAC  >  NC_000913/3553995‑3554181

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