Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,862,780 G→C 7.6% L79L (CTC→CTG)  ligA ← DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,862,7800GC 7.6% 799.8 / 43.2 237L79L (CTC→CTG) ligADNA ligase
Reads supporting (aligned to +/- strand):  ref base G (119/100);  new base C (11/7);  total (130/107)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 238 total aligned reads displayed.

AGCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCAG  >  NC_002947/4862631‑4862928
                                                                                                                                                     |                                                                                                                                                    
aGCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTg                                                                                                                                                      >  2:1858723/1‑150 (MQ=255)
       cGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATgg                                                                                                                                               <  1:1196041/150‑1 (MQ=255)
       cGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATgg                                                                                                                                               <  1:86523/150‑1 (MQ=255)
        gTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATggg                                                                                                                                              <  1:2862469/150‑1 (MQ=255)
         tCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGc                                                                                                                                             >  1:3726026/1‑150 (MQ=255)
            cgGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACt                                                                                                                                          <  2:364023/150‑1 (MQ=255)
            cgGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACt                                                                                                                                          >  2:3919648/1‑150 (MQ=255)
             ggTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCATGCTCAGCCTGGGCaacg                                                                                                                                         <  2:3238828/150‑4 (MQ=255)
              gTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCATGCTCAGCCTGGGCaacgc                                                                                                                                        <  1:3399472/150‑5 (MQ=255)
               tACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGTTCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCg                                                                                                                                       >  1:271056/1‑150 (MQ=255)
               tACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCg                                                                                                                                       <  2:1941521/150‑1 (MQ=255)
                aCAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGt                                                                                                                                      <  1:1653389/150‑1 (MQ=255)
                 cagcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTg                                                                                                                                     >  1:191472/1‑150 (MQ=255)
                  agcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTgg                                                                                                                                    <  2:2845401/150‑1 (MQ=255)
                   gcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGc                                                                                                                                   <  2:508436/150‑1 (MQ=255)
                         tCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCt                                                                                                                             >  1:2406593/1‑150 (MQ=255)
                            cGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGc                                                                                                                          <  2:3590920/150‑1 (MQ=255)
                             ggccaggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCt                                                                                                                         <  1:2313845/150‑1 (MQ=255)
                              gccaggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTg                                                                                                                        >  2:3836666/1‑150 (MQ=255)
                                 aggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAc                                                                                                                     <  1:2749091/150‑1 (MQ=255)
                                      aTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAcgccgc                                                                                                                >  2:1808859/1‑150 (MQ=255)
                                       tCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACgccgcc                                                                                                               >  1:2247193/1‑150 (MQ=255)
                                       tCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACgccgcc                                                                                                               >  1:3583130/1‑150 (MQ=255)
                                        cGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCa                                                                                                              <  1:4060374/150‑1 (MQ=255)
                                         gAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCaa                                                                                                             <  1:2568587/150‑1 (MQ=255)
                                         gAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCaa                                                                                                             >  1:2150304/1‑150 (MQ=255)
                                              tGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAg                                                                                                        <  1:3594020/150‑1 (MQ=255)
                                               gggTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGc                                                                                                       >  2:3405790/1‑150 (MQ=255)
                                                ggTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCt                                                                                                      <  1:436258/150‑1 (MQ=255)
                                                 gTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTc                                                                                                     >  1:559989/1‑150 (MQ=255)
                                                 gTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTc                                                                                                     >  2:286173/1‑150 (MQ=255)
                                                  ttCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTcc                                                                                                    >  2:704021/1‑150 (MQ=255)
                                                    cGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccac                                                                                                  <  1:629811/150‑1 (MQ=255)
                                                     gCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccacc                                                                                                 >  1:4450049/1‑150 (MQ=255)
                                                       aGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccaccca                                                                                               >  1:3308749/1‑150 (MQ=255)
                                                           gTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACg                                                                                           >  1:1734672/1‑150 (MQ=255)
                                                             aGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGcg                                                                                         >  1:666132/1‑149 (MQ=255)
                                                              gTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtg                                                                                        <  2:2614161/150‑1 (MQ=255)
                                                                 gACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggt                                                                                     <  1:982396/150‑1 (MQ=255)
                                                                 gACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggt                                                                                     <  1:3266823/150‑1 (MQ=255)
                                                                  aCCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtg                                                                                    >  2:2234760/1‑150 (MQ=255)
                                                                  aCCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtg                                                                                    >  2:3871223/1‑150 (MQ=255)
                                                                     gcgcCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGc                                                                                 <  2:4109629/150‑1 (MQ=255)
                                                                         ccGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGt                                                                             >  2:2159733/1‑150 (MQ=255)
                                                                             cTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCggg                                                                         <  1:4125498/150‑1 (MQ=255)
                                                                             cTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCggg                                                                         <  1:3566549/150‑1 (MQ=255)
                                                                                 tCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGtt                                                                     >  2:1395349/1‑150 (MQ=255)
                                                                                   gAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTAc                                                                   <  2:4374145/150‑1 (MQ=255)
                                                                                     ggCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCa                                                                 >  1:2006196/1‑150 (MQ=255)
                                                                                     ggCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCa                                                                 <  2:3288720/150‑1 (MQ=255)
                                                                                     ggCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCa                                                                 >  1:4274860/1‑150 (MQ=255)
                                                                                      gCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAg                                                                <  2:4171608/150‑1 (MQ=255)
                                                                                       ccTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAgg                                                               >  2:4283275/1‑150 (MQ=255)
                                                                                        cTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGt                                                              >  1:17971/1‑150 (MQ=255)
                                                                                         ttCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTg                                                             >  2:3152312/1‑150 (MQ=255)
                                                                                           cTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCg                                                           >  1:4028937/1‑150 (MQ=255)
                                                                                             aCCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCggg                                                         >  2:1772240/1‑150 (MQ=255)
                                                                                              ccccgcgCCGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGgt                                                        <  2:2762827/147‑1 (MQ=255)
                                                                                                aCGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGgtgt                                                      >  1:2632423/1‑150 (MQ=255)
                                                                                                  gcgACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTc                                                    <  2:1641820/150‑1 (MQ=255)
                                                                                                   cgACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCa                                                   <  2:4245156/150‑1 (MQ=255)
                                                                                                   cgACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCa                                                   <  2:3996669/150‑1 (MQ=255)
                                                                                                          cGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                            >  2:1535358/1‑150 (MQ=255)
                                                                                                          cGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                            >  2:2407729/1‑150 (MQ=255)
                                                                                                            aaCTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAgc                                          <  1:1379797/150‑1 (MQ=255)
                                                                                                               tCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGcc                                       >  2:918284/1‑150 (MQ=255)
                                                                                                               tCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGcc                                       >  2:2462879/1‑150 (MQ=255)
                                                                                                                 aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                     <  1:4142451/150‑1 (MQ=255)
                                                                                                                       gTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTc                               >  1:1433330/1‑150 (MQ=255)
                                                                                                                       gTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTc                               >  1:2490268/1‑150 (MQ=255)
                                                                                                                        tCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCg                              >  1:2090999/1‑150 (MQ=255)
                                                                                                                         cGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGt                             <  1:3152312/150‑1 (MQ=255)
                                                                                                                         cGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGt                             >  2:1126916/1‑150 (MQ=255)
                                                                                                                           gCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTg                           <  1:2573158/150‑1 (MQ=255)
                                                                                                                           gCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTg                           <  2:378403/150‑1 (MQ=255)
                                                                                                                            cttcttCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGa                          >  1:946897/1‑150 (MQ=255)
                                                                                                                            cttcttCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGa                          >  2:4017123/1‑150 (MQ=255)
                                                                                                                             ttcttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGac                         <  1:574239/150‑2 (MQ=255)
                                                                                                                             ttcttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGaa                         >  2:2551459/1‑150 (MQ=255)
                                                                                                                              tcttcGAAGGCGTTGCCCTGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCTAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAAc                        <  2:2357707/150‑1 (MQ=255)
                                                                                                                              tcttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAAc                        <  2:2726277/150‑1 (MQ=255)
                                                                                                                                ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAg                      >  1:1843811/1‑150 (MQ=255)
                                                                                                                                ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAg                      >  1:701758/1‑150 (MQ=255)
                                                                                                                                ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAg                      <  2:4450049/150‑1 (MQ=255)
                                                                                                                                 tcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAgg                     >  1:3649116/1‑150 (MQ=255)
                                                                                                                                    cacgaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCgg                  >  2:1281051/8‑150 (MQ=255)
                                                                                                                                     acgaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGt                 >  2:1073686/7‑150 (MQ=255)
                                                                                                                                       gaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCg               >  1:1034677/5‑150 (MQ=255)
                                                                                                                                       gCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCg               >  2:626884/1‑150 (MQ=255)
                                                                                                                                           tGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                            <  1:2036988/117‑1 (MQ=255)
                                                                                                                                           tGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                            >  2:2036988/1‑117 (MQ=255)
                                                                                                                                           tGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                     >  1:3087403/1‑124 (MQ=255)
                                                                                                                                           tGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                     <  2:3087403/124‑1 (MQ=255)
                                                                                                                                             cccATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTc                <  1:4484720/143‑1 (MQ=255)
                                                                                                                                             cccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCg         >  1:1120895/1‑150 (MQ=255)
                                                                                                                                               cAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGc       >  2:2514909/1‑150 (MQ=255)
                                                                                                                                                  gCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATc    <  2:3028334/150‑1 (MQ=255)
                                                                                                                                                   cTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTTAACAGGCGGTCGTATTCGGCATCa   <  2:270950/150‑1 (MQ=255)
                                                                                                                                                    tGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCAg  >  1:1537596/1‑150 (MQ=255)
                                                                                                                                                    tGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCAg  <  2:813668/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCAG  >  NC_002947/4862631‑4862928

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: ‑