Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,577,913 C→T 5.8% intergenic (+78/+31) PP_4057 → / ← glgB membrane protein/1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,577,9130CT 5.8% 400.0 / 9.9 137intergenic (+78/+31)PP_4057/glgBmembrane protein/1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (60/69);  new base T (5/3);  total (65/72)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 140 total aligned reads displayed.

CAGTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACC  >  NC_002947/4577769‑4578062
                                                                                                                                                |                                                                                                                                                      
cagTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCcg                                                                                                                                                  >  1:26760/1‑150 (MQ=255)
     atcgctggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccccttcgcgggtaaACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGaaa                                                                                                                                             >  2:2350860‑M1/88‑150 (MQ=255)
      tgtgGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAg                                                                                                                                            >  2:403052‑M2/1‑149 (MQ=255)
       gtgGATGACAATGCCCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCGCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgg                                                                                                                                           >  1:171219‑M2/1‑148 (MQ=255)
        tgGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAggg                                                                                                                                          >  2:2195990‑M2/1‑147 (MQ=255)
         ggATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggc                                                                                                                                         >  2:1289349‑M2/1‑146 (MQ=255)
          gATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggcc                                                                                                                                        <  2:403175‑M2/150‑6 (MQ=255)
           ggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccccttcgcgggtaaACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCa                                                                                                                                       >  1:896945‑M1/82‑150 (MQ=255)
           aTGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCa                                                                                                                                       >  1:1886934/1‑150 (MQ=255)
             gACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGt                                                                                                                                     <  1:1034324/150‑1 (MQ=255)
               cAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgg                                                                                                                                           <  1:6874‑M2/142‑3 (MQ=255)
                     ctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtaccggccccttcgcgggtaaACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCaa                                                                                                                             >  2:300221‑M1/72‑150 (MQ=255)
                     ccTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCaa                                                                                                                             <  1:2403812/150‑1 (MQ=255)
                     ccTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCaa                                                                                                                             <  2:1125273/150‑1 (MQ=255)
                      cTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAAc                                                                                                                            <  1:1711396/150‑1 (MQ=255)
                           cGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATc                                                                                                                       <  2:901527/150‑1 (MQ=255)
                            gggTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGAACACCAAACCgctac                                                                                                                      >  2:2352496/1‑146 (MQ=255)
                            gggTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTAGAGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCa                                                                                                                      >  1:1013572/1‑150 (MQ=255)
                             ggTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAg                                                                                                                     <  1:2352496/150‑1 (MQ=255)
                               ttGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggc                                                                                                                                         >  1:215604‑M2/1‑124 (MQ=255)
                               ttGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggc                                                                                                                                         <  2:215604‑M2/128‑5 (MQ=255)
                                      gCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAg                                                                                                                     >  1:216130/1‑141 (MQ=255)
                                      gCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGc                                                                                                                   >  2:588676/1‑143 (MQ=255)
                                      gCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGc                                                                                                                   >  2:588677/1‑143 (MQ=255)
                                      gCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGc                                                                                                                   <  1:588677/143‑1 (MQ=255)
                                      gCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGNCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGc                                                                                                                   <  1:588676/143‑1 (MQ=255)
                                          ccTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGa                                                                                                        <  2:1455258/150‑1 (MQ=255)
                                            ttGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATc                                                                                                      <  2:1197359/150‑1 (MQ=255)
                                               gggTGCGGATGCGGTGGTATGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAac                                                                                                   <  2:182535/150‑1 (MQ=255)
                                               gggTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAac                                                                                                   <  1:403175/150‑1 (MQ=255)
                                                 gTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGtct                                                                                                                                                       <  2:832516/97‑1 (MQ=255)
                                                      gATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCa                                                                                            >  2:1076141/1‑150 (MQ=255)
                                                           ggtggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggc                                                                                       >  2:341899/1‑150 (MQ=255)
                                                            gtggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCa                                                                                                           >  2:1769666/1‑129 (MQ=255)
                                                            gtggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCa                                                                                                           <  1:1769666/129‑1 (MQ=255)
                                                             tggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggcgg                                                                                     >  2:1607767/1‑150 (MQ=255)
                                                               gtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCa                                                                                                                      >  2:1986518/1‑115 (MQ=255)
                                                               gtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCa                                                                                                                      <  1:1986518/115‑1 (MQ=255)
                                                                tAGGCAGGGCCGGCCCCTTCGCGGGTCAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAg                                                                                  <  2:2317993/150‑1 (MQ=255)
                                                                tAGGCAGGGCCGGCCCCTTCGCGGGTAACCCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAg                                                                                  <  2:1886934/150‑1 (MQ=255)
                                                                tAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGTCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAg                                                                                  >  2:1330143/1‑150 (MQ=255)
                                                                   gCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGAGCCTGCGTAGGAACGGGTCTCCCGCGAAAGTGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGtt                                                                               >  2:270097/1‑150 (MQ=255)
                                                                    cAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTc                                                                              >  1:1532798/1‑150 (MQ=255)
                                                                      gggCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                            <  1:89088/150‑1 (MQ=255)
                                                                       ggccggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggc                                                                                       <  1:1312414/138‑1 (MQ=255)
                                                                       ggccggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAgg                                                                           <  1:2210722/150‑1 (MQ=255)
                                                                          cggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcat                                                                        <  2:2395323‑M2/150‑70 (MQ=255)
                                                                            gccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcatcc                                                                      <  2:2138193‑M2/150‑72 (MQ=255)
                                                                            gccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGCTGAACACACCCAACGTCGGCAGGGTCAGGGCCAg                                                                      >  1:748970/1‑150 (MQ=255)
                                                                            gccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAg                                                                      >  1:489317/1‑150 (MQ=255)
                                                                             ccccTCCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggccagccctgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcatcca                                                                     <  2:221586‑M2/150‑73 (MQ=255)
                                                                             caccTTCGCGGGGAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGc                                                                     >  2:958510/3‑150 (MQ=255)
                                                                               ccTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAGgggccggccatgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcatccaca                                                                   <  1:1446157‑M2/150‑75 (MQ=255)
                                                                               ccTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCa                                                                                                                                       >  2:278373/1‑82 (MQ=255)
                                                                                cTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGAt                                                                  <  2:489317/150‑1 (MQ=255)
                                                                                 tttcgcgggagacccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGAt                                                                  >  2:2350860‑M2/64‑150 (MQ=255)
                                                                                   cgcgGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAGgggccggccatgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcatccacatcac                                                               >  2:1446157‑M2/1‑72 (MQ=255)
                                                                                      gggTCAACCCGCGCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGc                                                            <  1:196141/150‑1 (MQ=255)
                                                                                       ggTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCt                                                           <  1:1076141/150‑1 (MQ=255)
                                                                                        gTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTg                                                          >  2:1964073/1‑150 (MQ=255)
                                                                                            ACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCg                                                      <  1:1443417‑M1/150‑1 (MQ=255)
                                                                                             cccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggaacgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGAAGGTTCAGGGCCAGCGATAATGGCTGCCCGt                                                     >  1:1078394‑M2/52‑150 (MQ=255)
                                                                                              ccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtg                                                    >  1:1277454‑M2/51‑150 (MQ=255)
                                                                                              ccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtg                                                    >  1:1278386‑M2/51‑150 (MQ=255)
                                                                                              ccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtg                                                    >  1:1100230‑M2/51‑150 (MQ=255)
                                                                                              ccGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaac                                                                                           <  2:823925‑M2/111‑51 (MQ=255)
                                                                                              ccGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaac                                                                                           >  1:823925‑M2/1‑61 (MQ=255)
                                                                                              CCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTg                                                    <  1:994128‑M1/150‑1 (MQ=255)
                                                                                              CCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTg                                                    <  1:1802390‑M1/150‑1 (MQ=255)
                                                                                              CCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTg                                                    <  1:222359‑M1/150‑1 (MQ=255)
                                                                                               CGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAGCGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTgg                                                   >  2:1013572‑M1/1‑150 (MQ=255)
                                                                                                    cTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                            >  2:1073798/1‑120 (MQ=255)
                                                                                                     tacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGcc                                             <  2:1272607‑M2/107‑1 (MQ=255)
                                                                                                     tacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGcc                                             <  2:1277454‑M2/107‑1 (MQ=255)
                                                                                                         cGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACtt                                         <  1:761345/150‑1 (MQ=255)
                                                                                                          gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcgg                                                                                        <  1:255380/102‑1 (MQ=255)
                                                                                                          gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcgg                                                                                        >  2:255380/1‑102 (MQ=255)
                                                                                                          gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAgg                                                                                 >  1:722219/1‑109 (MQ=255)
                                                                                                          gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAgg                                                                                 <  2:722219/109‑1 (MQ=255)
                                                                                                          gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATgg                                                             >  1:1268656/1‑129 (MQ=255)
                                                                                                          aggcaccgcacccggcggtgacccgcgtaggaacgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGa                                                                                                        <  1:1349054‑M2/48‑1 (MQ=255)
                                                                                                          aggcaccgcacccggcggtgacccgcgtaggaacgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGa                                                                                                        >  2:1349054‑M2/39‑86 (MQ=255)
                                                                                                                ccgccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACt                                  >  1:833199/1‑150 (MQ=255)
                                                                                                                 cgccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTg                                 <  1:341899/150‑1 (MQ=255)
                                                                                                                 cgccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTg                                 <  2:1597553/150‑1 (MQ=255)
                                                                                                                   ccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTCCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTc                                        >  2:1311574/1‑141 (MQ=255)
                                                                                                                   ccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTCCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTc                                        <  1:1311574/141‑1 (MQ=255)
                                                                                                                       gtgcggtgcCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCAcc                           >  1:1136665/1‑150 (MQ=255)
                                                                                                                                cTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCgaga                  >  2:2269025/1‑150 (MQ=255)
                                                                                                                                  gCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGAtt                >  2:1358915/1‑150 (MQ=255)
                                                                                                                                     tAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGcc             <  2:1532798/150‑1 (MQ=255)
                                                                                                                                        gAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcatccacatcactgctgtaatcggcgc                                                >  2:2336547‑M2/1‑19 (MQ=255)
                                                                                                                                        gAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcatccacatcactgctgtaatcggcgc                                                <  1:2336547‑M2/112‑94 (MQ=255)
                                                                                                                                        gAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGAc          >  1:424620/1‑150 (MQ=255)
                                                                                                                                         aaCGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACa         >  2:1276669/1‑150 (MQ=255)
                                                                                                                                          aCGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                            >  2:333207/1‑82 (MQ=255)
                                                                                                                                          aCGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                            <  1:333207/82‑1 (MQ=255)
                                                                                                                                              gtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGa    <  1:1064188‑M2/148‑1 (MQ=255)
                                                                                                                                              gtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGAAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGa    <  2:1078394‑M2/148‑1 (MQ=255)
                                                                                                                                                ctcCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGAcc  <  1:2269025/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
CAGTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACC  >  NC_002947/4577769‑4578062

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