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
A3 F11 I1 R1 56 79.5 3383150 98.4% 3329019 146.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 336,125 +T coding (210/237 nt) PP_0278 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC336,1241.T98.1% 174.1 / ‑3.2 52coding (209/237 nt)PP_0278hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (15/36);  total (15/37)
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.81e-01

ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAG  >  NC_002947_CJ‑RC/335980‑336264
                                                                                                                                                 |                                                                                                                                             
aCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCTTGTCTCGTTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATATATTCCTGGGGGCAGCCGAttac                                                                                                                                           <  2:280024/149‑4 (MQ=255)
aCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACATAAATTCCTGGGGGCAGCCGAttac                                                                                                                                           <  4:230320/148‑4 (MQ=255)
 cccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCTACGCGGCCTACGAGCATGTCTAGATCCTGATTGGCTACATCAAGCACCTGGTGCGCGAGGGCGATTTGGAGGACGACCATAAATTCCTGGGGGCAGCCGAttacc                                                                                                                                          <  3:274879/149‑5 (MQ=255)
   cATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCtg                                                                                                                                        <  8:392130/149‑1 (MQ=255)
      cccTCAATCTCCTACGCGAAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGc                                                                                                                                     <  2:331139/149‑1 (MQ=255)
       ccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGc                                                                                                                                     <  7:361164/148‑1 (MQ=255)
        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGa                                                                                                                                   >  1:209069/1‑149 (MQ=255)
            aTCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATg                                                                                                                               <  3:410103/149‑1 (MQ=255)
            aTCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATg                                                                                                                               <  4:384944/149‑1 (MQ=255)
               tcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGc                                                                                                                             >  2:257258/1‑148 (MQ=255)
                 cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGTGGGCAGCCGATTACCTGTGCGACATGGc                                                                                                                             <  8:74906/146‑1 (MQ=255)
                 cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGcc                                                                                                                            >  1:44612/1‑147 (MQ=255)
                 cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGc                                                                                                                             <  5:277639/146‑1 (MQ=255)
                  tACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAg                                                                                                                         >  8:263334/1‑149 (MQ=255)
                    cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGgc                                                                                                                       <  6:112598/149‑1 (MQ=255)
                       gTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCt                                                                                                                    >  6:286118/1‑149 (MQ=255)
                              ccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAAc                                                                                                             <  7:113005/149‑1 (MQ=255)
                                ggCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGa                                                                                                           <  7:139695/149‑1 (MQ=255)
                                    tCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                       <  2:385739/149‑1 (MQ=255)
                                    tCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                       <  1:257258/149‑1 (MQ=255)
                                          tCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATGAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATc                                                                                                 <  1:210714/149‑1 (MQ=255)
                                          tCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATc                                                                                                 <  4:245370/149‑1 (MQ=255)
                                               gcgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCaa                                                                                            <  1:278300/149‑1 (MQ=255)
                                                 ggCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAgg                                                                                          <  1:228439/149‑1 (MQ=255)
                                                    cTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                        >  1:228570/1‑148 (MQ=255)
                                                       cGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAg                                                                                     <  3:326871/148‑1 (MQ=255)
                                                        gAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgc                                                                                   >  5:122310/1‑149 (MQ=255)
                                                           cATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCa                                                                                >  7:61468/1‑149 (MQ=255)
                                                                ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                           <  6:43493/149‑1 (MQ=255)
                                                                 tcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                           >  7:133306/1‑148 (MQ=255)
                                                                    aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAtt                                                                       <  4:41841/149‑1 (MQ=255)
                                                                    aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAtt                                                                       >  5:360776/1‑149 (MQ=255)
                                                                    aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAtt                                                                       <  6:264765/149‑1 (MQ=255)
                                                                        tGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcg                                                                   >  7:405971/1‑149 (MQ=255)
                                                                               cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAt                                                            <  7:429400/149‑1 (MQ=255)
                                                                                    tCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAttt                                                          >  4:226173/1‑146 (MQ=255)
                                                                                    tCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAttt                                                          <  3:226173/146‑1 (MQ=255)
                                                                                    tCAAGCACCTCGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgtg                                                       <  4:92904/149‑1 (MQ=255)
                                                                                      aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGTCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                      <  6:292758/148‑1 (MQ=255)
                                                                                      aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                      <  7:263334/148‑1 (MQ=255)
                                                                                       aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAgg                                                                                          >  7:178536/1‑111 (MQ=255)
                                                                                       aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAgg                                                                                          <  8:178536/111‑1 (MQ=255)
                                                                                       aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGa                                                    <  6:207148/149‑1 (MQ=255)
                                                                                                   gcgAGGGCGATATTGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACCTGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCTTTAGGCATTGAAGCGGTCATTTGTGGGAGCGGCCTTGCGt                                        <  8:236986/149‑1 (MQ=255)
                                                                                                        ggCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGa                                   <  5:286118/149‑1 (MQ=255)
                                                                                                             tatGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                              <  6:360776/149‑1 (MQ=255)
                                                                                                               tGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGc                             <  2:271839/148‑1 (MQ=255)
                                                                                                                ggaggaCGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCt                            <  4:312248/148‑1 (MQ=255)
                                                                                                                    gacgacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATCTCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCaaa                       <  1:321056/149‑1 (MQ=255)
                                                                                                                     acgacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAag                      <  2:209069/149‑1 (MQ=255)
                                                                                                                       gacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAag                      >  1:293091/1‑147 (MQ=255)
                                                                                                                            tAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGcccc               >  7:103512/1‑149 (MQ=255)
                                                                                                                                      ggggCACCCTATTACCTGTGCGACATGCCCACGGCGCTGATCAACGATACCGCCATCCCCACGGGTAAGGCGAATTAGCCATTGCCGCGGTCATTTGGGGCAGCGGCCTTGCGGCGCTATTGGGCTGCAAAGCAGCCCCGGCGATTTCa     <  8:46273/149‑1 (MQ=255)
                                                                                                                                        ggCAGCCTATTCCCTGGGCGACATGCCCAAGGCGCTGATCAACGATATCCAGACCGCCAAGGGTAAGGCCCATTAGGCATTGACGCGGTCATTTGGGGGCGCGGCCTTGCGCCGCTAATGGGCTGCAAAGCAGCCCCGGCGATTTCaga   <  2:293091/149‑1 (MQ=255)
                                                                                                                                         gccgCCGATTACCGGTGCGACATGGCCACGGCGCTCATCAACGAAATCCAGATCGCCACGGGTAAGGCGCATTAGCCATTGACGCGGTCATTTGTGGCAGCGGCCTTGCGCCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCagag  <  2:44612/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAG  >  NC_002947_CJ‑RC/335980‑336264

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCT  >  NC_002947_CJ‑RC/335980‑336272
                                                                                                                                                  |                                                                                                                                                   
ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCTACGCGGCCTACGAGCATGTCTAGATCCTGATTGGCTACATCAAGCACCTGGTGCGCGAGGGCGATTTGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCT                                                                                                                                                 <  NB501016:282:H2V7GAFX2:2:21111:24323:14721/151‑1 (MQ=60)
  CCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACACCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAAATCCTGGGGGCAGCCGATtaccagt                                                                                                                                              <  NB501016:282:H2V7GAFX2:4:21505:22511:7739/151‑8 (MQ=60)
      CCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCG                                                                                                                                            <  NB501016:282:H2V7GAFX2:4:21511:16440:20287/150‑1 (MQ=60)
       CCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGAC                                                                                                                                          >  NB501016:282:H2V7GAFX2:1:11311:2769:13937/1‑151 (MQ=60)
           AATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGG                                                                                                                                      <  NB501016:282:H2V7GAFX2:2:21312:22427:17384/151‑1 (MQ=60)
               TCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCA                                                                                                                                   >  NB501016:282:H2V7GAFX2:3:11404:17181:13748/1‑150 (MQ=60)
                CCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCA                                                                                                                                   >  NB501016:282:H2V7GAFX2:1:11107:5439:20088/1‑149 (MQ=60)
                CCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCC                                                                                                                                    <  NB501016:282:H2V7GAFX2:3:21412:10747:3136/148‑1 (MQ=60)
                             CCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACG                                                                                                                    <  NB501016:282:H2V7GAFX2:4:11507:22104:3276/151‑1 (MQ=60)
                               CGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGAT                                                                                                                  <  NB501016:282:H2V7GAFX2:4:11511:10653:9477/151‑1 (MQ=60)
                                   GTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCG                                                                                                              <  NB501016:282:H2V7GAFX2:1:21107:25439:12725/151‑1 (MQ=60)
                                         GTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATGAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCG                                                                                                        <  NB501016:282:H2V7GAFX2:1:11311:5565:19381/151‑1 (MQ=60)
                                              CGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAG                                                                                                   <  NB501016:282:H2V7GAFX2:1:21111:20362:6059/151‑1 (MQ=60)
                                                CGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGG                                                                                                 <  NB501016:282:H2V7GAFX2:1:21102:19095:16945/151‑1 (MQ=60)
                                                   CCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTA                                                                                               >  NB501016:282:H2V7GAFX2:1:21102:10744:17293/1‑150 (MQ=60)
                                                      ACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGG                                                                                            <  NB501016:282:H2V7GAFX2:2:21210:3246:11475/150‑1 (MQ=60)
                                                       CGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCG                                                                                          >  NB501016:282:H2V7GAFX2:3:11509:24152:15858/1‑151 (MQ=60)
                                                          GCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCAT                                                                                       >  NB501016:282:H2V7GAFX2:4:11410:24656:15418/1‑151 (MQ=60)
                                                                CTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGC                                                                                  >  NB501016:282:H2V7GAFX2:4:11510:3059:8440/1‑150 (MQ=60)
                                                                   GATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTG                                                                              >  NB501016:282:H2V7GAFX2:3:21602:4309:9636/1‑151 (MQ=60)
                                                                       CTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGC                                                                          >  NB501016:282:H2V7GAFX2:4:21607:12520:15692/1‑151 (MQ=60)
                                                                              GCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATT                                                                   <  NB501016:282:H2V7GAFX2:4:21611:18153:16974/151‑1 (MQ=60)
                                                                                   ATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTG                                                                 <  NB501016:282:H2V7GAFX2:2:21102:23624:16458/148‑1 (MQ=60)
                                                                                     CAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGG                                                             <  NB501016:282:H2V7GAFX2:4:21408:7260:5359/150‑1 (MQ=60)
                                                                                      AAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGG                                                                                                 >  NB501016:282:H2V7GAFX2:4:11606:22358:2748/1‑113 (MQ=60)
                                                                                                       GGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAA                                          <  NB501016:282:H2V7GAFX2:3:21501:2103:12760/151‑1 (MQ=60)
                                                                                                               TGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGgggggggggggggggggggggg                                  >  NB501016:282:H2V7GAFX2:4:21605:6432:15301/1‑129 (MQ=60)
                                                                                                                GGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGggggggggggggggtgggggggg                                 >  NB501016:282:H2V7GAFX2:4:21508:13717:1988/1‑128 (MQ=60)
                                                                                                                   GGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATCTCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAG                              <  NB501016:282:H2V7GAFX2:1:21206:9030:11996/151‑1 (MQ=60)
                                                                                                                   GGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCcggggggataaacgaaaaagggcaaaa                               >  NB501016:282:H2V7GAFX2:4:21401:7696:17732/1‑123 (MQ=60)
                                                                                                                      CGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGC                             >  NB501016:282:H2V7GAFX2:1:21201:6862:17044/1‑149 (MQ=60)
                                                                                                                        ACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGGCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGgggggggggggggggggggggggcgggcgtc                         >  NB501016:282:H2V7GAFX2:3:11402:7106:12896/1‑120 (MQ=60)
                                                                                                                         CCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGgggggggggggggggggggggggggggggaac                        >  NB501016:282:H2V7GAFX2:2:11211:1773:13052/1‑119 (MQ=60)
                                                                                                                         CCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGgggggggggcggccggcgggggggacaacgcc                        >  NB501016:282:H2V7GAFX2:4:11509:5042:11422/1‑119 (MQ=60)
                                                                                                                          CATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGgggggggggggggggggggggcggggggggg                         >  NB501016:282:H2V7GAFX2:2:11206:12557:18658/1‑118 (MQ=60)
                                                                                                                           ATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGggggggggggggggggggggggggggggggggggggg                      >  NB501016:282:H2V7GAFX2:1:11310:3580:15188/1‑114 (MQ=60)
                                                                                                                           ATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCG                      >  NB501016:282:H2V7GAFX2:4:11505:11233:12131/1‑151 (MQ=60)
                                                                                                                                    TGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGggggggggggggggggggggggggggggtggggggggcgcggcggg             >  NB501016:282:H2V7GAFX2:2:11105:16364:6877/1‑105 (MQ=60)
                                                                                                                                          CAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGggggggggggggcgggggggggggaaaagcccccccccgtaatgtgatt       >  NB501016:282:H2V7GAFX2:4:21502:19360:13075/1‑102 (MQ=60)
                                                                                                                                              CGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGggggggggggggggggggggggggggggggggggcgggacggggcccgagcagg     >  NB501016:282:H2V7GAFX2:1:21311:8857:7666/1‑95 (MQ=60)
                                                                                                                                               gatTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGgggggggggggggcggggggggcggacaccaccccccccagagggagttatatg  >  NB501016:282:H2V7GAFX2:4:21610:2025:4400/4‑97 (MQ=60)
                                                                                                                                                  |                                                                                                                                                   
ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCT  >  NC_002947_CJ‑RC/335980‑336272

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