Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,754,949 T→A 5.0% E178V (GAG→GTG)  PP_2406 ← quinate/shikimate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,754,9490TA 5.0% 1436.2 / 36.8 417E178V (GAG→GTG) PP_2406quinate/shikimate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (216/178);  new base A (6/16);  total (222/195)
Fisher's exact test for biased strand distribution p-value = 1.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 422 total aligned reads displayed.

GATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGG  >  NC_002947/2754801‑2755098
                                                                                                                                                    |                                                                                                                                                     
gatgatCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTc                                                                                                                                                      <  2:8675200/150‑1 (MQ=255)
gatgatCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTc                                                                                                                                                      <  2:7825523/150‑1 (MQ=255)
gatgatCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTc                                                                                                                                                      >  1:2264318/1‑150 (MQ=255)
   gatCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCTACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCcgc                                                                                                                                                   <  1:1968032/150‑1 (MQ=255)
       tctGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGGAGCGCGCCg                                                                                                                                               <  1:5745447/150‑1 (MQ=255)
          gCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGaaa                                                                                                                                            <  1:8305410/150‑1 (MQ=255)
          gCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCACTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGaaa                                                                                                                                            <  2:4066884/150‑1 (MQ=255)
           ccacccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAAt                                                                                                                                           <  1:2127169/150‑1 (MQ=255)
             acccacCGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATgg                                                                                                                                         <  1:56590/150‑1 (MQ=255)
             acccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATgg                                                                                                                                         <  1:2872102/150‑1 (MQ=255)
               ccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGgt                                                                                                                                       <  2:188456/150‑1 (MQ=255)
                cacaGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGgtg                                                                                                                                      >  1:2478936/1‑150 (MQ=255)
                     gCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAg                                                                                                                                 <  2:3207835/150‑1 (MQ=255)
                      cGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAgg                                                                                                                                <  2:7973836/150‑1 (MQ=255)
                      cGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAgg                                                                                                                                >  1:6265963/1‑150 (MQ=255)
                          gATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGt                                                                                                                            >  2:7514003/1‑150 (MQ=255)
                             gCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTcca                                                                                                                         <  1:5317787/150‑1 (MQ=255)
                              caacaaCTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTccac                                                                                                                        <  2:8381683/150‑1 (MQ=255)
                              caacaaCTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTccac                                                                                                                        >  1:358918/1‑150 (MQ=255)
                               aacaacTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTccacc                                                                                                                       >  2:4506020/1‑150 (MQ=255)
                               aacaacTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTccacc                                                                                                                       >  1:2979384/1‑150 (MQ=255)
                               aacaacTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTccacc                                                                                                                       >  2:8202121/1‑150 (MQ=255)
                                acaacTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTccacca                                                                                                                      <  2:4289190/150‑1 (MQ=255)
                                 caacTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAACTGGGTGTTCAGGTTGTCCACCaa                                                                                                                     >  2:5453392/1‑150 (MQ=255)
                                 caacTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCaa                                                                                                                     >  2:6643551/1‑150 (MQ=255)
                                     tCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGcc                                                                                                                 >  1:6521516/1‑150 (MQ=255)
                                     tCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGcc                                                                                                                 >  2:768243/1‑150 (MQ=255)
                                      ccACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGTCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCt                                                                                                                >  2:541629/1‑150 (MQ=255)
                                      ccACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCt                                                                                                                >  1:1824685/1‑150 (MQ=255)
                                      ccACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCt                                                                                                                <  1:4133816/150‑1 (MQ=255)
                                       cACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTg                                                                                                               >  2:2934133/1‑150 (MQ=255)
                                            gCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCAc                                                                                                          <  2:3587756/150‑1 (MQ=255)
                                               gTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACggg                                                                                                       <  2:6199974/150‑1 (MQ=255)
                                                     tGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCAt                                                                                                 >  2:3007502/1‑150 (MQ=255)
                                                     tGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCAGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCAt                                                                                                 <  2:3910477/150‑1 (MQ=255)
                                                      gCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAATGCCTGCGCACGGGTCGCATc                                                                                                <  2:5316021/150‑1 (MQ=255)
                                                         cGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCAc                                                                                             <  1:1597028/150‑1 (MQ=255)
                                                          ggTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACt                                                                                            >  1:4487478/1‑150 (MQ=255)
                                                           gTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACtt                                                                                           >  1:6744738/1‑150 (MQ=255)
                                                            tAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTc                                                                                          <  2:6130609/150‑1 (MQ=255)
                                                               cTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGaa                                                                                       >  2:1523080/1‑150 (MQ=255)
                                                               cTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGaa                                                                                       <  1:4039341/150‑1 (MQ=255)
                                                               cTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGaa                                                                                       <  2:4129240/150‑1 (MQ=255)
                                                                 tGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACa                                                                                     >  1:3503068/1‑150 (MQ=255)
                                                                  ggCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAg                                                                                    >  1:5983101/1‑150 (MQ=255)
                                                                    ccATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCa                                                                                  >  1:3773960/1‑150 (MQ=255)
                                                                     cATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGTACGGGTCGCATCCACTTCGAACAGCAc                                                                                 >  1:8076906/1‑150 (MQ=255)
                                                                     cATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCAc                                                                                 <  1:7241399/150‑1 (MQ=255)
                                                                      aTGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCAcc                                                                                <  1:13592/150‑1 (MQ=255)
                                                                        gCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAg                                                                              >  2:3146547/1‑150 (MQ=255)
                                                                        gCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAg                                                                              >  2:5673807/1‑150 (MQ=255)
                                                                          ccACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGcc                                                                            >  2:8594643/1‑150 (MQ=255)
                                                                          ccACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGcc                                                                            <  2:2660881/150‑1 (MQ=255)
                                                                           cACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCa                                                                                  >  2:6841024/1‑143 (MQ=255)
                                                                           cACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCg                                                                           >  1:4559370/1‑150 (MQ=255)
                                                                                 cGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCGGCGGCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCTCGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTAc                                                                     <  1:4705961/150‑1 (MQ=255)
                                                                                 cGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTAc                                                                     <  1:311311/150‑1 (MQ=255)
                                                                                   tggtgTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTAccc                                                                   <  1:7719741/150‑1 (MQ=255)
                                                                                    ggtgTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTAcccc                                                                  <  1:8635593/150‑1 (MQ=255)
                                                                                      tgtTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCtt                                                                <  2:6049310/150‑1 (MQ=255)
                                                                                        ttGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCa                                                              <  2:2131492/150‑1 (MQ=255)
                                                                                           aCCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCAccc                                                           >  2:2395441/1‑150 (MQ=255)
                                                                                            ccAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCca                                                          >  2:2775869/1‑150 (MQ=255)
                                                                                             cAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCAc                                                                                             >  2:3716505/1‑114 (MQ=255)
                                                                                              aGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCcagc                                                        <  1:8149901/150‑1 (MQ=255)
                                                                                                cccGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCcagcag                                                      >  1:7118424/1‑150 (MQ=255)
                                                                                                 ccGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAgg                                                     >  1:2725250/1‑150 (MQ=255)
                                                                                                    tCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCa                                                  <  2:7242152/150‑1 (MQ=255)
                                                                                                      cGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCACATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATg                                                >  2:5878751/1‑150 (MQ=255)
                                                                                                       ggttgggCCTAGGCTGTGGGCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGATATGGGTGTTCTGGTTGTCCACGAACGCCTGCGCTGGGGTCGCATCGACTTCGATCCGGACGAGGGGTTCTTCCCCTTCACCCAGCAGGGCATgg                                               <  2:1841461/145‑1 (MQ=255)
                                                                                                       gCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATgg                                               <  1:4372047/150‑1 (MQ=255)
                                                                                                       gCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATgg                                               <  2:1266210/150‑1 (MQ=255)
                                                                                                       gCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATgg                                               <  1:5424007/150‑1 (MQ=255)
                                                                                                       gCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATgg                                               <  2:7019990/150‑1 (MQ=255)
                                                                                                         ttCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGc                                             >  2:114085/1‑150 (MQ=255)
                                                                                                          ctggCCAAGGCTGTGGCCAGGTCGGTTCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGcc                                            >  1:6908641/3‑150 (MQ=255)
                                                                                                            ggCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCAc                                          >  2:4936748/1‑150 (MQ=255)
                                                                                                            ggCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCAc                                          <  2:6471023/150‑1 (MQ=255)
                                                                                                             gCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACg                                         <  2:4176758/150‑1 (MQ=255)
                                                                                                              ccAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACgg                                        >  1:1565692/1‑150 (MQ=255)
                                                                                                               cAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGc                                       >  2:1626491/1‑150 (MQ=255)
                                                                                                                  ggCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGa                                    <  1:481942/150‑1 (MQ=255)
                                                                                                                  ggCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGa                                    <  2:7627454/150‑1 (MQ=255)
                                                                                                                     tgtgGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGcc                                 >  1:4622318/1‑150 (MQ=255)
                                                                                                                          ccAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCAcc                                                                                <  2:1012203/98‑1 (MQ=255)
                                                                                                                          ccAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCAcc                                                                                >  1:1012203/1‑98 (MQ=255)
                                                                                                                          ccAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAgcgc                            >  2:1365123/1‑150 (MQ=255)
                                                                                                                             ggTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAATAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCAc                         <  1:3194125/150‑1 (MQ=255)
                                                                                                                                  gTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCa                    <  2:4778007/150‑1 (MQ=255)
                                                                                                                                  gTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCa                    <  1:3380129/150‑1 (MQ=255)
                                                                                                                                   tGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCAc                   <  1:5946623/150‑1 (MQ=255)
                                                                                                                                     ccAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCAccc                 >  1:8238107/1‑150 (MQ=255)
                                                                                                                                      cAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCa                >  2:4585186/1‑150 (MQ=255)
                                                                                                                                       aaGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCAt               >  1:6113070/1‑150 (MQ=255)
                                                                                                                                       aaGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCAt               >  1:5103805/1‑150 (MQ=255)
                                                                                                                                        aGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATc              <  1:7514003/150‑1 (MQ=255)
                                                                                                                                           aCGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGc           >  1:1109183/1‑150 (MQ=255)
                                                                                                                                            cGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGca          >  1:7380523/1‑150 (MQ=255)
                                                                                                                                                  gcTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCt    <  2:2900693/150‑1 (MQ=255)
                                                                                                                                                    tCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCTAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTgg  >  2:6376383/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGG  >  NC_002947/2754801‑2755098

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: