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 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,391 +C coding (474/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3911.C98.6% 209.6 / ‑3.9 69coding (474/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (30/38);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458
                                                                                              |                                                                             
cccGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCACAAGGCCATGCGCCGGCAg                                                                                >  2:547012/3‑93 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                               <  1:2463423/93‑2 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAgc                                                                               >  1:3949962/1‑92 (MQ=255)
   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGGGGCGGTTCAAAAGGCCATGAGCTGGAAgcgg                                                                             >  2:573336/1‑90 (MQ=255)
   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                             <  2:2846331/93‑4 (MQ=255)
   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAAGCCATGCGCTGGAAgcgg                                                                             >  1:1683766/1‑90 (MQ=255)
         gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGCAGCGGTGGCG‑‑‑‑a                                                                  >  2:2584046/1‑93 (MQ=255)
         gCGAGCGCGTCAACGTCAGTTCCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                  <  1:3422101/93‑1 (MQ=255)
          cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                 <  1:3182967/93‑1 (MQ=255)
          cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                 >  2:1849334/1‑93 (MQ=255)
           gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                <  1:575495/93‑2 (MQ=255)
           gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                >  2:2756240/1‑92 (MQ=255)
           gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                <  1:3581944/93‑2 (MQ=255)
           gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                <  1:2674900/93‑2 (MQ=255)
            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  2:2329246/1‑91 (MQ=255)
            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  2:2283480/1‑91 (MQ=255)
            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  1:1052050/1‑91 (MQ=255)
              cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                             <  2:2369141/93‑5 (MQ=255)
                                              tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:2592961/93‑1 (MQ=255)
                                              tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 >  2:3254248/1‑93 (MQ=255)
                                              tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:4432993/93‑1 (MQ=255)
                                                tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:3632859/93‑3 (MQ=255)
                                                tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:3793677/93‑3 (MQ=255)
                                                 ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                              >  1:1240529/1‑90 (MQ=255)
                                                  gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:1051163/93‑5 (MQ=255)
                                                  gAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             >  1:3702178/1‑89 (MQ=255)
                                                        aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       >  2:1077404/1‑93 (MQ=255)
                                                        aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  2:229684/93‑1 (MQ=255)
                                                        aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  2:760103/93‑1 (MQ=255)
                                                         tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:3374702/93‑1 (MQ=255)
                                                         tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:1566229/93‑1 (MQ=255)
                                                         tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:1894017/93‑1 (MQ=255)
                                                          ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:401856/93‑1 (MQ=255)
                                                          ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCCGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:3408020/93‑1 (MQ=255)
                                                           gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  1:3117785/93‑1 (MQ=255)
                                                           gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  2:1467147/93‑1 (MQ=255)
                                                             gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:4425676/93‑1 (MQ=255)
                                                             gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:1119358/93‑1 (MQ=255)
                                                               gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                >  1:4457103/1‑93 (MQ=255)
                                                                cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               >  1:4576177/1‑93 (MQ=255)
                                                                cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  2:3526130/93‑1 (MQ=255)
                                                                 ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  1:2669860/93‑1 (MQ=255)
                                                                 ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  2:516542/1‑93 (MQ=255)
                                                                 ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  2:4523682/1‑93 (MQ=255)
                                                                 ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  2:4322943/93‑1 (MQ=255)
                                                                 ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  1:2958163/93‑1 (MQ=255)
                                                                  gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgggg             >  1:2074904/1‑90 (MQ=255)
                                                                  gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  2:2058391/1‑93 (MQ=255)
                                                                   cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  2:1382101/93‑1 (MQ=255)
                                                                   cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  1:1265185/93‑1 (MQ=255)
                                                                   cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            >  1:2192640/1‑93 (MQ=255)
                                                                   cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            >  2:4150932/1‑93 (MQ=255)
                                                                    gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  2:1470272/93‑1 (MQ=255)
                                                                    gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:1671891/1‑93 (MQ=255)
                                                                    gTTTCAAAAGGCCATGCGCAGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:4571105/1‑93 (MQ=255)
                                                                        caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:2234580/91‑1 (MQ=255)
                                                                        caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:3847674/91‑1 (MQ=255)
                                                                         canaaGNNCATGCGCTGGAAGNNGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:4559497/90‑1 (MQ=255)
                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:1013108/2‑93 (MQ=255)
                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:4466518/2‑93 (MQ=255)
                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      <  2:4362651/92‑1 (MQ=255)
                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:2030356/2‑93 (MQ=255)
                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:2669759/2‑93 (MQ=255)
                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:2475525/2‑93 (MQ=255)
                                                                          aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  1:4590471/1‑93 (MQ=255)
                                                                          aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  2:470735/1‑93 (MQ=255)
                                                                          aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  1:1050044/1‑93 (MQ=255)
                                                                            aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  2:1078745/93‑1 (MQ=255)
                                                                            aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  1:3864513/1‑93 (MQ=255)
                                                                            aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:2907437/93‑1 (MQ=255)
                                                                            aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:4450069/93‑1 (MQ=255)
                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:3115899/93‑1 (MQ=255)
                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:4654945/1‑93 (MQ=255)
                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:4631955/93‑1 (MQ=255)
                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:2559408/93‑1 (MQ=255)
                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:1507239/1‑93 (MQ=255)
                                                                              ggCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:3964472/93‑1 (MQ=255)
                                                                                              |                                                                             
CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCG  >  NC_000913/4008299‑4008461
                                                                                              |                                                                                
GCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCGCTAACTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAG                                                                                  <  SRR4302535.641192/94‑1 (MQ=60)
 CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGC                                                                                >  SRR4302535.3976930/1‑94 (MQ=60)
  TCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcg                                                                                <  SRR4302535.2481186/94‑3 (MQ=60)
   CGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAAGCCATGCGCTGGAAGcgg                                                                               >  SRR4302535.1695960/1‑91 (MQ=60)
       GGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAACGCCATGCGCTGGAAGcggtggc                                                                           >  SRR4302535.4250813/1‑87 (MQ=60)
        GGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcg                                                                          <  SRR4302535.485115/94‑9 (MQ=60)
          GCGAGCGCGTCAACGTCAGTTCCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AA                                                                    <  SRR4302535.3445904/94‑1 (MQ=60)
           CGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaag                                                                       <  SRR4302535.3204960/94‑12 (MQ=60)
            GATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaaga                                                                      <  SRR4302535.2694038/94‑13 (MQ=60)
            GATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaaga                                                                      <  SRR4302535.579534/94‑13 (MQ=60)
            GATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaaga                                                                      <  SRR4302535.3606788/94‑13 (MQ=60)
            GATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGcggtggcgaaga                                                                      >  SRR4302535.1059533/1‑82 (MQ=60)
              TCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcgggggcgaagaaa                                                                    >  SRR4302535.3556017/1‑80 (MQ=60)
               CGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCCTGCGCTGGAAGcggtgggggagaaac                                                                   >  SRR4302535.4128343/1‑79 (MQ=60)
               CGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGGggggggggagaaac                                                                  >  SRR4302535.1786769/1‑80 (MQ=60)
                GCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaact                                                                  <  SRR4302535.3044760/94‑17 (MQ=60)
                GCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaact                                                                  <  SRR4302535.1147712/94‑17 (MQ=60)
                GCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaact                                                                  <  SRR4302535.762486/94‑17 (MQ=60)
                GCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaact                                                                  >  SRR4302535.966024/1‑78 (MQ=60)
                 CGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactg                                                                 <  SRR4302535.4048510/94‑18 (MQ=60)
                  GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgg                                                                <  SRR4302535.3472692/94‑19 (MQ=60)
                  GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgg                                                                <  SRR4302535.2166428/94‑19 (MQ=60)
                  GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgg                                                                >  SRR4302535.4344823/1‑76 (MQ=60)
                  GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgg                                                                <  SRR4302535.1327585/94‑19 (MQ=60)
                   TCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactggg                                                               <  SRR4302535.3848928/94‑20 (MQ=60)
                    CAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGcggtggcgtagaaactgggg                                                              >  SRR4302535.1566742/1‑74 (MQ=60)
                     AACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggct                                                             >  SRR4302535.2113793/1‑73 (MQ=60)
                     AACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggct                                                             <  SRR4302535.4037668/94‑22 (MQ=60)
                     AACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGagggggcgaagacactggcga                                                             >  SRR4302535.2629781/1‑73 (MQ=60)
                        GTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctaca                                                          <  SRR4302535.2466397/94‑25 (MQ=60)
                        GTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctaca                                                          <  SRR4302535.4352304/94‑25 (MQ=60)
                        GTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctaca                                                          <  SRR4302535.3178487/94‑25 (MQ=60)
                         GCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacag                                                         <  SRR4302535.4169011/94‑26 (MQ=60)
                          CAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcgggggcgaagaaactgggctacagc                                                        >  SRR4302535.4097557/1‑68 (MQ=60)
                            GTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcct                                                      <  SRR4302535.2446853/94‑29 (MQ=60)
                            GTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcct                                                      <  SRR4302535.1356051/94‑29 (MQ=60)
                             TTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctg                                                     >  SRR4302535.3054709/1‑65 (MQ=60)
                             TTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctg                                                     <  SRR4302535.3086875/94‑30 (MQ=60)
                                CTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATACGCTGGAAGcggtggcgaagaaaatgggctacagcctgtgg                                                  >  SRR4302535.1195847/1‑62 (MQ=60)
                                 TCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcgggggcgaagaaactgggctacagcctgaagg                                                 >  SRR4302535.752788/1‑61 (MQ=60)
                                   CACCTTAACCTGTCTGGAACTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggat                                               <  SRR4302535.323955/94‑36 (MQ=60)
                                    ACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggatt                                              >  SRR4302535.4496868/1‑58 (MQ=60)
                                    ACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggatt                                              >  SRR4302535.2372067/1‑58 (MQ=60)
                                    ACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggatt                                              <  SRR4302535.2373920/94‑37 (MQ=60)
                                    ACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggatt                                              >  SRR4302535.1022596/1‑58 (MQ=60)
                                     CCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattg                                             <  SRR4302535.624054/94‑38 (MQ=60)
                                      CTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgt                                            <  SRR4302535.1939338/94‑39 (MQ=60)
                                      CTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgt                                            >  SRR4302535.664935/1‑56 (MQ=60)
                                       TTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgta                                           <  SRR4302535.4479187/94‑40 (MQ=60)
                                       TTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgta                                           <  SRR4302535.3600699/94‑40 (MQ=60)
                                        TAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtat                                          <  SRR4302535.194243/94‑41 (MQ=60)
                                         AACCTGTCTGGAAGTGATGGCGGGCGGCGGTTTAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactggggtatagcctgaaggattgtatt                                         >  SRR4302535.3721523/1‑53 (MQ=59)
                                         AACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagccttaaggattgtatt                                         >  SRR4302535.344694/1‑53 (MQ=60)
                                         AACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtatt                                         >  SRR4302535.1910445/1‑53 (MQ=60)
                                          ACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattg                                        <  SRR4302535.4534257/94‑43 (MQ=60)
                                          ACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattg                                        <  SRR4302535.1148870/94‑43 (MQ=60)
                                           CCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgc                                       <  SRR4302535.3527012/94‑44 (MQ=60)
                                           CCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgc                                       <  SRR4302535.2364382/94‑44 (MQ=60)
                                            CTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcg                                      <  SRR4302535.3237860/94‑45 (MQ=60)
                                            CTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcg                                      >  SRR4302535.3966727/1‑50 (MQ=60)
                                              GTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtt                                    <  SRR4302535.3125872/94‑47 (MQ=60)
                                              GTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcggt                                    >  SRR4302535.125403/1‑48 (MQ=60)
                                               TCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgttt                                   <  SRR4302535.4463369/94‑48 (MQ=60)
                                               TCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgttt                                   <  SRR4302535.2611590/94‑48 (MQ=60)
                                                 TGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttgg                                 <  SRR4302535.3819833/94‑50 (MQ=60)
                                                 TGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttgg                                 <  SRR4302535.3658036/94‑50 (MQ=60)
                                                 TGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttgg                                 >  SRR4302535.1249385/1‑45 (MQ=60)
                                                  GGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggt                                >  SRR4302535.3727819/1‑44 (MQ=57)
                                                   GAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtg                               <  SRR4302535.1058638/94‑52 (MQ=60)
                                                    AAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtga                              <  SRR4302535.1591254/94‑53 (MQ=60)
                                                    AAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtga                              <  SRR4302535.4195952/94‑53 (MQ=60)
                                                     AGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcgggggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtgcg                             >  SRR4302535.3132113/1‑41 (MQ=60)
                                                       TGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtgacgg                           >  SRR4302535.3700461/1‑39 (MQ=60)
                                                        GATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtgacggg                          <  SRR4302535.222542/94‑57 (MQ=56)
                                                        GATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtgacggg                          <  SRR4302535.1613095/94‑57 (MQ=56)
                                                        GATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtgacggg                          <  SRR4302535.3068576/94‑57 (MQ=56)
                                                                          caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         >  SRR4302535.4497158/3‑94 (MQ=18)
                                                                          caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         >  SRR4302535.2044988/3‑94 (MQ=18)
                                                                          caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         <  SRR4302535.2250476/92‑1 (MQ=18)
                                                                          caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT         >  SRR4302535.2493380/3‑94 (MQ=18)
                                                                          AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG        >  SRR4302535.4621718/1‑94 (MQ=19)
                                                                          AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG        >  SRR4302535.1057512/1‑94 (MQ=19)
                                                                            AAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGC       >  SRR4302535.3891097/1‑94 (MQ=21)
                                                                             AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT      <  SRR4302535.4480591/94‑1 (MQ=21)
                                                                             AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT      >  SRR4302535.4686452/1‑94 (MQ=21)
                                                                             AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT      <  SRR4302535.2927986/94‑1 (MQ=21)
                                                                              AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG     <  SRR4302535.4663384/94‑1 (MQ=23)
                                                                              AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG     <  SRR4302535.3137592/94‑1 (MQ=23)
                                                                               GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT    >  SRR4302535.4137075/1‑94 (MQ=25)
                                                                                 CCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCG  <  SRR4302535.3291791/94‑1 (MQ=28)
                                                                                              |                                                                                
GCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCG  >  NC_000913/4008299‑4008461

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