Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,862,776 G→T 7.6% L81M (CTG→ATG)  ligA ← DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,862,7760GT 7.6% 802.3 / 40.1 238L81M (CTG→ATG) ligADNA ligase
Reads supporting (aligned to +/- strand):  ref base G (124/95);  new base T (11/7);  total (135/103)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 239 total aligned reads displayed.

ACCAGCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGC  >  NC_002947/4862628‑4862923
                                                                                                                                                    |                                                                                                                                                   
  cAGCTGGCCGTCGCGGTACAGCGGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCt                                                                                                                                                  >  2:4148809/1‑150 (MQ=255)
  cAGCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCt                                                                                                                                                  >  2:2707384/1‑150 (MQ=255)
  cAGCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCt                                                                                                                                                  >  2:2707580/1‑150 (MQ=255)
    gCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCCTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTTCCAAGCGACGGCCGAACTCACGCAGTTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGa                                                                                                                                                <  1:2661760/150‑1 (MQ=255)
         ccGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATg                                                                                                                                           <  1:999751/150‑1 (MQ=255)
          cGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATgg                                                                                                                                          <  1:1196041/150‑1 (MQ=255)
           gTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATggg                                                                                                                                         <  1:815854/150‑1 (MQ=255)
               cgGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACt                                                                                                                                     >  2:3919648/1‑150 (MQ=255)
                ggTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCATGCTCAGCCTGGGCaacg                                                                                                                                    <  2:3238828/150‑4 (MQ=255)
                  tACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCg                                                                                                                                  <  2:1941521/150‑1 (MQ=255)
                     agcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTgg                                                                                                                               <  2:2737965/150‑1 (MQ=255)
                            tCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCt                                                                                                                        >  1:2273346/1‑150 (MQ=255)
                            tCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCt                                                                                                                        >  1:3369125/1‑150 (MQ=255)
                              aCGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTgg                                                                                                                      >  1:3584675/1‑150 (MQ=255)
                              aCGGCCAGGCCATCGAGCTTGGGTACGCAGCAGTAGTCGACCGCGCCCGGCTGATCGAGGCCATCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTgg                                                                                                                      >  2:2411867/1‑150 (MQ=255)
                               cGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGc                                                                                                                     <  2:3590920/150‑1 (MQ=255)
                               cGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGc                                                                                                                     <  2:2592519/150‑1 (MQ=255)
                                ggccaggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCt                                                                                                                    <  1:2313845/150‑1 (MQ=255)
                                 gccaggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTg                                                                                                                   >  2:3836666/1‑150 (MQ=255)
                                 gccaggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTg                                                                                                                   <  2:2322278/150‑1 (MQ=255)
                                    aggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAc                                                                                                                <  1:2749091/150‑1 (MQ=255)
                                    aggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAc                                                                                                                <  2:1802893/150‑1 (MQ=255)
                                     ggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAcg                                                                                                               >  1:866529/1‑150 (MQ=255)
                                            gAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCaa                                                                                                        <  1:2568684/150‑1 (MQ=255)
                                            gAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCaa                                                                                                        <  1:357619/150‑1 (MQ=255)
                                                ttGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCa                                                                                                    >  1:321470/1‑150 (MQ=255)
                                                  gggTTCGCAGCTGTATTCGACCGCGCCCGGCTGAACGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGCCTGAACGCCGCCAACGCAGc                                                                                                  >  2:1128540/1‑150 (MQ=255)
                                                  gggTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGc                                                                                                  >  2:3405790/1‑150 (MQ=255)
                                                     ttCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCGGCTcc                                                                                               >  2:703838/1‑150 (MQ=255)
                                                     ttCGCAGCTGTAGACGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGTGTACCTGGCTGAACGCCGCCAAGGCAGCTcc                                                                                               >  1:2546635/1‑150 (MQ=255)
                                                       cGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccac                                                                                             <  1:629751/150‑1 (MQ=255)
                                                        gCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccacc                                                                                            >  1:659893/1‑150 (MQ=255)
                                                        gCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccacc                                                                                            >  2:1563540/1‑150 (MQ=255)
                                                        gCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccacc                                                                                            >  2:4356727/1‑150 (MQ=255)
                                                             tGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCacac                                                                                       >  1:3700266/1‑150 (MQ=255)
                                                              gTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACg                                                                                      >  1:1734672/1‑150 (MQ=255)
                                                              gTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACg                                                                                      >  1:4456425/1‑150 (MQ=255)
                                                                 gTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtg                                                                                   <  1:4031130/150‑1 (MQ=255)
                                                                 gTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtg                                                                                   <  2:3626876/150‑1 (MQ=255)
                                                                  tCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgg                                                                                  >  2:4392846/1‑150 (MQ=255)
                                                                  tCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgg                                                                                  >  1:2624951/1‑150 (MQ=255)
                                                                     aCCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtg                                                                               >  2:3871223/1‑150 (MQ=255)
                                                                     aCCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtg                                                                               >  2:2234760/1‑150 (MQ=255)
                                                                        gcgcCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGc                                                                            <  2:4109629/150‑1 (MQ=255)
                                                                            ccGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGt                                                                        >  2:2159733/1‑150 (MQ=255)
                                                                            ccGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGt                                                                        >  1:4398201/1‑150 (MQ=255)
                                                                                     cGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTa                                                               >  2:583539/1‑150 (MQ=255)
                                                                                        ggCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCa                                                            >  1:2006196/1‑150 (MQ=255)
                                                                                         gCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAg                                                           <  1:670973/150‑1 (MQ=255)
                                                                                         gCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAg                                                           <  1:2047295/150‑1 (MQ=255)
                                                                                           cTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGt                                                         >  1:17971/1‑150 (MQ=255)
                                                                                            ttCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTg                                                        >  2:3152312/1‑150 (MQ=255)
                                                                                                aCCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCggg                                                    >  2:1665887/1‑150 (MQ=255)
                                                                                                   aCGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGgtgt                                                 >  1:4001332/1‑150 (MQ=255)
                                                                                                     gcgACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTc                                               >  2:747144/1‑150 (MQ=255)
                                                                                                      cgACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCa                                              <  2:4245134/150‑1 (MQ=255)
                                                                                                       gACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAg                                             >  1:4165653/1‑150 (MQ=255)
                                                                                                             cGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                       >  2:2407729/1‑150 (MQ=255)
                                                                                                             cGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                       >  2:1535358/1‑150 (MQ=255)
                                                                                                             cGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                       >  1:2185771/1‑150 (MQ=255)
                                                                                                               aaCTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAgc                                     <  1:1379797/150‑1 (MQ=255)
                                                                                                                  tCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGcc                                  >  2:918284/1‑150 (MQ=255)
                                                                                                                  tCACGCAGGTCGGCCTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGcc                                  >  2:2806754/1‑150 (MQ=255)
                                                                                                                      gCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGa                              <  1:4120175/150‑1 (MQ=255)
                                                                                                                      gCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGa                              >  1:961227/1‑150 (MQ=255)
                                                                                                                      gCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGa                              >  2:4266876/1‑150 (MQ=255)
                                                                                                                          gTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTc                          >  1:1433330/1‑150 (MQ=255)
                                                                                                                            cGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGt                        <  1:3152312/150‑1 (MQ=255)
                                                                                                                             ggCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTGCCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGtt                       <  1:3528890/150‑1 (MQ=255)
                                                                                                                             ggCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGtt                       <  1:1646363/150‑1 (MQ=255)
                                                                                                                              gCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTg                      <  2:378403/150‑1 (MQ=255)
                                                                                                                               cttcttCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGa                     >  2:4017123/1‑150 (MQ=255)
                                                                                                                                ttcttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGaa                    >  2:2551459/1‑150 (MQ=255)
                                                                                                                                ttcttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGaa                    >  2:133365/1‑150 (MQ=255)
                                                                                                                                 tcttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAAc                   <  2:2850900/150‑1 (MQ=255)
                                                                                                                                  cttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACa                  >  1:1959521/1‑150 (MQ=255)
                                                                                                                                   ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAg                 >  1:2108269/1‑150 (MQ=255)
                                                                                                                                    tcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAgg                >  2:2911041/1‑150 (MQ=255)
                                                                                                                                    tcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAgg                >  1:2487394/1‑150 (MQ=255)
                                                                                                                                      gAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCg              >  2:3954205/1‑150 (MQ=255)
                                                                                                                                       aaGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCgg             <  1:4343181/150‑1 (MQ=255)
                                                                                                                                       aaGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCgg             >  1:3632788/1‑150 (MQ=255)
                                                                                                                                         cgaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTc           >  2:324190/6‑150 (MQ=255)
                                                                                                                                          gaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCg          >  1:1034677/5‑150 (MQ=255)
                                                                                                                                           cGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGt         <  2:191472/150‑1 (MQ=255)
                                                                                                                                              tGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                >  1:3087403/1‑124 (MQ=255)
                                                                                                                                              tGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                <  2:3087403/124‑1 (MQ=255)
                                                                                                                                              tGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTAtt      >  1:2805501/1‑150 (MQ=255)
                                                                                                                                               gCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTc     >  2:1114599/1‑150 (MQ=255)
                                                                                                                                                cccATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTc           >  2:4485045/1‑143 (MQ=255)
                                                                                                                                                cccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCg    >  1:2725101/1‑150 (MQ=255)
                                                                                                                                                cccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCg    >  1:475224/1‑150 (MQ=255)
                                                                                                                                                cccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCg    >  2:184227/1‑150 (MQ=255)
                                                                                                                                                cccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCg    >  1:1120895/1‑150 (MQ=255)
                                                                                                                                                 ccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCgg   >  1:414196/1‑150 (MQ=255)
                                                                                                                                                 ccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCgg   >  1:566101/1‑150 (MQ=255)
                                                                                                                                                 ccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCgg   >  1:1872116/1‑150 (MQ=255)
                                                                                                                                                 ccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCgg   >  2:2019946/1‑150 (MQ=255)
                                                                                                                                                  cAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGc  >  1:330429/1‑150 (MQ=255)
                                                                                                                                                  cAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGc  >  2:2514909/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCAGCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGC  >  NC_002947/4862628‑4862923

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